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FOREWORD 


This  interim  report  was  prepared  by  the  Advanced  Engineering  and 
Technology  Programs  Department,  Aircraft  Engine  Group  of  the  General 
Electric  Company,  Evendale,  Ohio  under  the  joint  sponsorship  of  the  Air 
Force  Aero  Propulsion  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio 
and  the  Department  of  Transportation,  Washington,  D.  C.  under  Contract 
F33615-73-C-2031.  The  inclusive  dates  of  this  interim  program  activity 
were  December  1972  through  December  1973.  The  work  described  is  part  of 
a  program  to  define  and  control  the  noise  emission  of  aircraft  propulsion 
systems. 

Mr.  Paul  A.  Shahady  of  the  Air  Force  Aero  Propulsion  Laboratory 
(AFAPL/TBC)  was  the  Project  Engineer.  The  program  is  being  conducted 
under  Project  3066,  Task  14. 

The  principal  technical  contributors  to  this  report  are:  R.  Mani, 

C.  Merkle,  P.  Scott,  P.  Mossey,  R.  Kantola,  K.  R.  Bilwakesh,  J.  Wang, 

J.  F.  Brausch,  H.  S.  Ribner  and  S.  P.  Pao. 

This  document  covers  interim  results  of  theoretical  and  experimental 
investigations  necessary  to  reveal  the  basic  mechanisms  of  supersonic 
exhaust  noise  typical  of  present  and  future  military  and  commercial  super¬ 
sonic  aircraft  propulsion  systems.  The  program  was  conducted  by  the 
General  Electric  Company  under  the  direction  of  Dr.  Paul  R.  Knott.  The 
report  was  submitted  by  the  author(s)  on  15  February  1974. 

Publication  of  this  report  does  not  constitute  Air  Force  approval 
of  the  report's  findings  or  conclusions.  It  is  published  only  for  the 
exchange  and  stimulation  of  ideas. 

/?■  \ 

PAUL  A.  SHAHADY/" 

Project  Engineer 

FOR  TIE  COMMANDER 


'ERNEST  C.  SIMEON 
Director,  Turbine  Engine  Division 


ABSTRACT 


This  interim  report  summarizes  the  major  theoretical  and  experimental  task 
efforts  performed  at  General  Electric  during  the  first  year  of  a  two  yuar  follow- 
on  program  sponsored  jointly  by  the  Air  Force  and  the  Department  of  Transpor¬ 
tation  on  Supersonic  Jet  Exhaust  Noise.  The  overall  objective  of  the  Program 
is  to  develop  the  technology  to  significantly  reduce  supersonic  aircraft  pro¬ 
pulsion  system  noise  with  minimum  associated  performance  and  weight  penalties. 

To  reach  the  objectives  of  this  program  a  varied  and  comprehensive  research 
program  is  being  carried  out  to  develop  the  basic  rheory  and  experimental  methods 
for  understanding  and  quantizing  the  acoustic  characteristics  of  simple  super¬ 
sonic  jets  for  a  range  of  velocities  and  temperatures  typical  of  present  and 
future  military  and  commercial  supersonic  aircraft  propulsion  systems.  A 
comprehensive  aero-acoustic  model,  which  had  its  origin  in  the  first  phase  of 
effort,  relating  the  local  fluid  dynamic  properties  of  the  jet  to  the  jet 
acoustic  nature  is  continuing  in  its  development.  The  main  thrust  in  the  de¬ 
velopment  of  this  aero-acouatic  model  has  been  1)  to  refine  the  calculation  of 
the  aerodynamic  input  by  including  the  effects  of  shock  waves  on  heated  jet 
flow  properties,  and  2)  to  more  clearly  delineate  the  acoustic  model  to  account 
for  heretofore  parodoxlcal  jet  acoustic  observations.  It  is  shown  that  the 
shock  structure  and  turbulent  fluid  dynamic  properties  for  a  heated  supersonic 
exhaust  jet  can  be  accurately  predicted.  Further,  comprehensive  analyses  are 
presented  which  explicitly  account  for  the  non-classical  density  dependence  of 
jet  noise,  the  influence  of  mean  flow  on  the  real  radiative  efficiency  of 
moving  sources  and  some  fundamental  theoretical  questions  regarding  the  inter¬ 
play  between  convection  and  refraction  in  jet  acoustic  propagation.  Prelim¬ 
inary  results  of  parametric  far-field  and  near-field  acoustic  experiments  are 
presented.  Results  are  displayed  in  such  a  way  as  to  illustrate  many  of  the 
salient  features  of  the  velocity,  density  and  spectral  dependency  of  heated 
supersonic  jets.  Major  advances  in  developing  General  Electric's  Laser  Veloci- 
meter  for  performing  turbulence  spectra  measurements  in  heated  high  velocity 
jets  is  discussed.  Results  are  presented  which  clearly  show  GGeneral  Electric's 
Laser  Velociraeter  to  be  a  viable  non-contact  type  probe  cabable  of  performing 
turbulence  rms  and  spectra  measurements  in  a  heated  high  velocity  jet,  and  to 
perform  in-jet  to  far-field  acoustic  noise  source  location  measurements. 

Detailed  discussions  and  results  are  also  presented  on  the  evaluation  of  in-jet 
static  pressure  fluctuations  probes  source  location  in  high  velocity  jets. 
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INTRODUCTION  | 
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I 

With  the  advent  of  larger  and  more  powerful  military  and  commercial  air¬ 
craft  propulsion  systems,  it  is  increasingly  apparent  that  to  improve  the 
general  community  environment  greater  efforts  must  be  taken  to  reduce  jet 
engine  noise.  A  considerable  effort  has  been  spent  in  the  last  twenty  years 
in  the  development  of  jet  noise  suppressors,  but  because  of  a  lack  of  clear 
understanding  and  detailed  mathematical  specification  of  the  dominant  noise 
producing  sources  necessary  for  the  establishment  of  meaningful  prediction 
procedures  for  even  the  most  simple  nozzles,  only  partial  success  has  been  met 
in  reducing  the  noise  with  acceptable  jet  nozzle  performance. 


The  overall  objective  of  this  joint  Air  Force  and  Department  of  Transpor¬ 
tation  Supersonic  Jet  Exhaust  Noise  investigation  is  to  develop  the  technology 
to  significantly  reduce  supersonic  aircraft  propulsion  system  noise  with 
minimum  associated  performance  and  weight  penalties.  To  reach  the  Program 
objectives  a  varied  and  comprehensive  research  program  is  being  carried  out 
to  develop  the  basic  theory  and  experimental  methods  necessary  for  under¬ 
standing  and  quantizing  the  acoustic  characteristics  of  simple  supersonic  jets 
over  a  range  of  velocities  and  temperatures  typical  of  present  and  future 
military  and  commercial  supersonic  aircraft  propulsion  systems. 


This  interim  report  summarizes  the  major  theoretical  and  experimental 
task  efforts  performed  at  General  Electric  during  the  first  year  of  a  two  year 
.program.  This  two  year  program  is  a  result  of  an  initial  exploratory  research 
program  initiated  by  the  Air  Force  Aero  Propulsion  Laboratory  during  FY  71. 

During  the  initial  phase  of  work  the  general  framework  for  General  Electric's 
method  of  approach  was  established.  At  the  time  of  the  initial  phase  of  work, 
a  thorough  review  of  various  competing  mathematical  models  used  to  explain 
supersonic  jet  noise  generation  processes  was  performed,  and  preliminary 
experimental  instrumentation  and  procedures  were  demonstrated.  One  result  of 
General  Electric's  initial  efforts  was  establishment  of  a  comprehensive  turbu¬ 
lent  mixing  aero-acoustic  model  capable  of  computing  all  the  main  acoustic  proper¬ 
ties:  overall  sound  power  level,  power  spectra,  overall  sound  pressure  level, 
sound  pressure  level  spectra  and  the  jet's  detailed  directivity  characteristics. 
Additionally,  the  procedure  was  extended  for  acoustic  noise  source  location 
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predictions  such  as  axial  power  distributions  ior  subsonic  and  supersonic 
exhaust  jets,  acoustic  peak  frequency  distribution,  and  the  effects  of  initial 
turbulence  intensity  on  jet  noise.  Hie  primary  emphasis  in  thinking  was  to 
c  imputations! 1 v  link  the  detailed  mean  and  turbulent  iiow  aerodynamic  properties 
to  selected  turbulent  mixing  acoustic  models.  This  computational  scheme  was 
completed  and  a  great  deal  of  success*’-  was  achieved  for  shock  free  supersonic 
high  temperature  jets.  The  computational  scheme  was  so  designed  as  to  enable 
the  acoustic  predictions  to  be  based  on  aerodynamic  input  which  could  be 
predicted  or  measured,  thus  allowing  the  scheme  to  be  compatible  with  exhaust 
nozzle  suppressor  investigations  where  the  detaiied  aerodynamic  properties  cannot 
as  yet  be  predicted. 

To  complement  the  theoretical  investigations,  experimental  instrumentation 
was  developed  to  measure  the  detailed  in-jet  flow  properties  of  heated  super¬ 
sonic  exhaust  jets.  It  was  demonstrated  that  General  Electric's  laser  veloci- 
meter  was  an  ideal  in-jet,  non-contact  type  probe,  capable  of  measuring  the 
detailed  mean  and  turbulence  velocity  of  supersonic  high  temperature  jets,  thus 
offering  the  option  of  using  measured  flow  properties  as  the  aerodynamic  input 
to  the  acoustic  prediction  models  developed. 

It  was  founi  during  the  first  phase  of  effort  that  the  original  models 
developed,  which  where  based  of  the  concepts  of  Lighthill,  Ribner  and  Ffowcs- 
Williams,  had  certain  limitations  and  difficulties  with  regard  to  properly 
accounting  for  the  density  dependence  of  jet  noise,  the  influence  of  moving 
sound  sources  on  the  acoustics  of  the  problem,  the  interplay  between  convective 
and  refractive  coupling  and  the  assessment  of  the  degree  each  plays  in  under¬ 
standing  the  directivity  characteristics  of  jet  noise,  and  the  influence  of 
shock  turbulence  interaction  on  the  aerodynamic  and  the  acoustic  properties 
of  high  velocity  and  high  temperature  supersonic  jets. 

Section  I  of  this  report  deals  with  the  aerodynamic  input  which  serves  as 
the  starting  point  for  the  acoustic  calculations  of  General  Electric's 
Comprehensive  Model.  The  shock  free  solutions  for  the  aerodynamic  input  have 
been  dealt  with  in  detail  in  the  first  phase  report. ^  Here  Section  I  deals 
with  inclusion  of  shock  waves  into  the  general  aerodynamic  flow  model.  A 
comprehensive  account  is  given  of  the  theoretical  and  computational  foundations 
for  formulating  and  computing  the  shock  structure  and  the  mean  and  turbulent 
flow  properties  of  heated,  shocked  supersonic  jets. 
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Section  II  deals  with  the  acoustic  input  of  the  aeroacoustic  model.  Here, 
detailed  discussions  and  theory  development  is  given  which  will  be  used  in  the 
refinement  of  the  acoustic  input  for  the  comprehensive  aero-acoustic  model 
development.  Since  the  predictive  capacity  of  the  model  has  been  well  demon¬ 
strated  in  the  Phase  1  report,  only  the  newer  theory  developments  of  jet  noise 
will  be  discussed.  Particular  emphasis  !  directed  toward  solutions  for  moving 
sources  and  the  new  insights  this  model  gives  for  explaining  the  reasons  of 
increased  low  frequency  convective  amplification  at  shallow  jet  angles,  the 
density  dependence  of  turbulent  mixing  noise,  and  correct  velocity  scaling  of 
jet  noise.  The  Philiips-Pao  turbulent  shear-layer  model  is  also  discussed. 
Particular  attention  is  devoted  toward  formulating  this  mcdel  in  a  way  acceptable 
for  aero-acoustic  predictions,  and  the  way  this  theory  accounts  for  convective/ 
refractive  coupling,  source  radiative  efficiency  and  the  temperature  dependency 
of  jet  noise.  Additionally,  a  section  is  devoted  to  forming  a  generalized 
Green's  function  approach  from  the  Lighthill/Ribner  jet  noise  theory  point  of 
view  to  unify  the  concepts  regarding  moving  and  stationary  noise  source  models, 
ana  how  the  theory  can  be  used  to  establish  the  framework  for  acoustic  refrac¬ 
tion  studies  that  will  be  carried  out  in  the  next  year's  effort. 

Section  III  gives  preliminary  results  from  a  set  of  detailed  far  field, 
near  field  acoustic  parametric  studies,  and  some  results  from  photographic 
studies.  The  results  are  presented  in  a  form  which  will  be  useful  for  future 
theory /data  comparisons  and  noise  source  location  inver  tigations. 

Section  IV  reports  on  the  recent  advances  made  in  developing  instrumen¬ 
tation  for  performing  in-jet  noise  source  location  studies  on  high  temperature 
high  velocity  exhaust  jets.  Section  IV  discusses  the  theory  and  error  analysis 
necessary  to  construct  turbulence  spectra  using  General  Electric's  laser 
velocimeter.  Results  of  demonstration  experiments  for  measuring  turbulence 
spectra  in  subsonic  cold  and  heated  high  velocity  jets  are  given.  Preliminary 
comparisons  of  LV  measured  turbulent  flow  properties  with  hot  film  measure¬ 
ments  are  also  discussed.  Additionally  a  detailed  discussion  of  in-let 
pressure  probe  and  in-jet  to  far-field  acoustic  probe  correlation  investigations 
is  given  for  a  sonic  jet.  How  these  studies  compare  with  work  performed  at 
low  velocities,  and  the  usefulness  of  in-jet  pressure  probes  for  high  speed  flow 
investigations  are  reviewed. 

Section  V  reviews  the  work  planned  for  the  completion  of  the  Program 
act ivities . 


3 


REFERENCES 


).  M.J.  Benzakein  and  P.R.  Knott,  "Supersonic  Jet  Exhaust  Noise,"  Air 
Force  Aero  Propulsion  Laboratory,  AFAPL-TR-72-52 ,  July  1972. 

2.  P.R,  Knott  and  M.J.  Benzakein,  "Analytical  and  Experimental  Supersonic 
Jet  Noise  Research,"  AIAA  73-188  (1973). 

3.  R.  Mani,  "A  Moving  Source  Problem  Relevant  to  Jet  Noise,"  GE  Class  I 
Report,  May  1972. 


4 


I.  THEORETICAL  DEVELOPMENTS  OF  THE  AERODYNAMICS  OF  SUPERSONIC  JETS 


C.L.  Merkle 

There  ere  two  distinct  techniques  which  have  been  used  to  analyze  the 
aerodynamic  flow  field  in  a  supersonic  jet.  In  the  first  approach,  the  jet  is 
treated  as  a  viscous,  boundary  layer  flow.  The  resulting  flow  field  is  of  the 
type  depicted  in  Figure  1.  According  to  the  usual  boundary  layer  approxima¬ 
tions,  the  radial  velocity  components  are  assumed  small  in  comparison  to  their 
axial  counterparts,  and,  in  addition,  the  pressure  is  taken  to  be  constant 
throughout  the  whole  flow  field.  These  approximations  implicity  assume  that 
the  static  pressure  at  the  jet  exit  plane  is  identical  to  the  ambient  pressure 
and  that  Prandtl-Meyer  expansions  and/or  shock  waves  are  not  present  in  the 
flow  field.  Consequently,  this  viscous  boundary  layer  analysis  can  only  be 
applied  to  subsonic  jets,  or  to  supersonic  jets  which  are  ideally  (or  nearly 
ideally)  expanded. 

contrast  to  this  viscous  analysis,  the  second  traditional  technique 
for  analyzing  supersonic  jets  completely  ignores  the  effects  of  turbulent  mixing. 
In  this  second  (inviscid)  analysis,  the  full  two-dimensional  equations  of 
motion  are  used,  and  strong  radial  and  axial  pressure  gradients  can  occur. 

These  pressure  gradients  have  their  origin  at  the  nozzle  exit  plane  where  the 
static  pressure  is  generally  significantly  different  from  the  ambient  pressure. 

In  adjusting  to  the  ambient  pressure,  the  flow  field  generally  developes  a 
series  of  shock  waves  and  Prandtl-Meyer  expansions  in  a  nearly  periodic,  cell¬ 
like  fashion.  A  schematic  description  of  the  qualitative  features  of  a  jet 
described  by  this  two-dimensional  analysis  is  shown  on  Figure  2. 

As  indicated  above,  both  of  these  approximate  models  are  applicable  to 
the  analysis  of  a  certain  class  of  supersonic  jet.  However,  as  might  be 
expected,  neither  model  applies  to  all  supersonic  jets.  Thus,  for  example,  the 
effects  of  friction  can  never  be  entirely  removed  from  the  jet.  Further, 
supersonic  jets  are  seldom  uniform,  parallel  ideally  expanded  jets.  Consequently, 
in  order  to  obtain  an  acoustic  prediction  technique  which  is  applicable  to 
both  ideally  expanded  and  non-ideally  expanded  jets,  the  aerodynamic  model 
ni us t  include  both  two-dimensional  effects  and  viscous  mixing  effects.  During 
the  present  contract  reporting  period,  we  have  developed  such  an  aerodynamic 
model.  The  details  of  this  analysis  are  given  in  the  following  sections. 
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1.  DESCRIPTION  OF  THE  AERODYNAMIC  MODEL 


Exhaust  nozzles  of  most  contemporary  gas  turbine  engines  generally  operate 
near  their  ideal  expansion  ratios.  In  view  of  this,  General  Electric's  initial 
efforLs  at  the  prediction  of  the  sound  field  of  a  supersonic  jet  relied  on  an 
aerodynamic  analysis  which  included  only  the  effects  of  turbulent  mixing. 
Specifically,  this  aerodynamic  model  was  of  the  viscous,  boundary  layer  type 
described  above.  The  computerized  version  of  this  analysis  is  referred  to  as 
the  JETMIX  computer  program.  This  compiler  program  solves  the  time-averaged 
turbulent  boundary  layer  equations  using  boundary  conditions  which  are  appro¬ 
priate  for  free  jets.  The  turbulent  Reynold's  stresses  are  included  by  means 
of  a  turbulence  model  which  is  based  on  a  turbulent  kinetic  energy  concept. 

This  turbulence  model  is  based  on  those  developed  by  Rotta^  ,  Glushko^ 

(3) 

and  Spalding  .  Details  of  the  turbulence  model  and  the  JETMIX  computer 
program  as  well  as  extensive  comparison  of  the  predictions  with  experimental 
data  are  given  in  References  4  and  5. 

Starting  from  this  basic  viscous  analysis,  we  have  extended  our  aero¬ 
dynamic  model  to  enable  it  to  predict  non-ideally  expanded  jet  flow  fields. 

This  improved  aerodynamic  capability  allows  us  to  predict  the  effects  of 
non-ideal  expansion  on  the  acoustic  field  of  a  supersonic  jet.  In  particular, 
the  analysis  enables  us  to  estimate  the  difference  between  the  sound  which  is 
generated  by  an  ideally-expanded  jet  from  a  convergent-divergent  nozzle  and  a 
highly  underexpanded  jet  of  the  same  Mach  number  from  a  convergent  nozzle. 

The  method  which  we  have  used  to  include  the  two-dimensional  effects  which 
occur  in  non-ideally  expanded  jets  is  based  on  dividing  the  jet  into  an  inner 
region  and  an  outer  region  as  shown  in  Figure  3.  The  outer  region  of  the 
jet  contains  that  part  of  the  jet  in  which  the  effects  of  turbulent  mixing  are 
significant.  Near  the  nozzle  exit  the  outer  region  is  composed  of  a  narrow 
annular  portion  of  the  flow  field  on  the  outer  edge  of  the  jet;  downstream  of 
the  exit  plane,  the  thickness  of  the  outer  region  increases  until  eventually  it 
includes  the  entire  jet.  In  our  analysis,  this  outer  region  is  computed  by  our 
original  viscous,  boundary  layer  (JETMIX)  computer  program.  Now,  whereas  the 
outer  region  of  the  jet  is  dominated  by  the  effects  of  viscous  mixing,  the 
inner  region  of  the  jet  is  dominated  by  the  familiar  Prandtl-Meyer  expansions 
and  shock  waves  which  characterize  two-dimensional  supersonic  flow  fields.  In 
order  to  include  these  effects  in  our  aerodynamic  model,  a  new  computer  program 
has  been  written  to  handle  this  inner  region.  This  new  program  is  called  the 


Supersonic  Finite  Difference  (SSFD)  computer  program.  Urns  in  the  refined 
aerodynamic  moJel  ,  a  supersonic  jet  is  analyzed  in  two  parts  by  two  separate 
computer  programs.  The  outer  portion  of  the  jet  is  analyzed  by  the  viscous 
JFTMIX  analysis.  Hie  inner  portion  of  the  jet  is  computed  by  the  two-dimensional 
SSFD  analysis. 

As  indicated  in  Figure  3,  these  two  separate  parts  of  the  flow  field 
are  matched  along  the  sonic  line.  Thus  the  inner  portion  of  the  flow  field  is 
supersonic  while  t he  outer  flow  is  subsonic.  (In  actuality  some  constant  Mach 
number  line  which  is  slightly  supersonic  is  chosen  as  the  matching  line  rather 
than  precisely  the  sonic  line).  However,  it  must  be  noted  that  the  sonic  line 
appears  in  the  jet  because  viscous  effects  have  reduced  the  Mach  number  of  the 
formerly  supersonic  flow.  This  indicates  that  the  outer  edge  of  the  supersonic 
region  has  experienced  considerable  viscous  effects.  Thus  in  order  to  include 
the  two-dimensional  effects  in  as  large  a  region  as  possible,  and  in  order  to 
enforce  as  smooth  a  match  as  possible  between  the  inner  and  outer  solutions, 
the  effects  of  the  viscous  mixing  are  included  in  the  inner  (SSFD)  analysis 
as  known  ''right-hand-side"  terms.  The  magnitude  of  the  "right-hand-side"  terms 
is  estimated  from  the  viscous  JETMIX  computer  program  as  indicated  later.  This 
matching  technique  allows  the  total  pressure  to  vary  continuously  from  the 
outer  edge  of  the  jet  (where  the  flow  is  essentially  stagnated)  through  the 
sonic  line  and  all  the  way  to  the  jet  centerline  (where  the  flow  is  supersonic). 
Then  by  matching  the  static  pressure  at  the  sonic  line,  we  can  be  sure  that  all 
flow  properties  are  continuous  at  the  matching  line. 
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2 .  DERIVATION  OF  THE  EQUATIONS  USED  IN  THE  INNER  REGION 

As  indicated  above,  the  equations  used  in  the  inner  region  include  the 
viscous  effects  as  "right-hand-side"  terms.  In  order  to  obtain  the  form  of 
these  "right  hand  side"  terms,  the  equations  for  the  inner  region  are  obtained 
from  the  complete  Navier-Stokes  equations.  An  outline  of  the  derivation  follows. 

The  equations  of  motion  for  steady,  compressible,  viscous  flow  are 
V  -p  v  =  o 

p  (v»V)  v  +  Vp  =  V*  t 
p  v  • Ve  =  t:  Vv  -p7*v  -  V-q 

Two  vector  identities  which  are  useful  are 
v-(vV)  v  =  v  .  V  (v^/2) 

v (V  •  t)=  V.  (v  .7)-  t :  Vv 

If  we  dot  Equation  (4)  by  the  velocity  vector,  V,  and  use  identity  (4)  , 
we  obtain 


(1) 

(2) 

(3) 

(A) 

(5) 


P  vV  (v^/2)  +  v*Vp  =  v*  (V.  t)  (g 

Then,  combining  Equations  (1)  and  (3)  and  using  identity  (5),  the  energy 
equation  becomes 


p  v  .  Vh  -  v*Vp=V*(v  .t)-v*(V  .x)-V-q  (7) 

where  we  have  also  converted  from  internal  energy  to  enthalpy.  Then,  adding 
Equations  (6)  and  (7)  gives 

p  v*Vh°=  V*(v  • t)  -V*q  (8) 
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We  now  define  the  scaler  function,  Q,  as 

Q(x,y)  =  pv*Vh°  (9) 

and  by  Equation  (8) ,  we  also  have 

Q(x,y)  =  V-(vt)  -  V*q  (10) 

If  we  now  combine  Equation  (10)  with  Equation  (7)  and  use  the  thermodynamic 
equation  of  state,  the  energy  equation  becomes 

pT  (v-VS)  =  -v*(V*t)  +  Q  (x,y)  (11) 

Finally,  defining  the  scalar  function,  as 

<f>  (x, y)  =  v  (V* t)  (12) 

we  can  write  the  entropy  variation  along  a  streamline  as 

pT  v-VS  =  -<t>  (x,y)  +  Q  (x,y)  (13) 

Thus,  the  final  version  of  the  equations  of  motion  (1),  (2),  and  (3)  can  be 
re-written  as 

V’pv  =  o  (14) 

p  (v*  V)  v  +  Vp  =  R  (15) 

pT  v  VS  =  -vR  +  Q  (x,y)  (16) 

where  the  vector  R  is  defined  as 


R  =  V*  t 
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so  that 


if>  (x,y)  =  v-R  (18) 

These  are  the  equations  which  are  solved  by  the  SSFD  computer  program.  At 
present  the  function,  Q(x,y)  ,  is  limited  to  the  trivial  function 

Q(x,y)  =  0  (19) 

This  implies  that  only  flow  fields  which  have  uniform  total  temperature 
throughout  can  be  calculated.  The  extension  of  the  computer  program  to  include 
an  arbitrary  specification  of  the  stagnation  enthalpy  is  relatively  simple. 

Note  that  the  function,  y,  is  not  restricted;  it  can  (in  principle)  be  any 
function.  When  coupled  to  the  JETMTX  viscous  analysis,  the  SSFD  program 
automatically  determines  <j>  from  the  JETMTX  -  predicted  entropy  gain  due  to  the 
turbulent  stresses. 

Equations  (14)  ,  (15),  (16)  and  (19)  are  solved  numerically  using  a  finite- 
difference  algorithm  developed  by  MacCormick^  . 
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CALCULATION  OF  THE  TURBULENCE  FIELD 


The  previous  paragraphs  have  described  the  model  which  is  used  to  predict 
the  velocity  field  in  an  off-design  jet.  However,  before  the  acoustic  charac¬ 
teristics  of  the  jet  can  be  determined,  it  is  necessary-  to  know  something  of 
the  turbulence  field  in  the  jet.  As  indicated  above,  the  turbulence  model 
which  is  used  in  the  JETMIX  computer  program  is  based  on  a  turbulent  kinetic 
energy  approach.  For  ideally  expanded  jets,  the  magnitude  of  this  turbulent 
kinetic  energy  has  been  used  to  evaluate  the  source  terms  in  the  classical 
Lighthill  acoustic  equations  (5,7).  Once  these  source  terms  are  evaluated, 
the  acoustic  signature  of  the  jet  can  be  readily  determined.  Since  this 
acoustic  formulation  is  based  on  the  local  mean  and  fluctuating  properties  of 
the  jet,  it  should  also  be  directly  applicable  to  non-ideally  expanded  jets. 
(This,  of  course,  does  not  imply  that  the  model  would  predict  the  same  acoustic 
radiation  from  an  ideally  expanded  jets,  because  both  the  mean  velocity  field 
and  the  turbulence  field  depend  on  the  expansion  ratio  of  the  jet.)  Thus,  it 
remains  to  determine  the  turbulent  kinetic  energy  in  the  non-ideally  expanded 


The  conservation  of  turbulent  kinetic  energy  is  governed  by  the  balance 
between  between  the  production,  dissipation,  convection,  and  diffusion  of 
turbulence  energy  throughout  the  flow  field.  The  form  of  the  turbulent  kinetic 
energy  equation  which  is  used  in  the  JETMIX  analysis  is 


„  3k  ,  3k  13  3kx  .  3u 

pU  ta  +  DV  ' 717  (Cl0“u  P  +  Ul  s? 


For  clarity,  the  physical  meaning  of  each  of  these  terms  is  labeled.  As  used 
in  the  JETMIX  analysis,  this  equation  basically  applies  to  ideally  expanded 
jets.  However,  the  local  properties  of  the  turbulence  should  not  depend  on 
whether  or  not  the  jet  is  ideally  expanded.  Consequently  we  take  equation  (20) 
as  being  applicable  in  non-ideally  expanded  jets  also. 


Although  we  use  the  same  turbulence  conservation  equation  in  both  ideally 
and  non-ideally  expanded  jets,  there  is  one  flow  phenomenon  which  affects  the  tur¬ 
bulence  levels  and  which  is  unique  to  non-ideally  expanded  jets  that  is  not  in¬ 
cluded  in  the  turbulence  conservation  equation  (20).  This  phenomenon  is  the 
presence  of  shock  waves  in  the  flow  field.  These  effects  have  been  included  in 


i  t  unaly.-  is  by  means  ut  Ribnei's  Slunk Turbulence  interaction  theory  (8,9). 

Ribner's  analysis  starts  by  ile  compos  i  ng  the  turbulence  field  into  an 
infinite  number  of  elementary  vorticity  waves  of  all  wavelengths  and  orienta¬ 
tions.  Then  ior  any  one  of  these  elementary  waves,  he  calculates  the  manner 
in  which  the  vorticity  of  the  wave  is  altered  as  it  is  convected  through  a 
normal  shock  wave.  Hu-  results  of  his  calculation  show  that  the  magnitude  of 
the  vorti'  ,‘ty  is  increased  as  the  wave  goes  through  the  shock.  (Besides  the 
increased  vorticity,  two  new  wa.es  are  generated,  an  entropy  wave  and  an 
acoustic  wave.)  A  summation  over  all  wave  numbers  of  the  effects  of  the  shock 
on  each  individual  wave  then  yields  an  amplification  factor  for  the  turbulence 
as  it  is  transmitted  through  the  shock.  Conversion  from  turbulence  convected 
through  a  normal  shock  to  turbulence  convected  through  an  oblique  shock  is 
made  by  a  transformation  of  coordinates.  Finally,  it  should  be  noted  that 
although  Ribner's  analysis  strictly  applies  only  to  straignt  shocks,  it  can 
also  be  applied  to  curved  shocks  (such  as  occur  in  supersonic  jets)  as  long 
as  the  radius  of  curvature  of  the  shock  is  significantly  larger  than  the 
longest  wavelength  of  the  turbulence. 

Figure  4  shows  the  amplification  of  turbulence  by  a  shock  in  terms  of 
the  ratio  of  turbulent  kinetic  energy  in  front  of  and  behind  the  shock.  The 
turbulence  amplification  is  plotted  as  a  function  of  the  ratio  of  the  normal 
components  of  velocity  in  front  of  and  behind  the  shock.  As  can  be  seen,  the 
amplification  is  unity  at  a  velocity  ratio  of  unity  (shock  of  vanishing  strength) 
but  quickly  increases  to  a  maximum  of  some  20Z  amplification  for  moderate  shock 
strengths  (normal  component  of  incoming  Mach  number  about  1.5). 

In  our  computer  model,  the  turbulent  kinetic  energy  is  monitored  at  each 
point  in  space  by  means  of  equation  (20).  The  source  term  in  this  equation 
accounts  only  for  the  production  of  turbulence  by  the  viscous  shear  forces. 

When  a  shock  wave  is  encountered,  the  turbulence  amplification  is  determined 
from  Ribner's  theory  and  the  turbulent  kinetic  energy  is  increased  locally  by 
the  amount  of  turbulence  which  is  generated  at  the  shock.  The  resulting 
turbulence  energy  profiles  have  a  discontinuous  jump  across  the  shock.  The 
magnitude  of  this  jump  is  determined  from  Ribner's  theory. 
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4.  MATCHING  BETWEEN  THE  INNER  AND  OUTER  SOLUTIONS 

As  described  above,  the  equations  for  the  inner  region  require  that  the 
variations  in  entropy  and  stagnation  enthalpy  due  to  the  turbulent  shearing 
stresses  be  specified  along  each  streamline  before  the  solution  is  calculated. 
In  the  computational  procedure,  these  variations  are  first  estimated  from  a 
solution  of  the  viscous  boundary  layer  equations  for  the  entire  mass  flow  in 
the  jet  (i.c.,  both  inner  and  outer  flows).  This  viscous  calculation  also 
establishes  the  value  of  the  stream  function  at  which  the  flow  becomes  sonic, 
and  hence  the  location  of  the  boundary  between  the  inner  and  the  outer  regions. 
Then,  using  the  calculated  entropy /enthalpy  variations  as  a  first  approximation 
to  their  actual  behavior,  the  velocity  field  in  the  inner  region  is  re-calcu¬ 
lated  by  means  of  the  two-dimensional  (inner)  equations.  (The  velocity  field 
in  the  outer  region  is  left  unchanged  except  for  re-positioning  the  streamlines 
so  that  they  match  with  the  streamlines  in  the  inner  region.)  In  principle, 
this  sets  up  an  iterative  process  which  could  be  continued  by  using  the 
predicted  pressure  gradient  as  an  impressed  static  pressure  field  for  the 
outer  equations.  Thus  upon  convergence,  the  "exact"  solution  would  be  obtained 
(except  that  the  viscous  terms  would  be  included  only  to  the  boundary  layer 
approximation).  Note  that  the  iteration  would  proceed  by  assuming  that  a 
known  static  pressure  field  is  impressed  on  the  outer  (boundary  layer)  equa¬ 
tions.  The  boundary  layer  calculation  then  defines  an  entropy /enthalpy  field 
which  is  impressed  on  the  inner  (two-dimensional  "inviscid")  equations  which, 
in  turn,  re-define  the  static  pressure  field,  and  so  forth.  Nevertheless, 
for  the  problem  at  hand,  it  is  assumed  that  the  use  of  a  constant  pressure 
field  in  the  boundary  layer  equations  will  give  the  entropy  /enthalpy  field  to 
sufficient  accuracy  that  an  improved  approximation  need  not  be  determined. 

The  computational  procedure  is  described  schematically  on  Figure  5. 

The  boundary  conditions  along  the  sonic  line  complete  the  matching  of 
the  inner  and  outer  solutions.  The  viscous  solution  assumes  the  static 
pressure  is  constant  throughout  the  outer  region  and  equal  to  the  ambient. 

At  the  matching  (sonic)  line,  the  static  pressure  in  the  inner  region  is 
required  to  approach  the  ambient  pressure.  Thus  by  requiring  the  static 
pressure  to  be  continuous  across  the  sonic  line,  and  by  obtaining  the  entropy/ 
enthalpy  field  for  the  en  :ire  jet  from  the  outer  solution,  we  are  assured 
that  all  other  flow  and  . t.ermodynamic  properties  are  continuous  at  the  inter¬ 
face  also. 
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Now  since  the  entire  total  pressure  field  is  obtained  from  the  outer 
(\iscous)  solution,  and  since  the  static  pressure  is  required  to  be  continuous 
along  the  sonic  line,  then  all  other  flow  and  thermodynamic  properties  will 
be  i  ontinuous  at  the  matching  line  also  for  the  case  of  the  cold  jet.  To 
achieve  this  complete  matching  of  all  properties  in  a  hot  jet  case,  the  total 
temperature  along  each  streamline  would  have  to  be  superimposed  on  the  inner 
solution  (in  addition  to  the  total  pressure).  As  mentioned  above,  the  computer 
program  does  not  currently  have  capability  for  variable  total  temperature, 
although  all  the  necessary  mathematical  requirements  have  been  worked  out. 

Finally  it  should  be  pointed  out  that  the  viscous  boundary  layer  equations 
which  are  solved  by  the  JETMIX  computer  program  are  parabolic  in  the  axial 
coordinate,  x,  and  that  the  solution  can,  therefore,  be  obtained  by  a  matching 
process.  A  similar  matching  process  can  be  used  to  solve  the  inner,  two- 
dimesnionil  equations  because  these  equations  are  hyperbolic  (so  long  as  the 
flow  remains  supersonic).  Consequently  both  the  inner  and  outer  solutions  lend 
themselves  to  a  matched,  coupled  calculation  such  as  has  been  described. 

5.  SHOCK  REFLECTION  FROM  AXIS  OF  SYMMETRY 

As  a  shock  wave  in  an  axisymmetric  flow  field  approaches  the  centerline, 
the  shock  becomes  increasingly  steeper.  Because  of  this  steepening,  the 
axisymmetric  equations  will  not  allow  the  shock  to  reflect  from  the  symmetry 
axis  in  a  regular  fashion.  Instead,  seme  sort  of  "strong"  reflection  must 
occur,  As  a  result,  a  local  pocket  of  subsonic  flow  appears  behind  the  shock 
and  any  computational  procedure  which  relies  on  the  hyperbolic  character  of 
the  equations  becomes  invalid  and  has  to  be  terminated.  However,  experimental 
schlieren  photographs  show  that  this  subsonic  region  is  frequently  small  or 

(13) 

even  non-existent.  Thus,  although  the  reflection  shows  up  as  a  nearly 
normal  "Mach  disc"  or  "Riemann  wave"  in  some  cases,  an  apparently  regular 
reflection  takes  place  in  other  cases  when  the  shock  is  sufficiently  weak, 

(even  though  the  inviscid  equations  will  not  allow  this).  In  order  to  provide 
a  means  for  continuing  the  flow  field  calculation  beyond  the  location  at  which 
the  shock  first  reflects  from  the  axis  of  symmetry,  we  have  incorporated  two 
approximate  techniques  for  "calculating  through"  this  presumably  small, 
localized  subsonic  pocket.  First,  when  the  incoming  shock  is  weak,  a  "regular" 
reflection  procedure  is  used.  However,  for  stronger  incoming  shocks,  we  switch 
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over  to  1  "Mach  disc"  reflection  procedure.  Wlii^n  of  the  two  techniques  is 
to  be  use  1  i.iur.t  be  determined  by  ihe  problem  ;it  hand. 

Tli»  "icgular"  reflection  procedure  utilizes  a  suggestion  by  Oswatitsch 
(14)  th  ti  '  axis  of  svnmetry  be  "enlarged"  near  the  shock  impingement  point 
so  that  ti...  radial  coordinate  becomes  ^smal  1  but  still  remains  finite.  The 
radial  size  of  this  "enlargement"  is  determined  by  the  program  depending  on 
the  local  strength  oi  the  shock.  (Stronger  shocks  require  more  "fattening" 
of  the  axis  of  symmetry.)  It  is  emphasized  that  these  "enlargements"  generally 
encompass  less  than  one  half  of  one  percent  of  the  original  mass  flow  so  that 
they  are  scarcely  detectable  on  a  "blown  up"  plot  of  the  shock  locus. 


The  "Mach  disc"  reflection  technique  (which  is  considerably  more  compli¬ 
cated  than  the  "regular"  reflection  procedure)  involves  the  insertion  of  a 
triple  point  and  the  use  of  an  iterative  technique  to  determine  its  location. 
In  this  analysis,  a  triple  point  is  inserted  at  a  chosen  location  on  the 
shock,  and  the  oblique  shock  which  is  moving  radially  inward  is  forced  to 
branch  into  a  second  outward- running  shock  and  a  normal  shock  which  extends 
to  the  axis  as  shown  in  Figure  6.  The  normal  shock  represents  the  Mach 
disc.  A  slip  line  is  also  generated  at  the  "lambda"  intersection.  Downstream 
of  the  Mach  disc,  this  slip  line  serves  as  a  boundary  between  the  supersonic 
flow  and  the  subsonic  flow.  The  supersonic  flow  is  handled  by  the  standard 
SSFD  algorithm,  while  the  subsonic  flow  is  analyzed  by  a  one-dimensional 
approximation.  The  height  of  this  one- dimensional  channel  at  succeeding  axial 
locations  is  determined  by  requiring  the  pressure  to  be  balanced  across  the 
slip-stream,  and  by  requiring  the  supersonic  flow  to  be  tangent  to  the  slip 
line.  Tli is  matching  requirement  causes  the  Mach  number  in  the  one-dimensional 
stream  to  vary  as  it  flows  downstream.  The  axial  position  of  the  Mach  disc 
is  then  iteratively  determined  based  on  the  behavior  of  the  flow  in  this 
one-dimensional  channel.  The  Mach  disc  is  said  to  have  been  correctly 
positioned  when  the  slip  line  forms  a  "throat"  which  re-accelerates  the  sub¬ 
sonic  flow  through  sonic  velocity  in  a  smooth,  continuous  fashion.  This 
Mach  disc  model  is  very  similar  to  the  ones  used  by  Abbett  (15),  Averenkova, 
et  al.  (16)  and  (17).  Comparisons  between  this  Mach  disc  model  and  experi¬ 
mental  results  have  shown  reasonable  agreement,  but  the  iterative  procedure 
is  quite  expensive  (in  terms  of  computer  processing  time)  and  tends  to  be 
unreliable. 
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Calculation  of  Jet  Noise  From  Turbulent  Kinetic  Ener 


As  indicated  above,  once  the  turbulent  kinetic  energy  is  known,  the 
acoustic  radiation  from  an  off-design  jet  can  be  computed  directly  as  was 
done  for  ideally  expanded  jets  in  references  1-5  and  1-7.  Note  that  this 
calculation  automatically  takes  into  account  the  indirect  effect  of  the  shock 
on  the  jet  noise.  That  is,  the  local  turbulence  level  is  increased  by  the 
shock  and  so  the  local  acoustic  radiation  is  similarly  increased  by  the 
presence  of  the  shock. 

Nevertheless,  the  presence  of  shock  waves  in  a  non-ideally  expanded 
jet  will  not  necessarily  lead  to  a  higher  predicted  level  cf  noise  (as 
compared  to  the  corresponding  ideally  expanded  jet)  even  though  the  shock 
wave  acts  to  increase  the  turbulence  (and  hence  the  noise).  The  reason  for 
this  is  that  the  corresponding  change  in  the  mean  velocity  field  will  also 
have  an  effect  on  the  turbulence  levels  which  could  tend  to  off-set  the 
generation  of  turbulence  by  the  shocks. 

In  addition  to  this  indirect  effect  of  shocks  on  the  jet  noise,  there 
is  also  a  direct  effect.  In  addition  to  generating  additional  turbulence, 
the  shock  also  generates  acoustic  waves  directly.  These  shock-generated 
acoustic  waves  are  of  two  different  natures;  one  subsonic  and  the  other 
supersonic. 

Those  waves  which  are  subsonic  in  nature  decay  exponentially  with 
distance  behind  the  shock,  and  so  can  be  neglected  in  the  far  field.  The 
supersonic  waves  propagate  to  infinity  (as  plane  waves)  and  so  are  not 
negligible  in  the  far  field.  The  proportion  of  subsonic  and  supersonic  waves 
depends  on  the  Mach  number  of  the  incoming  flow  and  on  the  particular 
(vectorial)  Fourier  component  of  the  incoming  turbulence. 

As  an  estimate  of  the  magnitude  of  the  acoustic  intensity  which  is 

generated  by  a  shock-turbulence  interaction,  the  shock-generated  T  oise  for 

a  sonic  jet  which  was  under-expanded  by  the  ratio  P.  /P  ,  =■  2  has  been 
J  jet  ani 

calculated.  For  this  jet  (in  which  the  strength  and  location  of  the  shock 
was  determined  by  means  of  the  SSFD  computer  program)  the  acoustic  intensity 
of  the  noise  generated  by  the  first  cell  of  the  shock  (shock  wave  running 
from  the  near  the  outer  boundary  of  the  jet  in  towards  the  axis  of  symmetry, 
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and  then  reflecting  back  to  the  outer  edge)  was  73  decibels,  the  reference 

-12 

level  for  the  energy  being  10  times  the  total  flow  energy  in  the  jet,  i.e. 

Shock  Generated  -  io  1  pu  ^shock  ^acoustic 

Acoustic  Energy  3  . 

puJ  ''jet  x  10 

And  the  acoustic  intensity,  1  ,  is  given  in  terms  of  the  incoming  turbulence 

A 

intensity  via  Ribner's  theory  (10).  These  results  indicate  that  the  acoustic 
intensity  generated  by  shock-turbulence  interaction  is  small  compared  to  the 
overall  acoustic  energy  of  a  supersonic  jet,  but  that  it  is  sufficiently  large 
to  have  a  non-negligible  effect  on  the  acoustic  signature  of  the  jet.  Finally, 
it  should  be  noted  that  this  calculation  was  made  for  a  jet  having  an  initial 
turbulence  intensity  oi  10%  at  the  jet  exit. 

In  addition  to  this  indirect  effect  on  the  jet  noise,  the  shock  also  has 

a  direct  effect  on  the  noise  because  of  the  acoustic  waves  which  are  produced 

by  the  shock  turbulence  interaction.  This  direct  source  o*  noise  must  be 

"added  on"  to  the  turbulence-produced  noise  in  the  jet.  The  magnitude  of  this 

(10) 

shock-produced  noise  is  also  given  by  Ribner 

Finally,  it  should  be  noted  that  the  acoustic  radiation  which  is  predicted 
by  this  shock-turbulence  interaction  model  is  broadband  in  nature.  As  is 
well-known,  shock  waves  add  both  broadband  noise  and  narrowband  noise  to  the 
overall  jet  acoustic  pattem^’^ .  The  above  theory  does  not  include  a 
mechanism  for  estimating  either  the  magnitude  or  the  frequency  of  these  narrow- 
band  (screech)  tones. 


■jyjlli  j,  )|ii|Www^fPyT^jill|ii.i,)i,li  PU'l.'iJM!  1 1 
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6.0  RESULTS  FOR  SHOCK  STRUCTURE  COMPUTATIONS 
Inviscld  Calculations 

Some  typical  predictions  of  the  aerodynamic  flow  field  are  given  in 
Figures  7  through  11.  The  results  in  these  figures  have  been  obtained 
from  completely  inviscid  calculations.  Figure  7  shows  the  predicted  shock 
shape  and  outer  boundary  shape  based  on  the  inviscid  calculation.  The  shock 
originates  near  the  outer  edge  of  the  jet  due  to  coalescing  characteristics 
coming  from  the  curved  outer  boundary.  The  shock  moves  radially  inward  and 
eventually  reflects  from  the  axis  of  symmetry  and  returns  to  the  outer  boundary. 
The  "regular"  reflection  technique  has  been  used  in  this  case.  Figures  8  and 
9  show  the  composite  result  of  a  number  of  .omputations  similar  to  that  of 
Figure  7.  In  Figure  8  is  plotted  the  distance  from  the  nozzle  exit  to  the 
point  at  which  the  shock  first  crosses  the  axis  of  symmetry  as  a  function  of 
pressure  ratio,  p  These  results  are  for  both  "regular"  and  Mach  disc 

reflection.  Also  shown  on  Figure  8  is  a  line  representing  the  experimental 
data  -f  Love^^^  .  The  inviscid  predictions  agree  quite  well  with  the  experi¬ 
mental  results;  however,  this  is  to  be  expected  since  the  viscous  effects 
don't  start  to  have  significant  effects  on  the  shock  shape  until  after  it 
reflects  from  the  axis  and  nears  the  outer  boundary.  Figure  9  is  similar  to 
Figure  8  except  that  it  shows  the  height  of  the  Mach  disc  as  a  function  of 
pressure  ratio.  Again,  Love's  experimental  data  are  shown  for  comparison. 

The  predicted  Mach  disc  heights  are  in  only  fair  agreement  with  the  experimen¬ 
tally  observed  values.  Nevertheless,  the  qualitative  agreement  is  sufficient 
to  show  that  the  Mach  disc  model  can  be  used  as  an  artifice  to  allow  the 
two-dimensional  supersonic  flow  calculation  to  proceed  beyond  the  point  where 
the  shock  hits  the  axis. 

An  overlay  of  the  inviscid  shock  shape  prediction  of  Figure  7  with  a 
schlieren  photograph  taken  under  the  experimental  portion  of  this  contact  is 
shown  as  Figure  10.  The  agreement  between  the  computed  shock  and  the 
experimental  shock  is  excellent  except  for  two  points.  First  of  all,  the 
computed  shock  starts  considerably  closer  to  the  nozzle  than  does  the  experimental 
shock.  However,  the  computed  version  of  the  shock  represents  a  Mach  number 
jump  of  only  about  0.02  until  very  near  the  centerline.  A  shock  this  weak 
would  not  be  expected  to  show  up  on  a  schlieren  photograph.  Secondly,  the 
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predicted  shock  does  not  turn  normal  Lo  the  flow  near  the  edge  of  the  jet. 
This  difference  is  due  to  the  i teg, Led  of  the  viscous  effects  in  the  outer 
region  of  the  jet.  Some  calculations  which  do  include  the  effects  of  viscous 
nixing  are  shown  in  the  next  section. 


Turbulent  How  l  ic  Id  Cal  cuiat  i  ons 

Figures  .11  through  IS  present  the  results  of  calculations  based  on 

the  full  coupled  viscous  -  inviscid  analysis.  Figures  11,  12  and  13 

show  the  effect  of  varying  amounts  of  underexpansion  on  a  jet  plume.  In  all 

three  figures,  the  totai-to-ambient  pressure  ratio,  P„/P  ,  =  4.10.  The  static- 

1  amb 

to-ambient  ratio,  P,/P  ,  ,  is,  however,  different  in  each  figure.  Figure  11 
shows  the  radial  variation  of  both  the  total  pressure  and  the  static  pressure 
for  the  ideally  expanded  jet,  P  /P  ^  =  1.0.  Here,  the  pressure  is  constant 
(and  equal  to  the  ambient)  throughout  the  entire  jet.  Consequently,  both  the 
complete  inner-outer  analysis  and  the  purely  viscous  boundary  layer  analysis 
give  identical  results  for  this  case.  At  the  axial  locations  shown,  x/R  =  1.90 
and  2.65,  the  total  pressure  near  the  centerline  of  the  jet  has  remained  at 
its  original  upstream  value  indicating  that  the  inviscid  core  is  still  present. 
Near  the  outer  edge  of  the  jet,  the  total  pressure  falls  off  quite  rapidly  due 
to  mixing.  This  decrease  continues  until  the  total  pressure  approaches  the 
static  (ambient)  pressure  signifying  thac  the  velocity  has  dropped  to  zero. 


A  slightly  unde rexpanded  jet  (P./P  ,  =  1.6)  is  shown  in  Figure  12. 

This  flow  field  contains  a  weak  shock,  which,  at  the  axial  location  shown, 
x/R  =  2.65,  has  just  reflected  from  the  axis  of  symmetry  and  is  moving  back 
toward  the  outer  boundary.  Because  of  the  shock,  there  are  now  two  sources 
of  total  pressure  loss.  Since  the  shock  has  already  reflected  from  the  axis 
of  symmetry,  the  flow  in  the  center  of  the  jet  has  experienced  a  finite,  shock- 
induced  total  pressure  loss  as  shown  by  the  smaller  shaded  region.  Between 
this  region  and  the  outer  mixing-loss  region  (also  shown  shaded)  lies  a  portion 
of  the  gas  which  is  unaffected  by  mixing  and  has  been  traversed  by  only  a  very 
weak  shock  so  that  its  total  pressure  remains  equal  to  its  upstream  value.  The 
radial  variation  of  the  static  pressure  is  no  longer  trivial  in  this  case  as  it 
was  in  Figure  11.  The  pressure  near  the  center  is  relatively  high,  then 
drops  across  the  shock  to  a  below-ambient  value  and  finally  asymptotically 
approaches  the  ambient  value  at  the  interface  between  the  inner  and  outer  regions. 


! 
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The  location  of  this  interface  as  well  as  the  location  of  the  sonic  point  are 
also  shown  on  Figure  12. 

The  last  figure  of  this  series  represents  a  still  larger  degree  of  under¬ 
expansion  than  did  Figure  12.  Figure  13  corresponds  to  flow  from  a 
convergent  nozzle  with  sonic  velocity  at  the  exit.  The  pressure  ratio  is 
PT/P  ,  ■  2.1.  This  figure  again  shows  radial  variations  of  both  total  and 
static  pressures  at  each  of  two  axial  stations,  x/R  **  1.90  and  x/R  ■  2.65. 

The  rate  at  which  the  mixing  region  spreads  with  distance  from  the  nozzle  exit 
can  again  be  seen,  as  can  the  increasing  total  pressure  loss  due  to  shocks. 

Note  the  relatively  large  levels  of  static  pressure  variation  even  though  the 
underexpansion  is  still  mild.  Finally,  note  that  the  viscous  boundary  layer 
analysis  by  itself  would  predict  the  same  flow  field  for  all  three  jets  in 
Figures  12  and  13,  (assuming  the  impressed  pressure  were  taken  as  the 
ambient  pressure  in  all  cases).  Also  note  that  the  considerable  effect  of  the 
mixing-induced  total  pressure  loss  on  the  flow  field  would  be  Ignored  by  pure 
inviscid  analyses. 

The  last  two  figures  show  the  predicted  shock  wave  shapes  for  the  jets 
of  Figure  12  and  13.  Figure  14  shows  a  comparison  between  the  shock 
shape  which  is  predicted  by  the  complete  inner-outer  analysis  and  the  shape 
predicted  by  a  completely  inviscid  analysis.  The  inviscid  calculation  was  made 
by  specifying  the  entropy  to  be  constant  along  all  streamlines  (except  for  shock 
losses).  As  Figure  14  shows,  the  two  calculations  give  nearly  identical 
shock  shapes  before  and  immediately  after  the  reflection  of  the  shock  from  the 
centerline.  Indeed,  the  minor  differences  between  the  two  calculations  in  this 
region  is  more  due  to  small  errors  in  calculation  (stemming  mostly  from  undesired 
interactions  between  the  initial  part  of  the  mixing  layer  and  the  expansion 
fan  at  the  nozzle  lip)  than  from  the  physics  of  the  problem.  However,  the 
sharply  curved  portion  of  the  shock  near  the  outer  boundary  is  due  to  real 
effects.  This  sharp  curvature  comes  about  as  the  shock  enters  the  strongly 
rotational  flow  region  which  has  been  created  by  the  viscous  mixing.  As  the 
shock  traverses  this  mixing  layer,  the  Mach  number  in  front  of  the  shock 
approaches  unity  so  that  even  as  the  shock  turns  normal  to  the  flow,  its 
strength  decreases  until  it  eventually  fades  out. 
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Also  shown  in  Figure  14  are  the  outer  boundaries  of  the  inviscid  calcu¬ 
lation  (whiih  is  of  course  a  streamline)  and  the  "sonic"  matching  line  (M  =  1.1) 
which  was  used  in  the  coupled  analysis.  Finally,  the  mixing  region  is  shewn 
by  the  shaded  area. 

Figure  15  presents  a  final  comparison  with  experiment.  This  figure 
again  shows  the  predictions  of  the  coupled  analysis  (for  a  jet  of  pressure 
ratio,  Pj/Pamb  =  2.1),  this  time  superimposed  on  a  Schlieren  photograph  of  a 
jet  at  the  same  conditions.  As  can  be  seen,  the  agreement  is  excellent.  Note 
that  the  curved  portion  of  the  shock  in  the  mixing  region  agrees  quite  well 
with  the  schlieren  result,  and  that  the  predicted  size  of  the  "Mach  disc" 
agrees  with  the  photograph  (although  without  the  theoretical  prediction 
superimposed,  the  photograph  appears  to  show  a  regular  reflection).  Again, 
note  that  although  the  predicted  shock  starts  much  too  close  to  the  nozzle 
exit,  it  remains  extremely  weak  until  it  nears  the  axis  and  so  would  not  be 
expected  to  be  visible  on  the  Schlieren  photograph.  Finally  Figure  15  shows 
that  both  the  predicted  and  the  experimental  outer  boundaries  show  a  point  of 
inflection  at  about  the  axial  distance  from  the  nozzle  exit  where  the  shock 
reflects  from  the  centerline.  This  inflection  in  the  outer  boundary  is  caused 
by  the  displacement  of  the  viscous  mixing  region  by  the  inner  inviscid  core 
of  the  plume.  At  the  exit  plane,  the  inviscid  flow  turns  outward  through  an 
expansion  fan.  Then  the  axisymmetric  effects  force  this  flow  to  again  turn 
and  approach  the  axis  (see  outer  boundary  shape  of  the  inviscid  calculations 
in  Figure  7.  The  superposition  of  an  ever-widening  mixing  region  on  these 
curved  inviscid  streamlines  generates  the  inflection  in  the  boundary. 
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Figure  1  The  Flovfleld  of  an  In  viscid  Two-dimensional  Supersonic  Jet 
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Figure  2  The  Flovfield  of  an  Ideally  Expanded  Viscous  Jet 


Figure  3  Subdivision  of  the  Jet  into  Inner  and  Outer  Regions 
(Outer  Region  Shown  Shaded) 
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Figure  8  Intersection  of  Shock  Shape  with  Axis 


of  Symmetry 


Figure  9  Height  of  Mach  Disk 


PREDICTED  SHOCK  SHAPF. 


Figure  11  Cross-Strear  Variation  of  Predicted  Total  Pressure  and  Static  Pressure 

for  Ideally  Expanded  Jet  M£xlt  =  1.60,  Pamb/PTref  =  *242 
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Figure  13  Cross-Stream  Variation  of  Predicted  Total  Pressure  and  Static  Pressure 
for  Jet  from  Convergent  Nozzle,  Showing  Total  Pressure  Loss  Due  to  Both 

Shocks  and  Mixing.  M  =  1.0,  P  /P  =  .242 
exit  ’  amb  Tref 
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Figure  lL 


Comparison  Between  Predicted  Shock  Shapes  Using  Inviscid  Prediction 
Technique  and  Coupled  (Inner-Outer  Analysis)  Viscous  Technique. 
Shaded  Region  Represents  Flow  Which  Is  Signifigaady  Affected  by 
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II.  THEORETICAL  ACOUSTIC  MODELS  OF  TURBULENT  JETS 


1.  MOVING  SOURCE  MODELS  FOR  JET  NOISE 
R.  Man! 


The  present  section  considers  several  problems  of  the  sound  power,  power 
spectrum  and  directivity  produced  by  moving  acoustic  sources  shrouded  by  jet 
flows.  The  jets  are  assumed  (for  simplicity)  to  be  characterized  by  a  slug 
flow  or  top  hat  type  mean  velocity  profiles  in  most  cases.  The  sources  are 
simple  harmonic  in  their  wwn  frame  of  reference  and  are  assumed  to  convect 
with  the  same  velocity  as  the  jet  as  well  as  at  velocities  of  about  0.65  the 
jet  velocity. 

The  studies  are  all  motivated  by  one  notion,  namely,  that  Llghthill's 
original  idea  of  ascribing  jet  noise  to  convected  sources  radiating  freely  to 
the  ambient  needs  revision  to  allow  for  mean  flow  "shrouding"  effects.  The 
studies  explain  several  experimentally  observed  features  of  jet  noise  such  as 
the  failure  to  exhibit  convective  amplification  (particularly  at  high  frequencies 
and  shallow  angles  to  the  exhaust  axis)  and  associated  failure  of  peak 
frequencies  in  the  power  spectrum  to  shift  linearly  with  jet  velocity.  Impli¬ 
cations  for  the  jet  density  exponent  issue  for  heated  jets  are  also  considered. 
The  study  may  be  regared  as  moving  source  solutions  to  the  Phillips*  equation 
for  jet  noise  with  a  specific  velocity  profile,  namely  the  top  hat  profile. 

The  advantage  of  choice  of  a  simple  velocity  profile  is  to  obtain  solutions 
valid  for  arbitrary  frequencies. 


1.1  First  Model  Problem  (Figure  16) 


Consider  the  problem  of  determining  the  sound  field  due  to  a  fluctuating 
monopole  point  source  translating  at  a  uniform  subsonic  velocity  Me  (where 
M  <  1,  M  being  the  Mach  number  and  c  is  the  speed  of  sound).  The  source 
translates  along  the  axis  of  a  round  jet  whose  velocity  profile  we  assume  to  be 
a  slug/flow  velocity  profile.  Also,  the  jet  velocity  is  taken  equal  to  that 
of  the  source.  The  problem  is  illustrated  in  Figure  16.  The  monopole  source 
is  assumed  to  have  a  time  dependence  in  its  own  frame  of  reference  of  q0  cos 
(u0t).  The  mean  jet  density  and  temperature  are  assumed  to  be  the  same  as  that 
of  the  ambient. 
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«.  a  i  i  >  itfjgti 


Analytically,  we  wish  to  determine  an  acoustic  velocity  potential  p  which 
satisfies  in  region  I  (outside  the  jet) 


0 


(1) 


and  in  region  II  (within  the  jet) 


(i  -  M*)*xx  ♦  V’* 


—  ccs(w  t)5(x 
PC  ° 


Mct)6  (y)4(z) , 


(2) 


where  V ^  stands  for  the  Laplace  operator  in  the  y  -  z  plane.  At  the  jet 
still-air  interface,  i.e.,  at  r  =  a,  we  require  (a)  continuity  of  pressure, 
p,  where 


p  =  pQ 


t’ 


in  region  I, 


(3) 


and 


p  =  -p  (i  +  Mcj>  ),  in  region  II,  (4) 

O  L  X 

and  (b)  continuity  of  radial  acoustic  particle  displacement,  say  o,  where 
4>r  =  nt,  in  region  I,  (5) 

and 

< ; i  =  nt  +  Mcnx*  in  regi°n  II.  (6) 

An  elegant  procedure  of  solution  suited  to  the  above  problem  has  been  given 
2 

by  Morse  and  Ingard  and  we  follow  closely  their  method  of  solution. 

Let  $,  p,  etc.,  denote  the  Fourier  transforms  with  respect  to  time  of  the 
corresponding  physical  quantities.  Thus 

«• 

i  ■  JT  I  >  Jt  i  -  ™  (7) 

-  m 


and 
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W  -  AU'-LWI'  U.HMIU. 


Also,  we  write  cos(u)Qt)  =  — [exp (ju^t )  +  exp(-j(jQt)  ]  .  The  problem  for  the 
transforms  is 


V*£  ♦  kl$  ■  0, 

q«(y)«d) 

(1  -  Ml)$xx  ♦  Vs,*  -  2jk>l»x  *  kl«  -  — TSp-TT — 


in  region  I, 


d') 


j(k  -  k  )  x 

x  [exp ( - jjj — 2 — )  ♦  exp( 


j ♦  k  )x 


-)]. 


in  region  II; 


(2’) 


P  °  j<“P0$  . 


in  region  I, 


O') 


P  *  -P0(-jwi  ♦  Mc$x)  , 


in  region  II; 


(A') 


*r  ■  *jwn  , 

$r  -  -jwn  +  Men  , 


in  region  I, 


in  region  II. 


(5') 


(6’) 


Let  41  ~  <t>+  +  <j>  and  similarly  for  p  and  h  where  ij)+  corresponds  to  the  solution 
with  the  term  exp(j(k  -  k0)x/M)  in  equation  (2’)  and  <j>  to  the  term  involving 
exp  {j  (k  +  kQ)x/M}.  Note  that  k  =  u/c,  kQ  =  u)Q/ c,  etc. 

Consider  in  detail  the  problem  for  Intuitively,  it  is  clear  that 

<J>+,  p+,  h+  all  have  an  x-dependence  of  the  type  exp{j(k  -  kQ)x/M}.  "Factoring” 
this  dependence  out,  one  is  left  with  the  following  problem  in  the  y  -  z  plane: 


vjr  ♦  k-r  -  0, 


in  region  I; 


(8) 
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,.  ....  qn'SCyJ-sCz) 

VU  *  <  +  V  -  -2-| - n - 

*  4  ripple 


in  region  II; 


where 


k+*  ■  c—r-HCr-^i-^ck-i-^,)] 


and  is  >  0  only  if  k*  >  k  >  k ' ,  where  k+  »  k  /,!  -  M\  and  k"  *  k  /£  *  Nl^ 
—  o  —  —  o  oo\__/  oo  'O  ~s 


ic*2  *  — t  ((«*  -  x)(k  -  «")), 
o  o 


where 


k  =  k  (1  +  M),  k  =  k  (1— m) 
o  o  o  o 


Also,  let 


-  jJT  I(k  *  <)(k  ■ 


Note  that 

,  ♦  ♦  .  ,  - 
X  >  K  >  t._  >  k  >  K  . 
o  —  O  ~  O  —  O  —  0 

•f  2  +  “ 

The  fact  that  k  >_  0  only  if  k^  >_  k  >_  kQ  expresses  the  result  that  in  the  far 
field  the  moving  source  yields  a  frequency  spectrum  containing  frequencies  in 
the  range  <Wl  -  M  >  w  >  u>  / 1  +  M  which  is  what  we  expect  from  the  Doppler 
shift  formula.  We  restrict  our  attention  to  this  range  of  k.  The  matching 
conditions  for  equations  (8)  and  (9)  are  that,  at  r  =  a. 


(pressure  matching  conditions),  (13) 


$(r»a)»r-  5  ('  “  3  ) 
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♦  r  (r  "  '■>*)  ■  5*  (r  »  a ') 

o 


(transverse  particle  dispace-  (14) 
ment  matching  condition). 


To  solve  equations  (8),  (9),  (10),  and  (11)  in  the  range  k  <  k  <  k  and  with 
restriction  to  outgoing  waves  at  infinity,  in  the  range  kQ  k  <_  kQ  we  assume, 
for  in  regions  I  and  II, 


in  I: 


(15) 


-♦  ♦  m  ♦  Kk  -  k0)x 

♦  ■  A;  H^(k*r)exp  -  0 


IT 


in  II: 


jq0lli1J(<r). 


«  -  [A„  J0(<  r)  -  — r - ]e*P 

and  if  k+  <  k  <  k+,  in  region  II, 

O  —  o 


j(k  -  k0)x 


(16) 


*  ♦  ♦  %Ko^r),  Hk  *  Vx 

♦*  -  [At .  I0(£'r)  -  exP 


8a2pQMc 


t r 


(17) 
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(The  form  for  <>  in  region  I  is  independent  of  whether  k  ^  £  kq.) 

+2 

Note  that  the  change  of  sign  k  depending  on  whether  k  e[kQ(l-M0,  kQ (1— M)  ] 
is  associated  with  the  fact  that  if  the  jet  in  the  present  problem  were  of 
infinite  ra  lius  (i.e.,  the  moving  fluid  occupied  all  space)  the  Doppler  shifted 
frequencies  would  range  over  wq(1  -  M)  to  wq(1  +  M) .  In  other  words,  as  is 
well  known,  there  is  a  difference  in  the  Doppler  shift  frequencies  depending  on 
whether  the  observer  moves  towards  a  source  or  whether  the  source  moves  towards 
the  observer.  This  difference  will  be  seen  later  to  play  a  key  role  in 
suppressing  convective  amplification  at  high  frequencies. 

Equations  (15)  and  (16)  or  (15)  and  (17)  may  now  be  readily  solved  for 
A+  and  by  using  the  matching  conditions  (13)  and  (14)  .  Since  we  are 
interested  in  far  field  pressures  far  outside  the  jet,  we  only  give  the  result 
for  A*: 
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(a)  if 


+ 

'  > 

o 


*  qck  k0<+[Y0(.-4a)J1(,c*a)  -  Y1(<*a)Jo(<*a}) _ 

1  16»p0V.c[k1<*H^1JCk*a)J1(<*a)  -  k*k*J0C** s  (k*a)  ]  ' 


(18a) 


(b)  if 


+ 


<• 


-qk<  ko [ Ko  c<  a)Ij(e  a)  ♦  IQ(ic  a)Kj  a)] 


8*,P0V[k*i*i1(«*a)Hi1)  Va)  ♦  k*k* I0 (<* a) hJ1 1 (k*a) ] 


IT 


(18b) 


Equation  (18)  essentially  completes  the  formal  solution  to  the  problem.  The 
■ ar  field  pressure  and  the  radial  acoustic  velocity  may  be  computed  by  using 
p  =  jwpQC>  and  In  this  problem,  every  point  on  a  cylindrical  surface 

concentric  with  the  jet  experiences  the  same  pressure  time  history.  Morse 

2 

and  Ingard  have  discussed  thoroughly  the  problem  of  determining  the  power 
spectrum  and  total  power  radiated  by  the  source  and  their  concluding  result  is 
that  the  power  spectral  density  extends  over  a  frequency  range  [m  /(I  +  M)  ] 

<  u)  <  [iW(l  -  M)]  and  is  given  by 


(16irpoMcw)  |  Aj  | 1  -  I  (u)  . 


(19) 


The  total  power  is  given  by 


Wq/U  -  M) 

I  (u)  doi . 

u0/(l  *  M) 


(20) 
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Acutally  Morse  and  Ingard  consider  th-  case  of  a  monopole  point  source 
convecting  at  Me  in  free  space,  for  which  case 


A 


I 


(21) 


and  hence 


1(«) 


for 


r-r-R-  i  u  i  T-TT7 


and  the  total  power  is 


(22) 


P  - 


q,ui* 

Mo  o 


8itp0(l  -  M* )  2 C 


(23) 


Thus,  in  thj  case  of  a  convected  monopole,  the  convective  amplification  is 
as  (1  -  M2)"2. 

If  we  take  the  limit  as  kQa  -  0  of  equation  (18),  we  find  that  A*  tends 
to  (independent  of  whether  k  >  <+  or  k  <  k+) 


At  - 


-5%“ 


I  16np  Mclj  ■ 
o  n 


(24) 


so  that 


«ou 


'  IsfpTC^5, 


for 


o  .  .  o 

-nr  i  “  -  T~nr 


(25) 


and 


P 


q2u2  (1  ♦  M2) 

^0  0 _ 

8p0c:r(l  -  M2)' 


(26) 
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In  the  general  case,  A^,  i(^i)  and  P  are  given  by  equations  (18a),  (18b),  (19) 
and  (20),  a  J  specific  results  will  be  discussed  in  the  following. 

2  2 

The  t l  .  power  emitted  bv  such  a  source  nondi  „ensionalized  by  [q  u  /8trp 

~  O  O  0 


(1  -  M-)*Y1  and  exp 
Figure  17. 


ad  in  dB  is  plotted  as  a  function  of  (kQa)  and  M  in 


Shown  by  single  points  on  the  extreme  right  in  Figure  17  are  points  given 

2 

by  20  log  (I  -  M  ),  being  the  correction  if  there  were  no  convective  ampli- 

1  3 

fication  a.  ail  corresponding  to  Ssanady'.  suggestion.  The  portions  of  the 
cur  es  corresponding  to  corrections  >  odB  indicate  underestimates  of  convective 
amplification  as  estimated  from  a  freely  moving  source  model  and  conversely. 

Clearly,  such  curves  confirm  the  frequency  dependent  nature  of  convective 
amplification.  The  curves  flatten  as  we  move  to  Lhe  right  and  if  we  identify 
the  point  on  each  curve  (for  the  different  Mach  numbers)  at  which  the  correctio 
is  wrlthin  a  decibel  of  the  limit  as  (l^a)  *•  one  deduces  that  beyond  a  source 
Strounal  number  I(2f0a)/Mc]  of  D.5  there  would  be  no  significant  convective 

4 

..rp] if ication.  Figure  15  of  Lush’s  paper  indicates  lack  of  convective 
amplification  beyond  [2f0a/Mc]  of  about  0.3. 

Finally,  we  consider  the  implications  with  regard  to  Strouhal  scaling  of 
the  results  shown  in  Figure  17.  As  a  starting  point,  in  Figure  18  we  show 

under  the  curve  labelled  M  =  0.3,  one-third  octave  intensities  obtained  by 

4  0 

Lush  in  Figure  8  of  his  paper  for  a  jet  Mach  number  of  0.37  at  90  .  This 

curve  is  chosen  as  a  base  line  because  at  that  low  Mach  number  of  0.37  and 
location  (90°  to  jet  axis)  we  expect  little  convective  amplification  effects, 
’ibe  abc  issae  are  shown  in  Strouhal  numbers,  St  =  (2fa/Mc),  and  the  ordinates 
are  only  relative  decibel  levels. 

An  intensity  spectrum  at  90°  was  chosen  because,  in  addition  to  lack  of 
convective  amplification  effects,  the  90°  location  also  provides  a  very  good 
and  clean  measure  of  the  intrinsic  strength  of  the  sources  (their  frequency 
distribution).  This  is  because  that  location  is  largely  characterized  by 
"self  noise".  A  basic  assumption  of  the  process  used  in  deriving  Figure  19  is 
that  the  frequency  distribution  of  the  "intrinsic  source  strengths"  does  follow 
Strouhal  scaling  with  respect  to  velocity.  This  is,  of  course,  excellently 
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borne  out  by  Figure  8  of  reference  £ 4 J  where  Lush  shows  that,  at  the  90° 
location,  Strouhal  scaling  with  respect  to  velocity  was  obtained.  The  basic 
argument  of  what  follows  is  to  point  out  that  the  radiative  efficiency  of  the 
sources  is  frequency  dependent  and,  being  higher  for  the  low  frequencies  than 
for  the  high  frequencies,  causes  peak  frequencies  of  the  sound  power  spectrum 
to  scale  with  velocity  much  slower  than  a  first  power  (as  is  assumed  in 
conventional  Strouhal  scaling).  The  particular  low  Mach  number  datum  used 
to  establish  this  result  (taken  in  this  case  as  the  90°  intensity  spectrum  of 
Lush^)  is  not  the  main  issue  of  this  paper:  a  different  datum  would  lead 
to  the  same  qualitative  conclusions.  Ideally,  perhaps,  one  would  have  to  work 
out  spearately  the  "shear  noise  i nd  "self  noise"  portions  of  the  power  spectra. 

The  spreading  of  the  source  frequency  due  to  the  Doppler  shift  makes  it  a 
little  difficult  to  apply  Figure  17  directly.  However,  it  can  be  showr.  that 
the  Doppler  spreading  will  be  narrower  than  conventional  moving  source  results 
would  indicate  (see  Reference  5).  Further,  if  ve  are  interested  in  the  sound 
power  spectrum,  it  seems  reasonable  to  apply  Figure  17  to  Figure  18  as  follows. 
For  each  Strouhal  eumber  St,  and  Mach  number  M,  determine  a  source  frequency 
parameter  l^a  =  St  *  TO  and  then  determine  the  decibel  correction  from  Figure 
17.  Starting  with  the  curve  labeled  M  *  0.3,  such  a  frequency  dependent 
correction  procedure  was  applied  to  derive  the  curve  labelled  M  =  0.5,  M  =  0.7 
and  M  =  0.9  from  the  curve  labelled  M  ■  0.3.  As  expected,  one  observes  a  shift 
back  of  the  peak  frequency  (in  terms  of  the  Strouhal  numbers)  at  which  the  sound 
power  spectrum  peaks.  The  spectra  are  pretty  flat  as  is  typical  of  jet  noise 
but  an  attempt  was  made  to  estimate  the  peak  Strouhal  number  as  a  function  of 
Jet  Mach  number  and  the  results  are  shown  in  Figure  19.  Undoubtedly  by  a 
purely  fortuitous  coincidence,  the  curve  in  Figure  19  is  fitted  very  well  by 
a  relation  of  the  type  (St)p  =  (0.21)/M.  Since  the  Strouhal  number  itself  is 
given  by  (fpD/V) ,  Figure  II-5  suggests  that  the  peak  frequency  in  the  sound 
power  spectrum  is  independent  of  jet  velocity  being  given  (in  the  case  of 
Figure  19  by  [(0.21)c/Dj.  Such  a  tendency  for  the  peak  frequency  to  be 
independent  of  Jet  velocity  has  been  noticed  in  several  experiments. 
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The  suggestion  that  emerges  therefore  is  that  the  tenacious  adherence  of 
the  total  power  to  an  eigth  power  law  as  well  as  the  tendency  of  peak  frequency 
of  the  powc  r  spectrum  to  be  relatively  insensitive  to  jet  velocity  are  both 
manifestations,  of  the  same  result  indicated  by  Figure  17,  namely  the  inhibition 
of  convective  amplification  with  increasing  frequency  and  jet  velocity. 


1 . 2  Asymmetric  Line  of  Source  Convection 

To  ii  '  .r  the  effect  of  lines  of  source  convection  different  from  the  jet 
axis,  we  study  the  acoustic  output  of  a  line  acoustic  source  convecting  at  the 
jet  velocity  in  a  plane,  slug  flow  jet.  The  problem  is  two  dimensional  and, 
as  indicated  in  Figure  20,  the  line  source  is  allowed  to  convect  at  a  distance 
ch  from  the  jet  centerline  where  h  is  the  half  width  of  the  jet.  We  restrict 
o  to  0  <_  a  <  1  to  ensure  that  the  source  is  always  within  the  jet.  The  case 
of  a  =  0  obviously  corresponds  to  case  of  symmetric  or  centerline  source 
convection. 

The  analytical  problem  corresponding  to  Figure  20  may  be  described  in  terms 
of  the  acoustic  velocity  potential  4  as  follows.  We  wish  to  determine  <f>  such 
that  in  regions  I,  III  (outside  the  jet),  $  satisfies: 


r 


(27) 


Inside  the  jet  (in  region  II),  it  satisfies: 


(i 


i‘.Ct ) 


x  5  (y  -  n'.i)  cos  (u  L'  . 

o 


(28) 


At  both  interfaces  (y  =  +  h)  the  acoustic  pressures  and  transverse  particle 
displacements  must  be  continuous.  Also  in  regions  I  and  III  only  outgoing 
waves  are  permitted. 


*-  :  3-Vo 
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As  in  Subsection  1.1,  the  method  of  Fourier  transforms  will  be  applied 
since  the  problem  is  basically  a  transient  problem. 

Let  <j>  be  the  Fourier  Transform  of  <f,,  i.t- 

00 

♦  =  sr  |  ♦  oj,Jt  at  (j  .=  /=Z) 


where 


* 


w 

-/ 


i  du 


Also  we  write  cos<wot)  =  A  [exp^t)  +  e*p(-  j^t)  ) 
The  problem  for  <j>  is: 

,2  2  2  - 

^  ^  +k  ♦  =  0  in  regions  I  and  III; 


(29) 


(27*) 


a  -  m!Hxx  +  iyy  +  2ikt:Jx  + 


qo  6  (y  -  oh) 

4tt  p  Me 

o 


(exp  (- 


j  (k  -  k  )  x 

- o_: 

M  ‘ 


j  ( k  +  k  )  x 

+  exp  ( - 0 

M 


)  }  in 


"eg ion  II. 


(28*) 


The  transforms  p  and  n  of  the  acoustic  pressure  and  transverse  acoustic 
particle  displacement  are  related  to  $  by: 


P  in  I,  III 


-  -p  ^  in  II 


(30) 
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<*>  -  - j a: p.  in  I,  ill 


«?  ,  =  ’  3->;i 

y 


:  Kcr  in  II 
x 


Let  J  =  (|,+  +  (♦,“  and  similarly  for  p  and  rj  where  $+  corresponds  to  the  solution 
with  the  texn  exp  {-J  ~^k°^x}  in  (28*)  and  $  to  the  term  exp{^  in 

(28').  Note  that  k  =  cj/c0>  kQ  =  w  /cq,  etc.  Consider  in  detail  the  problem  for 
I+.  Intuitively  it  is  clear  that  $+,  h+  and  p+  all  have  an  x-dependence  of  type 
exp { j (k  -  k  ) x/M} .  "Factoring"  this  dependence  out,  one  is  left  with  the  following 

^  — f- 

one-dimensional  problem  for  $  : 


d23  + 


+  k+2?>+  =  0  in  I,  III 


where 


q  6  (y  -  oh) 

M-  +  k  +  2a+  =  - - - -  in  II 

dy^  9  4 up  he 


,4-7  (1  -  M  )  t  U _ v\  IV  -  - - - 1 

*  =  - -  1  (X  -  M  KMK  1  +  M 


is  >  0  only  if  k+  >  k  >  k  ,  where  k+  =  k  /  (1  -  M)  and  k  =  k  /  (1  +  M) 
—  o  —  —  O  O  o  oo 


k+!  ■  JF  <<<  -  k)[k  -  ‘o’1 


where  x  =  k  (1  +  M)  and  x  =  k  (1  -  M) .  Also,  let 
o  o  oo 


.  +  2  =-  _k 


W  1(k  -  co)(k 


*o>! 


IWgP'WWWflU'.' 


.  'I'wwn 


t«  ~trj4vzv-r#*vitr-i<>  ■</  *• 


| 

| 


Note  that  k+  >  k+  >  k  >k  >k. 

o  —  o  —  o  —  o  —  o 


The  fact  that  k+2  >  0  only  if  k  >  k  >  K 
—  o  —  —  o 


expresses  the  result  that  in  the  far  field  the  moving  source  yields  a  frequency 
spectrum  containing  frequencies  in  the  range  u>q/  (1  -  M)  >_  w  >.  u)q/(1  +  M)  which 
is  what  we  expect  from  the  Doppler  shift  formula.  We  restrict  our  attention 
to  this  range  of  k.  The  matching  conditions  for  equations  (32,  33)  are: 


and 


$+{y  =  t  h  outside  jet} 


(y  =  ±  h  inside  jet)... 


(37) 


±  h  outside  jet) 


k_  d$* 
k  dy 


(y  =  ±  h  inside  jet)... 


(38) 


To  solve  (32),  (33)  in  the  range  kQ  <_  k  k+  subject  to  the  matching  conditions 
(37),  (38)  is  actually  a  somewhat  more  tedious  problem  than  the  problem  studied 
in  Section  1.1.  Assume  for  <}rt-  the  following  forms  in  regions  I,  II  and  III: 


j(k  -  k  )x 

A  exp { - - - }  exp [ jk  (y 


h)]... 


(39) 


In  region  if  k  <  k  <  k  : 

o  —  —  o 


{B  sin[tc+(y  -  h) } 

+  C  sin  [k+  (y  +  h) } 


j  qQ  exp[j 

8tt  pq 

j(k  -  k 
exp[ - rr— 


<+|y  -  ah | ] 
Me  K  + 


(AO) 
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+  ,  + 

on  the  other  hand  if  k  <  k  <_  k  ; 


4>+  ~  (B  sinh[ir+(y  -  h)3  -  — — 

8tt  p  Me  K  + 

o  — 

+  j  (k  -  k  )x 

+  C  smh[<  (y  +  h )  ]  }e>:p  [ - ~ — 2 — j 


q  exp  [-ic+  |y  -  oh  M 


(40b) 


In  region  II: 


<>  =  D  expj- 


j  (k  -  kjx 

~M 


-]exp[-j  k  (y  +  h)] 


(.1) 


Evaluating  and  —  in  regions  I,  II  and  III  and  applying  the  matching 
conditions  (37,  38)  yields  a  set  of  four  simultaneous  equations  for  A,  B,  C.  D. 
Solving  for  them  yields  (for  A,  D):  if  kQ  <_  k  kq: 


(A,  D) 


k  k  cos  (k+  oh) 

{ - ° - 

[k2  tc+  sin(<+  h)  +  j  k*  k+  cos  (i:+  h)  ] 


k  k  sin(<  oh) 
o 


■)  q° 


[k2  k+  cos  (k+  h)  -  j  k!  k+  sin(<+  h )  ]  071  po  Mc 


(plus  sign  goes  with  D,  minus  sign  with  A) 


and  if  <+  <  k  <  k*: 

o  —  —  o 


(42) 
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A,  D) 


-q 


k  k  cosh  (k  oh) 

Or  o  — 

Me  , ,  5  _+ 


8tt  pq  Me  j j, 2  <  +  sinh(^+  h)  -  j  k+  cosh  (£+  h)  ] 


sinhU  oh)k  k 


+ 


{k  £  cosh(<  h)  -  j  k  k^  sinh  (<_  h)] 


(43) 


The  power  radiated  across  tvo  planes  shown  dotted  in  Figure  20  may  be  obtained 
from  the  solutions  (42,  43)  by  integrating  the  power  spectral  density  which  is 


4tt  pq  V  k+  u>  [  | A |  2  +  jD|  2  ] 


(44) 


((44)  may  be  derived  by  a  calculation  very  similar  to  that  outlined  in 
Reference  3)  over  a  frequency  range  u)Q  (1  +  M)  ^  <_  u>  <_  io0(l  -  M)  .  If  there 


were  no  jet  at  all,  one  would  determine  the  acoustic  field  by  solving  for 


d2P  + 


dy 


+  k+2  tp+  - 


Sn  pc  V 


<S  (y)  .  .  . 


(45) 


the  solution  to  which  is: 


-+  -j  q0  e*?tj  **ly|l 

V  “ - z - 


(46) 


8ir  p  V  k 
o 


48 


The  power  in  this  case  is: 

U 


1  -  M 


OJ  Q-J 


a  w 
*o  o 


8r  p  V 
o 


0) 


ep  a  -  m2  ) 

o 


2,372 


1  +  H 


(47) 


2  -3/2 

Thus  cne  conv  ction  amplification  factor  for  line  velocity  sources  is  (1  -  M  )  . 

2  -2 

(The  corre  ;ponding  resuit  for  point  velocity  sources  is,  of  course,  (1  -  M  )  .) 

In  the  limit  as  kQh  -*•  o,  (42)  and  (43)  give 


-j  qn  k 

A  or  D  =  - 2 _ 


87i  p  k  k 

o  o 


(48) 


As  in  Section  1.1,  the  moving  source  in  a  jet  of  limitingly  small  thickness  does 
not  reduce  to  a  freely  moving  source  and  one  finds  that  the  convection  amplifi¬ 
cation  factor  for  a  line  displacement  source  (following  the  terminology  of 
Morfey  and  Tanna^  is  (1  +|m2)  (1  -  M2)  ^2,  i.e.,  stronger  than  the 
freely  moving  line  source  (or  line  velocity  source)  by  a  factor  (1  +  yM2) 

(1  -  M  )  .  (The  corresponding  factor  for  point  sources  is  (1  +  M  ) 

2  -2 

(1  -  H)  ) . 

More  general  calculations  require  incorporation  of  (42)  or  (43)  in  (44) 

and  an  integration  over  the  frequency  range  u0(l  +  M)  *  m  <_  wQ(l  -  M)  ^  to 

obtain  the  total  power.  Such  calculations  have  been  carried  out  for  M  ranging 

from  0.5  to  0.9,  (kQh)  ranging  from  0.01  to  1.0  and  a  from  0  to  0.99.  In  all 

cases  (results  are  shown  in  Figure  21,  the  total  power  Is  normalized  by  (47), 

2  2  3/2 

i.e.,  by  [qQ  w0/8p0(l  -  M  )  ].  The  interest  in  Figure  21  is  really  in  how 

much  the  results  in  the  range  0  <  a  <  1  deviate  from  the  a  =  0  results  (case 
of  centerline  source  convection).  As  is  observed  from  Figure  21,  the  results 
for  total  power  are  essentially  unaffected  by  the  exrct  location  of  the  line 
of  source  convection  so  long  as  it  Is  within  the  jet.  This  result,  while 
deduced  for  a  plane  problem,  lends  confidence  to  the  notion  that  the  results  of 
Cectirn  1.1  will  have  a  wider  range  of  applicability  than  just  to  sources  con- 
vecting  along  the  jet  centerline.  The  results  of  Figure  21  are  physically 
tantamount  to  the  observation  that  it  is  the  "total  extent"  of  shrouding  to 
which  the  moving  source  is  exposed  which  determines  it  radiative  efficiency. 
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1.3  Effect  of  Finite  Shear  on  Model  Problem  Studies  In  Section  1.1 


I 


In  a  recent  study  Eversman^  has  Indicated  an  Interesting  procedure  (in 
connection  wit^i  lined  duct  propagation  problems)  by  which,  in  case  the  transi¬ 
tion  from  the  mean  velocity  to  the  ambient  value  of  zero  velocity  occurs  in  a 
shear  layer  of  thickness  6  (the  slug  flow  model  assumes  this  transition  to  take 
place  in  a  shear  layer  of  zero  thickness) ,  calculations  correct  to  order 
(6/A)  or  (6/a)  (A  ■  wavelength  of  sound,  a  ■  jet  radius)  may  be  carried  out. 
This  section  will  briefly  summarize  the  application  of  Eversman's  approach  to 
the  model  problem  of  Section  1.1.  To  set  the  stage  of  how  Eversman's  method 

needs  to  be  adapted  for  the  current  study,  the  application  to  the  problem  of 

8  9 

Section  1.1  is  prefaced  by  an  application  to  the  Ribner-Mlles  ’  problem  of 
reflection  of  plane,  acoustic  waves  by  a  velocity  discontinuity. 

In  terms  of  Figure  22.  the  problem  is  to  determine  the  reflection 
coefficient  of  a  plane,  simple  harmonic  wave  by  a  shear  layer  of  thickness  6 
across  which  the  velocity  changes  from  a  value  MpC  to  zero. 

Eversman  used  the  terminology  of  inner  and  outer  expansions  to  develop 

his  method,  but  (in  hindsightl)  his  approach  can  be  illustrated  without 

recourse  to  such  terminology  as  under.  In  general,  the  governing  equation  for 

the  linear,  lnviscid  propagation  of  sound  in  a  parallel  sheared  flow  involves 

a  third  order  differential  equation  but,  in  terms  of  Figure  22,  if  interest  Is 

restricted  to  pressure  waves  of  type  P(v)exp[j  (kKx  -  tot)]  where  k  ■  —  (this 

c 

class  of  waves  is  all  that  is  involved  in  the  Ribner-Miles  problem) ,  we  may 
show  that  P(y)  is  governed  by: 


d  <• _ 1 _  dP,  _  f  K2 

dy  (1  -  M(y  )K)‘r  dy;  (1  -  M  (y )  K> 1 


-  1}  k2P 


Now  integrate  (49)  form  y  *  0  to  y  ■  5  to  obtain  first: 


( - - - )  2  — ] 

'1  -  M  K'  dvJ 


dP] 

dy 


y  =  6  y  =  0 

6 

k  I  *[i  -  m(v)k]z  "  p(y)dy 
0 


(49) 


(50) 
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As  it  stands  (50)  is  just  an  integro-dif ferential  version  of  (49)  but  (antici 
pacing  a  step  to  follow  in  (52)  by  integrating  (22)  twice)  to  zeroth  order  in 
O/A)  P(y)  =  P(0)  or  P (6 )  so  long  as  0  <  y  5  6.  Correct  to  order  O/A)  then 

we  may  assert  that: 


■ft  -  V>‘  dy  v  =  6 


(dP) 

Vdy ' 


=  0 


+  k2  [P  (0)  or  P  (5)  ]  I-j, 


(51) 


where 


h  r 


Wi’.vftrT-r  -  U«*y 


By 


integrating  (50)  twice,  again  to  order  O/A),  one  finds  that: 


P  (y  =  <5 ) 


P(y 


0)  + 


(52) 


where 


6 

’ 

(1  -  M(y)K)  2  dy 

« 

0 
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It  is  rather  important  to  observe  that  in  (52)  we  may  substitute  (if  we 


choose)  (^)  ^  by  (^) 


- - - —  and  not  by  (j^) 

6  (1  -  M  Kr  y  y  =  5 

o 


In  other 


words  to  zeroth  order  in  (6/A),  the  quantitites  conserved  across  the  shear  layer 

1  dP 

are  P(y)  and  - 9  .  As  Eversman  has  pointed  out,  conservation  of 

(1  -  M(y)K) L  dy 


(1  -  M(y)K) 


amounts  to  continuity  of  acoustic  particle  displacement. 


In  case  of  the  Ribner-Miles  problem  of  Figure  II-7,  K  =  [cos<j>/(l  4-  M  cos<j>)] 

-1  ° 

and  the  angle  of  emission  8  is  given  by  cos  (K) .  Assume  for  M(y)  a  profile 

of  type  M  sin  Cr  ^)  for  0  <  y  <  6.  Then  we  may  show  that 
o  l  o  —  — 


<t~,v  2K.,i  -i,  {M  K  + - i - 

511  Mo  *  1  ° 

o 


M  K 

+  tan"1  ( - — - 


)]}-!] 


/I  -  M‘KZ 

o 


[tan 


,  /I  -  M  K 

"  i\/- _ 

V  1  +  M  K 


x2  .  6(1  +  -I-  -  4  -f-> 


These  results  now  suffice  to  determine  a  reflection  coefficient  R  defined  as 

follows.  The  incident  pressure  wave  is  denoted  by  esp [j  {k(xcos<{>  -  (y — _ 

r  J  J  1  +  M  cos<{)  J 

mt}].  R  depends  on  <j>,  H  and  (6/A)  where  A  =  [27i/k].  The  case  (6/A)  =  0 

°  8  9 

corresponds  to  the  case  calculated  by  Ribner-Miles  ’  . 

Some  calculations  of  |r|  and  arg(R)  are  tabulated  in  Appendix  II-l  for 
Mq  =»  0.5,  0.7  and  0.9,  for  0  <_  £  0.18.  (The  value  of  0.18  corresponds  to 

(k6)  -  1  and  it  was  felt  that  the  0(5/A)  calculation  would  not  be  meaningful 
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for  (k6)  >>  1).  The  tabulation  in  Appendix  11-1  shews  that  in  roost  cases  |r| 
decreases  with  increasing  (6/A)  in  accord  with  the  physical  expectation  that 
the  reflection  coefficient  should  be  weaker  for  a  profile  with  continuous  shear 
as  compared  to  the  case  of  a  step  function  change  in  velocity.  In  five  of  the 
thirty  cases  tabulated  (marked  by  an  arrow  in  Appendix  II-l),  J R [  records  an 
increase  with  increasing  (6/A).  Evidently  these  correspond  to  cases  of  some 
unusual  constructive  interference  of  waves  reflected  continuously  from  the  shear 
layer.  This  conjecture  is  supported  by  the  observation  that  arg(R)  in  these 
cases  exhibits  an  unusual  trend  (with  increasing  (6/ A))  being  either  the  reverse 
of  the  usual  trend  of  variation  of  arg(R)  with  (6/A)  or  exhibiting  unusually 
large  changes  of  arg(R). 

Adaptation  of  a  procedure  similar  to  the  above  to  the  cylindrical  shear 

flow  geometry  of  Figure  23  yields  the  following  matching  conditions  between 

dp  ,  rv  .  do 

p,  at  r  =  (a  -  o)  and  p,  •  at  r  =  a. 
r  dr  r  dr 

,c!n.  _ ,do. 

dr  (1  -  M  K) '  'dr' 

r  -  a  r  =  (a  -  5) 


+  k2  c  (r  =  (a  -  S)  or  r  -  a)  .l3 


(55) 


and 


p  (r  =  a)  =  p  (r  =  a  -  5) 


(— ) 
'dr ' 


r  =  (a  -  6) 


(1  - 


MoK) 


(56) 


where 


I 


3 


I 


a 


6 


K2  dr _ 

(1  -  "k  (r )  kT* 


6 


(57) 
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and 


(1  -  M(r)K)2  cr 


(58) 


The  integrals  13*  I4  are  again  evaluated  assuming  a  half  sine  profile  for  the 
transition  from  M0  to  the  ambient  as  sketched  in  Figure  23.  The  calculations 
of  Section  1.1  were  revised  using  (55) ,  (56)  in  place  of  the  p,  n  matching 
employed  therein.  In  Figures  24  and  25,  calculations  for  the  total  power 
normalized  by  that  of  the  freely  moving  source  are  shown  for  (l^a)  in  the 
range  0.01  to  1,  and  (6/a)  in  the  range  0  to  0.2.  (The  case  6/a  ■  0  corres¬ 
ponds  to  the  previous  calculation  of  Section  1.1).  As  can  be  seen  there  is  a 
small  (1-2  dB)  effect  on  the  total  power  results. 

1.4  Jet  Density  Exponent  Issue  for  Noise  of  Heated  Subsonic  Jets 

The  subject  matter  of  the  present  section  is  the  question  of  how  the 
sound  power  of  a  jet  of  constant  exit  velocity  would  vary  as  the  jet  exit 
density  is  varied.  Changes  in  jet  exit  density  would  inevitably  be  accomplished 
in  a  real  experiment  by  changes  in  speed  of  sound  (temperature)  of  the  jet  so 
that  both  effects  must  be  considered  simultaneously.  The  point  of  view  advanced 
at  the  end  of  the  section  is  that  experimentally  observed  results  in  this  area 
seem  to  admit  of  an  explanation  based  how  the  radiative  efficiency  of  moving 
acoustic  sources  is  affected  by  the  shrouding  effect  of  a  jet  flow  whose 
velocity,  temperature  and  density  differ  from  those  of  the  ambient.  This 
change  of  efficiency  is  calculated  with  the  aid  of  a  simple  model  problem  as 
follows.  We  determine  the  acoustic  power  output  of  a  convected  monopole 
source,  simple  harmonic  in  its  own  frame  of  reference,  moving  along  the  ax4s 
of  a  slug  flow,  round  jet  whose  velocity  is  taken  to  be  the  same  as  that  of  the 
source.  The  jet  is  doubly  infinite  and  the  source  is  assumed  to  be  of  infinite 
lifetime.  The  jet  density  and  temperature  are  allowed  to  be  different  from 
that  of  the  ambient  though  the  specific  heat  ratio  of  the  jet  fluid  is  assumed 
to  be  the  same  as  that  of  the  ambient.  The  requirement  of  equality  of  the 
static  pressure  inside  and  outside  the  jet  then  calls  for  a  certain  restraint 
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on  how  the  jet  density  and  temperature  must  vary.  For  a  specific  value  of  jet 
exit  velocity,  the  variation  of  acoustic  power  with  the  ratio  of  jet  density 
to  the  ambient  density  along  with  a  simple  assumption  on  how  the  source  strength 
would  vary  with  jet  density  is  employed  to  theoretically  deduce  the  "jet 
density  exponent"  for  jets  which  are  subsonic  with  respect  to  the  ambient  speed 
of  sound.  The  jet  density  exponent  is  found  to  depend  both  on  the  jet  Mach 
number  and  even  more  strongly  on  a  source  frequency  parameter.  The  theoretical 
results  are  compared  to  some  experimental  studies  of  this  problem.  Encouraging 
agreement  is  obtained  both  for  the  detailed  observed  effects  on  the  power 
spectrum  and  the  exponent  for  the  overall  power. 


The  model  problem  is  sketched  in  Figure  26.  We  wish  to  determine  the 

sound  field  due  to  a  fluctuating  monopole  point  source  translating  at  a  uniform 

subsonic  velocity  U  (where  U  <  cq,  cq  being  the  speed  of  sound  of  the  ambient). 

The  source  translates  along  the  axis  of  a  round  jet  whose  velocity  profile  we 

assume  to  be  a  slug  flow  velocity  profile.  Also  the  jet  velocity  is  taken 

equal  to  hat  of  the  source.  The  source  is  assumed  to  have  a  time  dependence 

in  its  owr  frame  of  reference  of  q^  cos(uq  t) .  The  mean  jet  density  and  speed 

of  sound  are  taken  as  c  ,  c,  while  those  of  the  ambient  are  denoted  by  p  c  . 

A  1  2  o  o 

Now  the  static  pressure  insdie  the  jet  is  given  by  p  =  c^/y^  and  similarly  of 

the  ambient  by  pQ  c^/y^.  Since  the  static  pressures  inside  and  outside  the  jet 

nust  be  equal,  if  we  assume  that  -  YQ  (a  reasonable  assumption  for  the  heated 

jet  situation  but  less  valid  if  foreign  gases  such  as  Freon,  etc.,  are  used  for 

2  2 

the  jet  fluid),  then  we  must  have  p^  c^  =  p^  to  balance  the  static  pressures. 

This  implies  a  coupling  between  the  density  ratio  and  speed  of  sound  ratio 

_2 

(i.e.  =  (c^/c^)  )  which  is  always  employed  in  the  current  study. 

Analytically  we  wish  to  determine  an  acoustic  velocity  potential  <|>  which 
satisfies  in  region  I  of  Figure  26  (outside  the  jet) 


V?0 


0 


(59) 
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and  in  region  II  (inside  the  jet)  of  Figure  26 


and  <|>r  =  (nt  +  U  nx)  in  region  II.  (64) 

Outside  the  jet,  i.e.  in  region  II,  the  velocity  potential  <|>  is  also 
subject  to  a  radiation  condition  which  states  that  only  outgoing  waves  be 
emitted  by  the  moving  source. 

Before  proceeding  further,  it  is  necessary  to  discuss  the  implications 
of  employing  a  slug  flow  or  top  hat  velocity  profile  which  is  known  to  be 
unstable  when  excited  by  certain  wave  number-  frequency  combinations  of 
longitudinally  traveling  waves.  The  general  procedure  for  examining  the 
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instability  issue  (e.g.  Batchelor  and  Gill  consider  the  incompressible 

case)  is  t  consider  the  unforced  jet  eigenvalue  problem.  In  other  words,  a 

solution  of  type  R(r)  exp(j(kx  -  u-t))  is  assumed  for  $  whence,  in  the  absence 

of  any  ouri o  term  driving  the  systems,  specifying  real  w  determines  k  as  a 

lunctiu.i  of  ..  (spatial  stability  analysis)  or  alternatively  specifying  real 

k  determines  as  a  function  of  k  (temporal  stability  analysis).  The  imaginary 

parts  oj  k  or  ...  respectively  determine  the  regimes  of  instability.  In  the 

present  problem  a  source  of  type  cos(.  t)o(x  -  lit)  imposes  a  rather  specific 

t  of  traveling  wave  disturbance  on  the  jet  column.  By  Fourier  decomposi- 

f  ion  of  the  source  excitation,  the  traveling  waves  turn  out  to  be  of  type 
— 

esp  j  [  (u  +  -  at]  with  real  w.  Thus  formal  consideration  of  the  infinite 

lifetime  source  problem  yeilds  simple  traveling  wave  excitation  of  the  jet 
column  which  produces  either  a  propagating  or  a  decaying  sound  field  outside 
the  jet  depending  on  whether  the  wave  speed  parallel  to  the  jet  axis  which  is 

II  1  1  +  ^—1  ^1  exceeds  c  or  not.  This  then  leads  to  the  obvious  result  that 

...  1  o 

oiteciive  acoustic  power  is  produced  in  the  far  field  over  a  frequency  range 

-j  w  2 

: — — — ■ —  ■  a  <  - — —  where  M  where  M  =  U/c  (Morse  and  Ingard  as 
1  +  M  -  o  —  1  -  M  o  o  o  & 

o  o 

predicted  by  the  Doppler  shift  formula). 

It  seems  likely,  however,  that  if  one  sought  the  solution  to  the  problem 

as  the  limit  of  an  initial  value  problem,  i.e.  assumed  a  source  strength  of 

type  ‘(y)i(z)  cos(nfl  t)4(z)  cos(wo  t)‘(x  -  Ut)  H(t  -  tQ)  where  H(t)  is  the 

unit  step  function  and  then  studies  the  limit  of  the  solution  as  t  +  one 
r  o 

would  find  that  the  "starting  up"  process  of  the  source  triggers  the  Kelvln- 
Helmlioltz  instability  at  the  jet  still  air  interface.  (If  the  source  is  assumed 
switched  on  only  at  the  time  t  =  tQ,  it  would  be  required  that  <p  =  0  (or 
constant)  for  t  <_  tQ.)  To  ensure  this  requires  that  in  the  Fourier  integral 
representation  of  4>  the  path  of  integration  in  the  w-plane  be  specified  in  a 
certain  manner.  Presumably,  then,  in  deforming  that  path  of  integration  on 
to  the  real  axis  of  u(-“>  <  u  <  <*>) ,  unstable  pole  contributions  corresponding 
to  the  excitation  of  the  instability  modes  would  be  picked  up  in  addition  to 
tlit-  contribution  from  integration  over  the  real  w-axis  which  alone  is  discussed 
herein. 
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This  aspect  of  the  problem,  while  undoubtedly  a  difficulty  with  the  slug 

flow  profile  or  indeed  any  mean  flow  velocity  profile  that  is  inflectional, 

is  ignored  in  the  present  study  based  on  the  following  physical  argument. 

Practical  jets  (at  sufficiently  high  Reynolds  numbers)  do  represent  a  stable 

flow  situation  though  characterized  by  high  turbulence  levels  (the  rms 

turbulence  level  can  often  be  as  high  as  15Z  of  the  jet  velocity).  The  real 

jet  flow  then  represents  a  flow  with  a  distribution  of  mean  velocity  and 

turbulence  levels  which  is  manifestly  stable  to  source  excitations  of  the  type 

that  lead  to  jet  noise.  The  high  turbulence  level  in  the  jet  itself  could  be  a 

stabilizing  agent  by  a  mechanism  of  eddy  viscosity  as  has  been  noted  by 

several  previous  studies  of  turbulent  shear  flows.  For  example,  Landahl^ 

surmised  that  even  for  turbulent  boundary  layers  (generally  characterised  by 

lower  turbulence  levels  than  a  jet  flow)  the  eddy  viscosity  seesmd  about 

12 

80  times  as  effective  as  the  molecular  viscosity.  Similarly  Bishop  ,  et  al., 
specifically  suggested,  with  respect  to  high  speed  jets,  the  substantial 
lowering  of  the  effective  jet  Reynolds  number  due  to  the  eddy  viscosity.  The 
justification  for  the  use  of  a  slug  flow  velocity  profile  then  rests  on  the 
fact  in  attempting  to  infer  the  effect  of  the  more  complicated  mean  velocity 
profile  of  the  pure  jet  on  the  radiative  efficiency  of  a  moving  source,  one 
may  exploit  the  relatively  low  frequency  nature  of  jet  noise  sources  to  argue 
that  the  precise  nature  of  the  true  velocity  profile  need  not  be  retained. 

Since  we  know  that  the  true  jet  flow  is  stable  and  it  is  understood  that  the 
slug  flow  profile  is  only  employed  as  an  analytical  artifice  to  conveniently 
assess  the  shrouding  effect  of  the  flow  on  the  radiative  efficiency  of  the  source, 
we  may  then  reject  the  unstable  excitation  of  the  slug  flow  profile  owing  to 
the  start  up  process  of  the  source  (which  appears  in  the  rigorous  analytical 
solution  of  (59)  and  (60)  when  a  proper  initial  value  problem  is  posed)  aa  not 
germane  to  the  real  physical  problem.  In  what  follows  then  we  will  deal  only 
with  the  stable  and  bounded  solution  to  (59)  and  (60)  subject  to  the  matching 
conditions  of  (61)  -  (64)  assuming  a  source  of  Infinite  lifetime.  The  problem 
posed  by  (59)  and  (60)  is  a  transient  one  and  the  required  bounded  solution  may 
be  obtained  by  formally  applying  the  Fourier  integral  method  taking  w  real. 

A  similar  difficulty  arises,  of  course,  in  the  calculations  of  Section  1.1 
through  the  arguments  for  ignoring  the  stability  issue  were  not  spelled  out 
in  as  much  detail  as  above.  From  a  fluid  mechanics  rather  than  a  mathematical 


. 
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point  of  view  the  case  for  both  the  present  calculations  as  well  as  those  of 
the  earlier  study  (Section  1.1)  must  perhaps  be  judged  on  the  degree  to  which 
the  results  are  in  accord  with  physically  observed  features  of  jet  noise. 


Tlie  analysis  that  follows  parallels  closely  that  of  Section  1.1  and  is 
given  below  mainly  for  completeness  and  to  clearly  delineate  the  new  elements 


that  arise  due  to  p^,  c^  not  being  equal  to  p^,  cq. 


Let  p,  etc.,  denote  the  Fourier  transforms  with  respect  to  time  of 
the  corresponding  physical  quantities.  Thus 


"fe  I 


<f>  e-^*1  dt 


j  =  /=1 


(65) 


and 


■  J 


5  e-iut  du. 


Also,  we  write  cos(u>o  t)  =  -(exp(jui^  t)  +  exp(-jiao  t)J.  The  problem  for  the 


transforms  is 


V2$  +  k =  0, 


(ko  =  u/co> 


in  region  I, 


(59') 


q  <5(y)5(z' 

U  -  K? )  «xx  +  Vji  ♦  2jk1K1Sx  +  — 


(ezp  ( 


j  (oJ  -  u  )  x 


-° — )  +  exp(- 


j  (a)  +  to  )  x 


in  region  II;  (60') 


U 
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■A 

■t 


1 

I 


1 
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I 

f 

c 


where  and  =  u/c1 


P  =  jwP0$, 

P  -  "P1(tjw?  +  U$  ), 

$r  = 

$r  ”  -jwh  +  Un^, 


in  region  I, 
in  region  II. 
in  region  I, 
in  region  II. 


(61 ») 
(62' ) 
(63') 
(64’) 


^  ’  '  ,+  +  r  Sl"Uarly  £ot  “  *“■*  n  --re  *+  corresponds  to  the  solution 

wit  t„e  ter.  ,„J  (.  -  in  e,uatio„  (60- ,  end  to  the  ter.  involving 

exp{j (u  +  u  )x/U}.  8 


Consider  in  detail  the  proble.  for  Intuitively,  it  is  clear  that 

’  V  n  a“  hav'  “  x-‘‘epe"‘ie“<:e  of  the  type  exp(j(.  -  .)*/„}.  ■•Factoring" 
this  dependence  out,  one  is  left  lmfk  e  n  0  8 

,  one  is  left  with  the  following  problem  in  the  y-2  plane: 


vlV  +  k+27>+  = 


in  region  I; 


qnMy)Mz) 

Vf$  +  k  +  2o  =  -2 _  , 

A  rr  r\  II  9 


4ir  p  U 

o 


in  region  II; 


wliere 


1  —  M*  u)  /c 

k+!  =  <(r^-)=0)(>to 


to  /c 

— o  0)1 

1  +  M 

o 


o  u 

and  is  >_  0  only  if  j U  * 


:+!  ■  Ij  "ax11  +  V  ki)(ki  -  sx(1  -  H1>»! 
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H  »  Owr  • '  MWIW 


„ _ „  -  — <  — 


Also  let 


+  2  +2 
K  ~  —  K  - 


“  S[('  ki  -  4'1  +  "i’X*!  '  57(1  -  Hl>)> 


(l)  (i) 

We  restrict  our  attention  to  the  frequency  ran^e  ~z - —  >  gj  >  —  - —.  The 

L  —  M  —  —  1  +  M 

matching  conditions  for  equations  (68)  and  (66)  are  that,  at  r  =  a, 


<>+  (r  =  a+)  = 


p,  u  . 

A.  _-°-  (Ji  (r  =  a  )  (pressure  matching  (70) 


PD  “ 


<j>+(r  =  a+) 


5L-  o+(r  =  a") 
u  vr 
o 


condition) , 


(transverse  particle  (71) 
displacement  matching 
condition) . 


To  solve  equations  (65),  (66),  (67),  and  (68)  in  the  range  r — , ■  <  ft  <  - — - 

1  T  Mq  X  —  Mq 

and  with  restriction  to  outgoing  waves  at  infinity,  in  the  range  we  assume, 
for  J+  in  regions  I  and  II, 


in  I  : 


,  i  (u  -  f.  )  x 

r  r,  a*  II  ])  (h'rjox?  - - XT"” 

I  o 


in  II:  l  =  [Att  J  (<  r)  - 

II  o 


jqo"o:)  (i;Tr) 
16r p^U 


1  exp 


j  (w  -  cj  )>: 
o 


if  u  f  [us  (1  -  n.  )f  «  (1  +  M,  )] 
o  1  o  1 
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Otherwise  in  region  II, 


q  K  (<  r)  j(“>  “  wJx 

♦+  -  Ki  V« r>  -  ~f;°p^u  '  e!!p - r2- 


Equations  (72)  and  (73)  or  (72)  and  (74)  may  now  be  readily  solved  for 
amd  A  by  using  the  matching  conditions  (70)  and  (71).  Since  we  are 
interested  in  far  field  pressures,  far  outside  the  jet,  we  only  give  the 
result  for  A*: 


f  w  C  [uo<l  -  Mx)  ,  wo(l  +  M  )]  , 


q  k+[Y  (<+a)  J1  (tc+a) 

w _ O _  _ -4- _ 

16t,p  k+  H<1}  (k+a)J.  (<+a) 

O  10  O  1 

o 


Y1(<+a)Jo(<+a)] 
k+  (►c+a)H^1)  (k+a)]U 

p  OJ  O  1 

o 

(75a) 


(b)  otherwise. 


-q  c+[K  (<  a)I,(<+a)  +  I  (<+a) K, (<+a) ] 

o  —  O  —  j.  — _ o  —  1  —  _ 


8ti2p  U[  — 
o  w 


I  (ic+a)H*1)(k+a)  +  —  k+  I  (tc+a)H(1)  (k^a)  ] 

1—0  0)  0—1 

o 

(75b) 


Equation  (75)  essentially  completes  the  formal  solution  to  the  problem.  The 
far  field  pressure  and  the  radial  acoustic  velocity  may  be  computed  by  using 
£  =  jwp  $  and  In  this  problem,  every  point  on  a  cylindrical  surface 

concentric  with  the  jet  experiences  the  same  pressure  time  history.  Morse  and 


lngard  have  discussed  thoroughly  the  problem  of  determining  the  power 
spectrum  and  total  power  radiated  by  the  source  and  their  concluding  result  is 
that  the  power  spectral  density  extends  over  a  frequency  range  [w0/(l  +  Mq)  ] 

<  w  <  [u  /(l  -  M0) J  and  is  given  by 

(U..P0Uw)|A*|J  »  1M  (76) 

The  total  power  is  given  by 

u  /(I  -  M  ) 

to  o 

P  =  I(t*i)du  (77) 

uo/(l  +  Mq) 

1 . 5  Computed  Results  and  Inferences 

The  power  calculated  in  (77)  is  nondimensionalized  first  by  the  power  of 

o  2  2  2 

a  freely  moving  source  which  is  q*  iWStt  pQ  cQ  (1  -  Mq)  .  Fixing  Mq  and 
(m  a/rU),  the  nondimens ional  power  (say  P')  against  log(pj/p0)  by  a  least 
squares  fit  an  exponent  n'  is  determined  for  each  Mq  and  (u>0a/irU).  (ti>0a/irU) 
may  be  termed  a  source  Strouhal  number  (e.g.  as  in  Lush  (1971)).  Now  the 

source  strength  qQ  itself  will  vary  linearly  with  jet  density  whether  one  uses 

13  14 

the  quadrpole  model  of  Lighthill  or  the  fluid  dilatation  D;odel  of  Ribner 

Since  p'  itself  varies  as  p"  and  P'  is  the  power  normalized  by  q^,  with 

qQ  -  p^«  one  would  then  expect  the  actual  power  to  vary  with  density  p^  as  with 

an  exponent  (n'  +  2)  =  n. 

This  theoretically  deduced  exponent  n  is  plotted  in  Figures  27,  28  and  29 
for  Mq  =  0.5,  0.7  and  0.9  and  a  range  of  source  Strouhal  numbers  in  range  0.1 
to  1.0.  For  source  Strouhal  numbers  in  excess  of  1,  owing  to  the  high 
frequencies  involved  the  slug  flow  model  would  be  less  adequate. 

If  we  take  the  limit  of  75(a)  or  (b)  at  very  low  frequencies  one  can 
readily  show  that  the  index  n  would  be  expected  on  the  basis  of  this  model  to 
tend  to  zero.  In  general,  then,  for  subsonic  Mach  numbers,  the  present  model 
predicts  that  n  +  0  as  the  frequency  parameter  approaches  zero,  and  it's  then 
negative  for  a  range  of  frequencies  and  finally  starts  increasing  monotonically 
with  frequency.  It  is  not  possible  to  analytically  extract  a  high  frequency 
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limit  from  75(a)  and  (b)  but  if  one  used  the  argument  that  at  high  frequencies 

the  source  output  is  determined  only  by  its  own  immediate  ambient,  the  exponent 

2  -1/2 
should  depend  on  how  (qQ/PiC^)  varies  with  pi  and  since  q0  ~  pj_,  and  c^  ~  , 

for  a  monopole  source  model,  the  exponent  would  tend  to  1.5  at  high  frequencies. 

The  intrinsic  source  distributions  generating  jet  noise  do  exhibit  Strouhal 

scaling  with,  respect  to  velocity  (this  is  confirmed  either  by  in-jet  measurements 

or  by  looking  at  the  90°  point  far  field  data  where  convective-refractive 

effects  are  absent)  so  that  high  speeds  do  go  with  high  true  source  frequencies 

and  vice  versa.  In  Figure  32,  the  exponents  for  Mq  =  0.5,  0.7,  and  0.9  are 

shown  as  a  function  of  a  true  frequency  parameter  (w0/acQ).  There  is  a  general 

trend  to  exponents  of  value  zero  as  (u^a /cQ)  -*■  0,  followed  by  a  region  of 

negative  exponents  and  a  tendency  for  n  to  attain  values  of  1.5  for  high  values 

of  (u0a/c0)  almost  independent  of  jet  Mach  number.  In  view  of  the  Strouhal 

scaling  with  respect  to  velocity  exhibited  by  the  source  distributions,  one 

would  expect  (in  terms  of  Figure  32)  the  higher  jet  velocities  to  go  with 

higher  values  of  (uQa/cQ).  (ui0a/c0  is  (hMq)  times  the  source  Strouhal  number). 

In  other  words  Figure  32  indicates  that  even  in  terms  of  jet  velocities  one 

would  expect  a  changing  exponent  (say  for  the  total  power)  starting  off  at 

zero  at  the  lowest  velocities,  then  being  negative  and  then  finally  increasing 

monotonically  with  velocity.  The  present  calculations  are  of  course  limited 

to  subsonic  jet  velocities. 

In  Figure  31,  we  show  first  the  empirical  result  obtained  for  the  exponent 
n  for  the  total  power  obtained  by  Hoch^,  et  al.  As  indicated  earlier, 

Hoch,  et  al.,  find  that  n  is  a  junction  of  jet  Mach  number.  In  order  to  compare 
the  present  analysis  with  the  data  of  Hoch,  et  al.,  it  is  necessary  to  estimate 
source  Strouhal  numbers  representative  of  the  total  power.  Based  on  jot  noise 
at  low  Mach  numbers  (wherein  refractive,  shrouding  and  Doppler  shift  effects 
should  be  negligible),  it  was  felt  that  a  source  Strouhal  number  somewhere 
between  0.3  and  0.6  would  represent  a  "typical"  source  Strouhal  number  for 
assessing  an  exponent  for  the  total  power.  Shown  in  Figure  31  for  1^,  in  range 
0.5  to  0.95  are  results  fo  n  of  the  present  study  for  source  Strouhal  numbers 
of  both  0.3  and  0.6.  Except  at  the  lowest  Mach  number  of  0.5,  the  predicted 
values  of  the  exponent  for  source  Strouhal  numbers  of  0.3  and  0.6  bracket  the 
experimental  values  of  Hoch,  et  al.  quite  well. 


It  appears  from  the  present  study  that  the  differences  in  noise  between 
a  heated  and  a  cold  jet  do  adnit  largely  of  an  acoustic  explanation  being 
attributable  to  the  effect  on  the  radiative  efficiency  of  a  moving  source  due 
to  the  mismatch  of  velocity,  density  and  temperature  within  and  outside  the 
jet.  The  l.ict  that  Figure  30  taken  from  Koch,  et  al.,  is  in  accord  with 
Figures  27-29  appears  to  be  the  most  impressive  evidence  of  this  for  it  is 
diffucult  to  conceive  of  other  explanations  based  on  entropy  fluctuations, 
jet  mixing,  etc.,  that  would  explain  the  tendency  of  heating  to  raise  the  low 
frequency  end  of  the  power  spectrum  while  depressing  the  high  frequency  end. 

As  noted  earlier,  extensions  to  higher  order  multipoles  and  to  sources  of 
finite  lifetime  are  undoubtedly  needed  but  the  least  the  present  study  may  be 
said  to  achieve  is  to  indicate  the  profitability  of  pursuing  such  analyses. 

1 . 6  Directivity  of  Subsonic  Jet.  Noise 

In  the  last  Subsection,  we  will  calculate  the  directivity  of  a  fluctuating, 
simple  harmonic  point  source  embedded  in  a  slug  flow,  cold  round  jet  with  the 
source  converting  along  the  jet  axis  at  a  velocity  different  from  the  jet 
velocity.  It  is  required  that  the  convection  velocities  of  the  source  or  eddy 
be  subsonic.  In  the  calculations  of  directivity  to  be  shown  later  the  eddy 
convection  velocity  is  always  assumed  to  be  0.65  times  the  jet  velocity. 

Initially,  we  start  with  the  case  of  a  moving  source  of  iime  dependence 
e^o11  in  a  stationary  medium.  The  objective  is  to  calculate  its  directivity 
as  opposed  to  its  power  spectrum.  Since  interest  is  in  directivity  rather 
.ban  a  power  spectrum  it  is  obvious  that  rather  than  a  transform  in  time  one 
in  space  needs  to  be  employed. 

The  velocity  potential  $(x,r,t)  satisfies  (cf.  Figure  33): 


1  9_ 
r  9r 


/  3*. 

lr  3?’ 


q  ju  t 

-9.  0  °  6  (x  -  Met)  6  (y)  <5  (z) 

P0  (78) 
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As  noted  earlier,  at  the  lower  Mach  numbers  (and  assc  ated  low  true 
frequencies),  the  present  analysis  predicts  the  exponent  to  tend  to  zero.  This 
result  is  apparently  at  variance  with  the  SNECMA-NGTE  study  of  Hoch,  et  al. 
There  are  two  points  to  note  in  this  regard: 


1)  The  present  study  needs  extension  to  higher  order  multipoles  as  well 
as  to  sources  convecting  at  velocities  different  from  (less  than)  the  jet  speed. 
So  also  the  extension  to  finite  source  lifetimes  is  also  needed.  It  is  not 
clear  how  much  these  extensions  will  alter  the  theoretical  predictions  of  the 
exponent. 


2)  It  is  in  the  lower  velocity  r^nge  that  isolation  of  the  jet  density 
exponent  associated  with  pure  jet  noise  becomes  most  difficult  from  an  experi¬ 
mental  point  of  view.  This  is  because  cc  the  ever  present  danger  of  internal 
noise  sources  termed  "parasitic"  noise  by  Hoch,  et  al.)  such  as  valve  noise, 
combustion  noise,  etc.  One  can  easily  show  that  the  effect  of  tuch  a  combustion 
noise  source  will  be  to  lower  the  effective  index  from  its  value  for  pure  jet 
noise  in  an  experimental  situation.  In  Hoch^,  et  al.,  it  is  pointed  out 
that  the  NGTE  group  worked  in  the  low  velocity  end  while  the  SNECMA  group 
worked  at  the  high  velocity  end.  Hoch,  et  al.,  have  cited  the  good  agreement 
between  the  results  of  the  two  groups  in  the  region  of  overlap  of  velocities  as 
one  indication  of  the  internal  cleanliness  of  their  facility.  However  Hoch, 
et  al.,  do  point  out  that  the  region  of  overlap  extends  from  Mq  =  0.6  on  up  so 
that  "pure  jet  noise  is  being  measured,  at  least,  above  jet  velocities  of 
200  m/s"  (Mq  =  0.6).  We  may  note  that  an  earlier  experimental  study  by 
Rollin^  did  conclude  that  the  density  exponent  was  zero. 


Finally,  in  Figure  30,  we  show  the  detailed  effects  on  the  power  spectrum 
due  to  heating  observed  by  Hoch,  et  al.,  at  two  velocities  corresponding  to 
Mo  =  0.6  and  Mq  =  1.2.  It  is  observed  that  they  find  that  heating  increases 
the  low  frequency  portions  of  the  spectrum  while  depressing  the  higher 
frequencies.  This  is  fully  in  accord  with  Figures  27-29  wherein  (as  indicated) 
indices  greater  than  zero  correspond  to  portions  of  the  frequncy  spectrum 
lowered  upon  heating  and  indices  less  than  zero  to  portions  raised  by  heating. 
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so  that 
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«  .-3«2 


Then 


I  L.  U  il) 

r  Sr  1  Sr' 
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1  S2$ 
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j»  t  jeVt 
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°-e  °  c  6  (y )  6  ( 2 ) 


(79) 


(80) 


(81) 


where  Me  =  V. 

Clearly  depends  on  time  as  exp{j(u>  +  aV)t}.  Thus  the  radial  dependence  of 
$  is  to  be  determined  from: 


1  d_ 
r  dr 


(r 


+  (1  -  M*) ( 


1  -  M 


c)  (a  + 


k 

o 

1  +  K 


)$ 


%_ 

2ve o 


6  (y)  5  (r.) 
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(83) 


Let  (Figure  33)  r  »  R’  sin  0*.  and  (x  -  Vt)  •  R'  cos  6*.  For  large  R*  and 
6’  i  0,  t,  r  will  also  be  large  and  $  will  be 


The  exponential  terns  depending  on  a  inside  the  integral  aay  be  written  as: 
exp[jR'  h(a)J  where 

h(a)  *  -  [a+  sin0*  +  acosS'J  ...  (89) 
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am. 


and  the  integral  for  $  nay  be  written  as: 


$  - 


(i  n/4) 
c  3  q  e 

J  'O 


iid  t 

J  o 


br:p 


o 


w 

/2A7  J 


jK'h(a) 


da 


Let  denote  the  value  of  u  for  which  h(a)  is  stationary,  i. 
Then  by  the  ..ethod  of  stationary  phase  (cf.  e.g.  A.  Erd^lyi, 
Expansions,"  p.  51),  for  large  R' , 


$ 


C  cr.ol 


*>  * 


h  (=0) 


/a’  iaQ) 


2tt 

h"  (a0)J 


( tiiis  form  assumes  that  a^  is  a  minimum  of  h(a),  i.e.  that  h,f 
C  is  just 


j  w/4 


j  q  e 


j  03  t 
J  O 


8-<p 


/2/«r  ... 


h(a)  *  -  [a  sin©'  +  acos©’ ]  ... 


is  stationary  when 


cin  (6 '  ) 


+  cosG' 


0  or  when 


(1  -  K* ) {  ( 


o 


1  -  M 


-a)  -  a  - 


1  +  M 


}  sin9'  +  2a  cos0 


(90) 

.  h' (a)  -  0. 
Asymptotic 

(91) 

(oq)  >  0)  where 

(92) 

(93) 

*  =  0  (94) 
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•  v  nu,  ■  *«yr*  .  •  **.■>* 


or  {2M  k 

o 


2a  (1  - 


M2)  } 


sin6*  +  2a  cosO1 


0.  Let  6  =  k  +  aM.  Then 
o 


(MB  -  c.)  sinG'  =  -cosO'  f$v  -  a2 


(95) 


After  some  painful  algebra  we  find  that 


cosO ' ( /I  -  M‘  sin‘6’  +  M  sin20']k 
_ _ _  o _ 

/I  -  K*  sin' 6'  { /I"  -  M'  sin*  6  '  -  M  cosG} 


(96) 


6',  R’  denote  (in  Figure  33)  the  angle  from  line  of  surce  motion  and 
distance  from  source  location  of  the  far  field  point  at  current  time.  Denote 
by  0,  R  these  same  quantities  measured  from  the  source  location  at  the  time 
that  the  radiation  reaching  the  observer  at  current  time  was  emitted.  6,  R 
are  shown  in  Figure  33.  Straightforward  kinematics  leads  to: 


tan (0' ) 


sin6 

cos9  -  M 
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or 


s  in  ( 6 '  ) 


_ sir.  9 

►1  +  -  2M  cosQ 


(98) 


and  cos (6 ' ) 


(cos6  -  M) _ 

•  T  +  M  -  2M  cos1- 


In  terms  of  9, 

k  cor.  0 

o _ 

o  (1  -  M  cos3j 

A  tedious  calculation  of  h''(a  >  from  (93)  shows  that. 


+  h'  ’  (ao) 


(1  -  M  cosa)  3 _ _ 

k  /I  4  M;  2M  CCS *3  sin7  0 
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(99) 


(100) 
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confirming  that  h"  (aQ)  >  0  (for  0  <_  0  <_  it).  Also  corresponding  to  (100) 


a+(a  ) 


k  sind 
o 


o  (1  -  M  cos0) 


(102) 


a  M  +  k  =  — - — - — 

o  o  (1  -  M  cos6) 


(103) 


h(a  ) 


-k  cos (6  -  0') 
o 


o  (1  -  M  cos0) 


(104) 


Since  p  »  -p  from  (87)  we  may  deduce  that, 

O  dt 


p  = 


qo  e 


ju  t 
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87TP 


f  (wo  +  aV)H<J2)  (a+r)e"ja(x  “ 


Vt) 


da  (105) 


The  point  of  stationary  phase  for  p  is  clearly  the  same  as  for  $  and  hence  in 
the  far  field  p  may  be  ap- -oximated  by 
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3  %  lo 


4  "T  R  ( 1 


joi  t  -ik  R 
o  _  *  o 


-  K  cost?) 


(106) 


This  result  is  in  complete  accord  with  the  far  field  part  of  equation  (11.2.15) 
of  Morse  and  Ingard's  "Theoretical  Acoustics"  (p.  724). 


1.6.1  Extension  to  a  Slug  Flow,  Cold  Jet  (Figure  34) 

In  this  case,  in  region  1,  the  transform  '‘(r.t)  satisfies: 

13  crg.-.’o-.A. 
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_£L_  e^o1  eJ“Vct  6(y)6(lt) 
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(107)  can  be  rewritten  as: 


13  ,  3$.  /r, 

?  5T  (r  3?>  +  (!ko  -  “(M 


V1’ 


-  et2  }$ 


2irp 


6  (y )  5  (z) 


(109) 


for  0  _<  r  <  a. 

For  r  >  a,  (109)  can  be  rewritten  as: 


p  jjp  <r  |~)  +  [  (kQ  +  aMc) 2  -  a2]v>  =  0  for  r  >  a 


(110) 


As  before  interest  is  only  in  range  of  a  such  that, 


k 


<  a  < 


1  +  M  -  -  1  -  M 

c  c 


-k 


The  solution  of  (109)  for  I  a  <  k/U  -  M  +  M)  will  be  assumed  of 

form  c 


!^.+  r)  +  A  J  (a+r)  }  vh^re  a+  =  /{(k  -  a  (M  -  M  )  ] 7  -  a'  ; 

Cr  D  O  —  I  O  —  —  O  C 
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a  (M  -  M 
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(113a) 


3  5l  (a  a)qo 


a  (*1  -  M  ) 

Brr  p  (1 - - - £-) 


(113b) 


The  solution  to  (113a,  b)  for  is: 


A  (a)  =  - r-. - - - 

II  7  qM  . _ . 

4*  p  {  (1  +  (a  a)  J,  (o  a) a  a/{l 

O  iC  _  O  X 

o 


a  (M  -  M  ) 
c 


+  +  (2)  -*•  a  (M  -  M  ) 

a  a  J  (a  a)H.v  ;  (a  a)  (1  -  - i r - — 

_ o  —  1 _ _ 


(1  + 

k 

o 


(114) 


Now  consider  Wd  -  Hq  +  M  <  a  <  kQ/ (1  -  Mc).  The  matching  condition* 
(111),  (112)  yield: 


i 
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Equations  (114)  and  (116)  give  the  complete  solution  to  A  (a)  for  —  Ij°  <  Q 
(T-  MCK  The^following  limits  are  to  be  noted.  As  a  tends  to  -^/(iVm^)  or 

k°/(l  "  a  tends  to  2ero  and  iL  turns  out  that  A^c)  also  tends  to  zero 

(as  1/log  (c)  as  e  tends  to  zero).  Also  as  u  tends  to  kQ/(l  +  M  -  M*.)  a+  and 
±  tend  to  zero  and  A^  tend  to: 


-qnd  +  m) 


4tt“  pQ  (a  ’  a)  2 }  (a+a) 


(117) 


where  a  a  =  k^/M  +  2M  (1  +  M  -  Mfi).  For  r  >  a  then,  one  has: 


)u  t  ,  . 

i  °  A  v  “3a(x 

J  Ajj (a)H  (a  r)e 


V> 


(118) 


P  "  'po  3l  (ror  r  >  a)  = 


Du  t 


•J  e 


(a)  H  ^ 


(Q+r)'“0  +av  ,e'ja(X- 


vct, 


(119) 

As  before  to  get  the  far  field,  directivity  is  obtained  by  an  asy^totlc 
expansion  of  (119)^by  the  method  of  stationary  phase  for  large  values  of 

[(x  Vct)  +  r2]  .  The  point  of  stationary  phase  for  the  Integral  in  (119) 

is  clearly  the  same  as  for  the  Integral  of  (105).  Thus  p  is  given  in  the  far 
field  by: 
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(120) 


2  e 


i  ( u.1  t.  —  k  R )  „ 
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where  .t  -  k  cosG/(l 
o  o 

1  he  limit  of  A1];(uo)  as 


M  costi ) . 

(k  a)  tends  to  zero  (low  frequency 
o 


limit) 


is : 


3  % 


8tt 


-  M 


(121) 


so  that  the  low  frequency  expression  for 
point  source  in  a  slug  flow,  cold,  round 


the  far  field  of  a  moving  harmonic 
subsonic  jet  is: 


P 


j  -  V° 


-i  q 

J  Mo  o 


4 it  (1  -  Mc  cos 


n)2  R(1  -  M  COS0) 


(122) 


Contras 


t  this  with  (106)  wherein  p  is  given  by: 


•i  q  ip  e 
J  o 


j  (uj  t  -  k  R) 
J  o  o 


4trR  (1  -  Mc  cos0 ) 


b(ilI  tile  directivity  of  a  freely  moving  point  harmonic  source  (i.e.  unencumbered 
by  the  presence  of  a  shrouding  jet  flow). 

The  general  case  is  given  by  (120)  with  An<a)  given  by  <1U>,  C1U). 
a  given  by  (100)  (substituting  Mc  for  M) . 
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1.6.2  Results 


In  Figures  35-37,  we  show  first  results  of  the  present  calculations 

4 

compared  with  data  from  Lush  .  A  range  of  source  Strouhal  numbers  from 
0.03  to  1.0  and  a  range  of  jet  velocities  from  124  m/s  to  300  m/s  (corresponding 
to  jet  Mach  numbers  from  0.366  to  0.878)  was  covered  by  Lush.  In  the  theoretical 
predictions  it  is  assumed  that  owing  to  the  finite  eddy  life  time  as  well  a3  the 
measurements  being  several  diameters  away  from  the  jet  nozzle,  the  angle  from 
the  jet  axis  will  correspond  to  the  angle  measured  from  the  position  of  emission 
which  appears  in  the  theoretical  results. 

The  predictions  and  data  of  Figure  35  agree  fairly  well.  In  Figures 
36  and  37  however  we  notice  that  the  plug  flow  model  consistently  overestimates 
the  refractive  effect  especially  at  the  higher  source  frequencies.  For  example 
the  angle  at  which  the  peak  radiation  occurs  is  consistently  overestimated  by 
anywhere  from  0  to  20°  or  so.  If  we  make  an  allowance  for  this  by  one  empirical 
adjustment  whereby  we  force  the  peaks  of  the  predicted  and  measured  patterns  to 
coincide,  the  present  calculations  turn  out  to  be  in  rather  good  agreement 
with  Lush's  data.  This  is  shown  in  Figures  38-40  wherein  as  indicated  the 
theoretical  predictions  have  been  shifted  or  translated  towards  the  jet  axis  by 
varying  amounts  (the  amounts  are  shown)  so  that  the  peaks  of  the  measured  and 
predicted  patterns  coincide.  Now  except  for  one  case  of  Mj  =  0.878  and  source 
Strouhal  number  =  1.0,  there  is  extremely  good  agreement  between  predicted  and 
measured  directivity  patterns.  In  the  case  of  Mj  =  0.878  and  source  Strouhal 
number  =  1.0  the  radiation  aft  of  the  peak  angle  is  still  nicely  predicted  but 
the  refractive  dip  is  again  overestima Led . 

Tims,  allowing  for  the  fact  that  ai  overestimation  of  the  peak  angle  of 
radiation  occurs  which  can  be  corrected  for  empirically  by  adjusting  (trans¬ 
lation  towards  jet  axis)  the  theoretically  predicted  directivity  so  that  its 
peak  coincides  with  the  peak  of  the  experimental  pattern,  the  calculations  of 
this  last  section  succeed  in  explaining  the  directivity  of  cold,  subsonic 
jet  noise  over  five  octaves  of  source  Strouhal  numbers  (0.03  to  1.0)  and 
velocities  ranging  from  124  m/s  to  300  m/s. 
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2.0  AERODYNAMIC  NOISE  EMISSION  FROM  TURBULENT  SHEAR  1  AYERS 

S.P.  Pao 

In  this  section  the  Phill ips-lJao  convecLeJ  wave  equation  is  employed  to 
study  aerodynamic  noise  emission  processes  in  subsonic  and  supersonic  shear 
layers.  Particular  •’‘tention  is  directed  toward  applying  the  theory  to  study 
the  effects  of  jet  density  on  sound  emission,  the  effects  of  refraction  and 
convection,  and  how  this  theory  can  be  incorporated  into  GE  comprehensive 
aero-acoustic  prediction  methods. 

The  Phiilips-Pao  analysis  is  based  on  the  convected  wave  equation 
first  introduced  by  Phillips  in  1960.  The  convected  wave  equation  itself  is 
derived  through  the  basic  principles  of  fluid  mechanics,  and  it  is  a  natural 
extension  of  the  Lighthill  equation  of  aerodynamic  noise.  The  lineari 
version  of  the  general  equation  has  the  form  of  a  simple  wave  equa 
I.ugrangian  coordinates.  The  right-hand  side  of  this  equation  con^j..  un 

terms:  a  turbulent  quadrupole,  shear  flow  and  turbulence  interaction,  c.  ^  py 

fluctuations  and  viscous  effects.  If  the  flow  field  is  free  of  shocks,  the 
acoustic  pressure  fluctuation  can  be  assumed  to  be  decouple:1  from  the  entropy 
fluctuations.  It  is  tacitly  assumed  in  the  analysis  that  all  terms  are  the 
right-hand  side  of  the  wave  equation  are  known  quantities  and  the  contributions 
of  individual  terms  can  be  considered  as  independent  of  each  other. 

2.1  The  Effects  of  Jet  Density  on  Sound  Emission 

Many  recent  experimental  investigations  have  reported  on  the  variation 

of  sound  power  and  frequency  characteristics  with  jet  density.  Of  particular 

interest  is  that  the  sound  power  depends  on  jet  density  following  a  power  law 

2 

dependence  different  from  the  classical  Lighthill  theory  p^.  dependency.  The 
actual  experimental  index  has  been  observed  to  change  significantly  with  the 
convection  Mach  number.  The  main  objective  here  is  to  set  forth  a  model  within 
the  framework  of  the  Phillips-Pao  theory  to  account  for  such  variation. 

In  the  Philiips-Pao  theory,  the  first  order  effects  on  convection,  refraction, 
and  their  coupling  are  included  in  the  formulation.  Hence,  the  only  empirical 
elements  in  the  theory  are  the  structure  of  turbulence  and  the  mean  flow 
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properties.  In  previous  numerical  work,  it  was  found  that  the  assumption  of 
I, j  'the  scale  of  the  turbulence)  and  a  (the  eddy  decay  parameter)  influences 
strongly  the  character  of  sound  emission.  In  the  present  GE  jet  noise  projcet, 
the  range  of  jet  velocity  is  within  200  m/s  to  1000  m/s,  or  0.6  <  M  <  3.0. 
Experimental  evidences  indicated  that  the  value  of  may  vary  from  0.4  to  2.0, 
and  the  value  of  a  is  somewhere  between  0.3  and  0.6.  However,  there  it.  little 
further  information  on  the  structure  of  turbulence.  Hence,  an  empirical  model 
is  proposed.  This  model  is  based  on  three  assumptions: 

(a)  For  Mj  less  than  one,  tue  spatial  scale  remains  constant.  The 
temporal  scale  is  ar  sumed  to  vary  as  a  function  of  the  local 
speed  of  sound.  Since  frequency  of  sound  with  a  fixed  wave 
length  increases  with  the  speed  of  sound,  the  time  scale  is 
expected  to  decrease  with  an  increase  of  the  local  speed  of 
sound.  Analytically,  this  assumption  is  represented  by 


a 


a  A 
o 


0.75 


where  A 


-2 


p./p  and  a 
J  o  o 


constant 


The  index  of  0.75  is  chosen  empirically. 

(b)  For  Mj  greater  than  one,  the  scale  ratio  is  assumed  to  be  a 
function  of  the  local  true  Mach  number.  In  the  classical 
approach,  the  scale  ratio  is  assumed  to  follow  the  convection 
Mach  number. 

This  assumption  can  be  expressed  by  the  following  relation: 
A-1 

a  =  a  A 
o 


The  above  formula  follows  from  a  comparison  between  the  classical 

and  the  modified  concepts  of  the  scale  ratio: 

L,  L.  a 

— -  =  aQ  Mc  (classical);  —  ^  =  ~  (Modified) 

o  t  o  t 


(c) 


The  spatial  scale,  L  ,  is  directly  proportional  to  the  local  speed 
of  sound: 

* 
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This  assumption  permits  not  only  a  precise  modeling  o£  the  sound 
power  variation  with  jet  density,  hut  also  a  correct  shift  In 
frequency  spectrum  as  compared  with  results  given  in  Reference  15. 


The  assumption  in  (a)  and  (b)  can  be  represented  by  a  single  formula: 

A0*75  /  U  A1-7'!'1 

a  =  a  A  {a+bA  } 
o 


where 

4  4 

a  =  1-b;  b  =  {2.5  +  M* } 

The  parameters  a  and  b  are  chosen  to  represent  the  change  of  range  in 
the  convection  Mach  number.  No  other  physical  meaning  is  associated  with 
the  definition  of  these  parameters. 

Based  on  this  empirical  model  of  turbulence,  sound  emission  calculations 
are  made  according  to  the  analytical  solutions  of  the  Phillips-Pao  theory. 

T'.  nty-four  cases  of  typical  jet  noise  radiation  in  the  range  of  Mj  from  0.6 
to  2.0,  and  the  range  of  density  ratio  from  1.0  to  4.0.  It  should  be  noted 
that  the  density  ratio  is  represented  in  the  Phillips-Pao  theory  by  the  speed 
of  sound  ratio.  For  a  perfectly  expanded  jet,  these  two  quantities  are  related 
via  a  simple  formula: 

p./Po  -  A'2 

* 

In  these  calculations,  a  =0.55  and  L,  =  0.72  which  are  consistent  with 
the  current  accepted  values  for  these  parameters.  Only  one  typical  slice  of 
the  jet  near  the  end  of  the  potential  core  is  used  in  the  numerical  work. 

This  slice  of  jet  is  further  subdivided  into  nine  annular  segments.  Noise 
radiation  from  each  segment  is  calculated  separately,  and  the  resultant  sound 
power  and  intensity  levels  are  summed  to  obtain  the  overall  noise  radiation 
from  the  jet.  From  experiences  of  previous  numerical  modeling  of  the  noise 
radiation  from  an  entire  exhaust  flow  (Reference  17),  the  total  sound  power  of 
the  entire  jet  is  contributed  mainly  from  the  transition  region  of  the  jet. 
Futhermore,  the  noise  radiation  characteristics  in  the  initial  region  of  the 
jet  is  analytically  similar  to  the  typical  slice  as  chosen  above.  In  most  cases, 
this  typical  slice  can  represent  approximately  70%  of  the  sound  power  from  the 
entire  jet.  Since  a  large  number  of  cases  are  now  considered,  this  simplified 
approach  leads  to  a  tremendous  saving  in  terras  of  computer  time. 
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The  computed  sound  power  is  given  in  Figure  41.  In  this  figure,  only 
relative  decibel  \ulues  are  given.  The  value  of  Mj  is  chosen  to  be  0.6, 

0.8,  1.0,  1.25,  1.50,  and  2.0.  The  speed  of  sound  ratio  A  is  chosen  to  be 
1.0,  1.26,  1.59,  and  2.0.  These  values  of  A  correspond  to  dt.isity  ratios 
of  1.0,  0.63,  0.395,  and  0.25,  respectively.  Since  the  results  in  Reference 
15  are  given  on  logarithmic  scales,  the  above  value  of  density  ratio  are  chosen 
for  convenience.  In  the  SNEGtA-NGTE^  report,  the  sound  power  for  supersonic 
jets  are  normalized  according  to  the  throat  area.  The  sound  power  dependence 
on  f-L  is  slightly  different  from  the  velocity  dependence  laws  which  normalizes 
the  sound  power  according  to  the  jet  exit  diameter.  Figure  42  illustrates  how 
this  prediction  compares  with  recent  GE  measurements.  The  exponential  index 
for  the  density  dependence  can  be  obtained  directly  from  Figure  41  by  measuring 
the  slopes  of  the  curves.  The  results  are  given  in  Figure  43.  The  agreement 
of  the  empirical  modeling  with  the  SNECMA-NGTE  and  GE  results  are  very  close. 

It  is  also  very  interesting  to  note  that  the  predicted  trends  of  the  power  index 
in  the  high  range  are  also  similar  to  the  experimental  results. 

A  normalized  sound  power  curve  has  also  been  obtained  from  the  numerical 
results.  Figure  44  shows  that  the  normalized  sound  power  follows  at  a 
rate  which  is  slightly  above  Uy-law  in  the  high  velocity  range.  This  increase 
in  velocity  dependence  is  caused  partially  by  the  normalization  with  respect 
to  the  nozzle  throat  area. 

In  a  high  velocity  jet,  the  frequency  characteristics  is  controlled  mainly 
by  the  spatial  scale  of  the  turbulence.  This  is  in  direct  contrast  with  low 
speed  jet  noise  where  the  temporal  scale  of  turbulence  controls  the  frequency. 

In  the  present  model,  the  spatial  scale  is  assumed  to  increase  with  A. 

Therefore,  the  entire  noise  spectrum  will  shift  to  lower  frequencies  as  the 
temperature  of  the  jet  increases.  This  is  exactly  the  trend  as  observed  in 
the  SNECMA-NGTE  experiments.  Hence,  with  the  present  assumptions  concerning 
the  turbulence  structure,  the  numerical  results  can  represent  accurately  the 
sound  power  and  the  frequency  characteristics  of  jet  noise  over  a  wide  range 
of  jet  velocity  and  density. 
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2.2  The  effects  of  Refraction  and  Convection 


In  the  Phillips  theory,  the  etfects  of  refraction  and  convection  are 
closely  coupled.  In  addition,  the  aero ! , namical  coupling  between  the  sound 
source  and  the  mean  flow  created  an  internal  attenuation  effect  which  is  not 
detected  in  the  classical  theory  of  aerodynamic  noise.  These  effects  are 
best  demonstrated  by  one  of  the  solutions  of  the  Phillips-Pao  theory: 


„  pr  3/2  2  4  4  4  4 

— 3/2  7!  y  v  M  M  ci 
*  or 

«(y)  =  - 


22, 

r  A  L. 


q. 


— T^r  cos4  6  { (1  -  M  cos9)2+a2M2cos20} 
I  (o )  c  c 


-5/2 


This  equation  is  applicable  for  noise  radiation  in  the  downstream  direc¬ 
tions,  where  the  turbulence  volume  as  responsible  for  noise  emission  is  located 
below  a  transition  point.  In  the  above  equation,  the  attenuation  factor  is 
represented  by  the  function  F(b)/F(0).  The  definitions  of  F(b)  and  its  argument 
b  can  be  found  in  Reference  18.  The  genuine  refraction  effect  is  given  by  the 
ratio  q  /q  .  The  convection  factor  is  different  from  its  classical  form. 

Although  some  refraction  effect  is  represented  by  q^/q^,  the  coupx  ng  of  convec¬ 
tion  and  refraction  is  so  strong  such  that  the  evaluation  of  the  pure  refraction 
factor  alone  becomes  an  academic  excercise.  In  the  present  study,  the  effect  of 
convection-refraction  coupling  is  evaluated  as  a  single  quantity.  If  the 
attenuation  factor  is  removed  from  the  above  equation,  the  remainder  of  the 
equation  is  equivalent  to  the  classical  result.  Two  special  cases  are 
investigated.  The  convection  Mach  number  is  chosen  to  be  1.0  and  2.0,  while 
the  speed  of  sound  ratio  is  1.0.  Results  from  the  Phillips-Pao  theory  are 
compared  with  the  Lighthill  theory.  The  results  oi  the  numerical  calculations 
are  given  in  Figures  45  and  46.  In  the  Mach  number  range  of  interest  to  the 
present  study,  the  coupled  effects  of  convection  and  refraction  reduce  the 
sound  emission  intensity  in  the  downstream  direction  by  a  significant  amount. 

As  a  consequence  of  these  interactions,  the  peak  noise  direction  is  moved 
outward  to  the  45°  to  60°  area.  The  effect  of  attenuation  reduces  further 
the  total  sound  power  in  the  directions  close  to  the  jet  axis.  The  net  effect 
of  refraction  and  convection  can  be  obtained  by  comparing  the  directivities 
of  the  Lighthill  theory  and  the  Phillips-Pao  calculation  without  the  attenuation 
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effect.  The  results  are  given  in  Figures  47  and  48.  The  effect  of  attenua¬ 
tion  is  also  given  on  these  figures.  In  the  downstream  directions,  the  effect 
of  refraction-convection  is  comparable  with  the  attenuation  effect.  However, 
there  is  no  attenuation  effect  in  the  transverse  and  upstream  directions. 

The  coupled  effects  ot  refraction  and  convection  generally  increase  the  noise 
intensity  in  the  upstream  direction.  At  90°  from  the  jet  axis,  the  predictions 
of  the  Lighthill  theory  and  the  Phlllips-Pao  theory  are  identical. 

It  should  be  noted  that  the  classical  convection  effect  is  completelv 
nullified  by  its  coupling  with  refraction  in  upstream  directions.  This  can 
be  seen  clearly  in  the  governing  solution  in  this  region: 


<W(y) 


3  ✓Tit 


3/2  2  4  4  4  4 
Y  v  M  M  a 
o  c 

2  a2  , 

r  A  L, 


(1-M^cosO) 


2  2  2 
(AZ  +  a  M  } 
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The  commonly  recognized  convection  factor  is  reduced  to  a  constant: 


C 


5 


2 

{k  + 


Hence,  the  directivity  in  the  upstream  direction  is  determined  entirely  by 
the  refraction  factor  and  a  Doppler  shift  factor,  as  indicated  in  the  above 
formula.  As  a  consequence  of  this  change  of  mechanism,  the  pre^cted 
directivity  pattern  follows  apparently  a  {(1-M^cose)^  +  o^M^}  law. 

Based  on  these  studies  of  refraction  and  convection,  there  are  several 
possible  experimental  approaches  which  can  serve  to  verify  the  Phillips-Pao 
theory : 

(a)  The  attenuation  effect  is  restricted  to  the  downstream  direction. 

The  amount  of  attenuation  is  a  function  of  frequency.  The  magnitude 
of  attenuation  as  measured  in  dB  is  directly  proportional  to 
frequency.  The  integrated  attenuation  effect  over  the  entire  spectrum 
is  given  by  F(b)  as  previously  shown.  Hence,  the  difference  in 
directivity  for  jet  noise  at  different  frequencies  can  be  compared 
to  the  analytical  predictions. 
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(b)  In  upstream  directions,  the  directivity  pattern  is  determined  entirely 

by  refraction  and  Doppler  effects.  Hence,  the  upstream  directivity 

,  2  2  2,-3/ 2 

appears  to  follow  a  { (1— M  cost))  +  u  Mc }  pattern  instead  of  the 

.  c  0  -5/2 

classical  {(1-Mccos0)^  +  a^M-}  .  Hence,  there  should  be  no  conflict 

c 

between  the  experimental  measurement  and  the  theoretical  predictions 
of  the  convection  effects. 

(C)  In  the  theory,  the  correspondence  between  the  source  and  the  sound 
field  is  precisely  defined.  Hence  correlation  measurements  between 
the  two  should  be  an  effective  way  to  verify  the  theoretical  model. 


2 . 3  Phillips-Pao  Model  for  Detailed  Predictions 

The  purpose  here  as  well  as  in  most  of  the  theory  development  in  this 
program  is  to  formulate  the  theories  in  a  fashion  applicable  for  incorporation 
into  General  Electric's  comprehensive  aero-acoustic  predictive  schemes.  To 
compute  the  sound  pressure  spectra  via  the  Phillips-Pao  theory  the  following 
sound  pressure  spectrum  equations  are  necessary: 


— ?  M,'  ,  f  * 

.  *  —  Tiv  M  L  q  (u>  L.  ) 

(So)  p  ,  ,  o  t  o  V 

■Hr  (r,u)  =  - s-t —  —  - 

p  32r2A6L,  qo  (1-M  cosQ) 


exp  [-|  {(k^)2  +  (wo  Lt)2}] 


tt2v^M^L  A2(o)  tan9  (w  L  cos0)^^ 
(S1)a  V  (r.u.)  -  *  °  C  1  °  1 


8r2A2  LJ/3  H1/3 


(1-Mccos0) 


16/3 


exp  [-  |  {(k1L1)2  +  (co^)2}] 


,2 


(Sl)b  (r ,w)  = 


2  4  8- 

it  v  ML 
o  t 

2  2 

32r  A  L, 


-  -20  <“oLiCOS0) 

-  e  - 5  x 

o  (1-M  cos0) 


exp  [-  ((kjL^)2  +  (wQLt )  2 }  ] 
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\ 

\ 


where 


t  ujcib 

m  -  2rS  (1-M  cosO) 
o  t  c 

k^  =  -MucosO 

k  =  -w  M/A 
o  o 

The  characteristic  Strouhal  number  of  self  noise  and  shear  noise  for  each 
volume  of  turbulence  may  be  computed  via  the  following  expressions: 

For  self  noise 


For  shear-noise 


where 

a^  =  1/8  {(1-M.cos0)2  +  a^M^cos^Q}  L? 

c  c  u 

d  =  MW(r,0)cos0 

M  =  Ui/a 
t  o 

Preliminary  checks  of  this  model  have  been  performed  and  a  complete  check-out 
of  the  method  is  underway.  However,  some  results  can  be  shown.  Comparisons 
of  Pao's  model  with  two  hot  supersonic  jet  measurements  have  been  run.  The 
comparisons  are  made  at  the  important  90°  location.  At  this  location  the 
effects  of  convection  amplification  and  refraction  are  minimal  and  this  location 
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is  a  base  location  for  any  initial  comparison  of  exact  theories.  Figure  49 
shows  a  theory  data  comparison  of  OASPL  for  cases  of  Mj  =  1.095  and  1.464. 
Figure  50  shows  a  set  of  SPL  spectra  comparisons  for  these  two  cases  as  well. 
For  the  case  where  Mj  =  1.095,  close  agreement  between  data  and  computation 
is  observed.  However,  the  comparison  is  rather  poor  for  Mj  =  1.464.  The 
difference  in  SPL  is  more  than  5  dB  for  a  significant  portion  of  the  spectrum, 
whereas  the  computed  spectrum  does  not  demonstrate  the  correct  trend  for 
the  SPL  variation  as  a  function  of  frequency.  It  is  important  to  point  here 
that  the  theoretical  results  as  stated  at  the  beginning  of  this  section 
describes  the  noise  radiation  from  a  single  correlated  volume  of  turbulence. 
Hie  structure  of  turbulence  is  assumed  to  be  Gaussian  and  isotropic  in  a 
moving  frame  of  reference.  For  noise  radiation  from  a  complete  jet  exhaust 
flow,  the  observed  sound  field  is  a  collection  of  contributions  from  many 
individual  source  volume  elements.  A  vastly  complex  synthesis  is  involved 
in  building  a  noise  calculation  model  such  as  the  GE  Aero  Acoustic  prediction 
methods.  In  addition,  a  more  realistic  model  of  turbulence  (vs.  Gaussian 
and  istropic)  may  be  required  for  spectral  predictions.  The  development  of 
this  model  is  currently  performed  by  the  engineering  staff  at  GE,  and 
undoubtedly  it  will  continue  to  receive  great  attention  so  that  this 
predictive  method  will  be  checked  out  completely. 
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3.  THE  QUESTION  OF  CONVECTION  AND  REFRACTION  COUPLING 
IN  JET  NOISE  TURBULENT  MIXING  THEORIES 
H.S.  Ribner 

Lighthill's  theory  of  jet  noise  employs  a  nonconvected  wave  equation 
forced  by  a  certain  source  term.  In  its  exact  form  this  term  accounts  fully 
for  the  co-moving  fluid,  including  convection  and  refraction  effects.  In 
practice  it  is  necessary  to  hold  the  density  factor  constant  in  the  source 
term,  and  this  has  the  effect  of  suppressing  the  refraction  effect  and  throwing 
some  doubt  on  values  of  the  convection  effect.  Alternate  formulations  via 
convected  wave  equations  such  as  are  discussed  in  Section  1  and  2  above  are 
gaining  interest.  An  attempt  is  made  in  this  section  to  bring  some  coherence 
into  the  subject. 

In  addition,  some  supporters  of  the  convected  wave  approach  have  asserted 
that  convection  and  refraction  effects  are  intimately  related  and  must  be 
dealt  with  jointly  in  the  solution  of  the  convected  wave  equation.  On  the 
other  hand,  the  position  taken  here  is  that  convection  and  refraction  can  be 
obtained  as  separable  multiplicative  factors.  In  the  present  work  an  analysis 
is  carried  out  to  indicate  the  degree  of  separation  of  these  two  effects  that 
can  be  obtained. 


3.1  Does  the  Llghthill  Equation  Account  for  Refraction?  Yes  and  No 

Lighthill's  wave  equation  for  an  inviscid  nonheatconducting  fluid  may 
be  written: 


2 

Pv  v 


8x.  x . 

i  1 


+ 


(1) 


On  expanding  the  right-hand  side,  with  use  of  the  equation  of  continuity, 
etc: 


J_£  _ 
at2 


v2P 


l 


2 


c 

o 


JL  ,  n2  3v. 

JLe  _  i_  +  D  —L 

.2  2  2  T  p  ax. 

at  c  Dt  j 


De  2  1  3p 

p  Dt  p  3xi 
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Shifting  terms  from  the  RHS  to  the  LIIS  yields  a  convected  wave  equation 


2  _  8v  9v. 

— ^  -  V  p  -  "extra  terms"  -  p  - —  -r— 

.  2  dx ,  3x , 


by  virtue  of  the  convective  operator 


D_ 
Dt  ’ 


(2) 


Thus  density  derivatives  in  the  expansion  of  Lighthill's  source  term, 
the  RHS  of  (1),  lead  to  a  convected  wave  equation,  hence  to  prediction  of 
refraction  effects  (cf.  work  of  e.g.  Schubert)  arising  from  the  jet  velocity 
gradients. 

In  the  normal  application  of  Lighthill's  equation  (1)  the  density  p  is 
taken  constant  on  the  RHS.  This  has  the  effect  of  eliminating  the  additional 
terms  which  led  to  the  convected  wave  equation  (2);  that  is,  setting  p  ° 
constant  in  effect  suppresses  the  refraction  effect 


3 . 2  Separation  of  Convection  and  Refraction  Effects  by  Use  of  Green's  Function 

We  can  restore  the  missing  refraction  effect  by  solving  (2)  Instead  of 
(1).  Solution  of  (1)  by  the  method  of  Ribner  yields  the  mean  square  pressure 
as  a  product  of  a  "basic  directivity"  times  a  "convection  factor".  The  basic 
solution  of  (2)  yields  a  Green's  function  for  a  stationary  monopole  point 
source.  This  is  then  manipulated  to  yield  a  "basic  directivity"  times  a 
"convection  factor"  -  which  are  the  same  as  for  equation  (1)  -  multiplied 
by  a  "refraction  factor".  The  latter  is  developed  from  the  Green's  function 
according  to  the  following  procedure. 

We,  in  effect,  extract  a  single  frequency  uuO  by  taking  Fourier  transforms 
of  both  sides  of  (2).  Then  we  inject  a  stationary  harmonic  point  source  into 
the  flow  at  £  by  replacing  the  RHS  by  6  (£-&).  The  sound  pressure  p  (x)  is 
the  Green's  function 


sc  (3) 

which  may  be  obtained  by  calculation  (cf.  Schubert)  or  by  the  experimental 
injected  point  source  technique  (cf.  Stvars  et  aL) .  The  refraction  effect  is 
embodied  In  gc  as  it  varies  with  the  direction  of  '•he  observer  vector  x^_ 
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The  solution  of  the  narrow  band  convected  wave  equation  (Fourier  transform 
of  (2))  is  then 


Pc  (x,w) 


_1  r  -  m  f  source  strength  at  yld3y 
♦it  '  ®c  *  [  (narrow  band)  J 


The  source  strength  is  the  F.T.  of  the  KHS  of  (2);  however,  equivalent 
results  can  be  obtained  with  the  Proudman-Lighthill  source  strength  in  narrow 
band, 

a2  2 
d  V 

F.T.  *-2  (5) 

c  3t 
o 

with  p  approximated  by  p  (in  this  case  the  refraction  is  not  suppressed). 

On  squaring  (4)  and  manipulating  the  results  according  to  established 
proceedures  for  random  functions,  there  results 


*  (x,w)  =  }j  g  (x|x';u>)  4>  (y ' ,  y";m)d3y_'d\" 


f  0  J  J  P  V**  |  J  )  w/  T  VJ  »  /  u  i.  / 

c  16tt  c  q 

.2  2 

P  X 

where  4>  is  a  two-point  cross-spectral  density  of  the  source-strength  — 2  — 7 — 

T  „  q  co  3t 

Let 


=  y’+y" 


and  ?  =  -  z.’ 


3  3  3  3 

Then  d  y'd  y"  -*■  d  f;  d  y  in  (6).  If  now,  we  specialize  the  LHS  to  apply 

3 

to  a  unit  volume  at  y,  rather  than  the  entire  jet,  the  integral  over  d  y 
may  be  dropped,  and  there  results 


1  lr--*  3 

(x|x;w)  =  - 2  J  sc  gc  'J’qU.ZJuOdC  at  y 


gc  and  gc*  are  evaluated  to  a  sufficient  approximation  at  rather  than 
y1,  y"  respectively.  This  approximation  reflects  the  known  insensitivity  of 
the  Green's  function  to  source  position  which  was  found  experimentally. 

In  the  far  field  gc  will  decay  like  1/X  and  will  have  a  directional 
factor  and  a  phase: 
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Green's  Function 


S  F(^Y>^)  i* 


(9) 


Thus  (8)  will  have  the  form 


Srctcf  ♦  (»k:»)  ■  /  «  (C.  i.  »)e14Y(T*0li>“) 

Density  vc  -IJL  16t,2x2  '  q  ^  d3r 


(10) 


where  t*  =  £  •  X/CQX,  and  a  weak  dependence  on  azimuth  angle  y  is  neglected 
in  F  and  in  * . 

The  corresponding  result  of  solving  the  ordinary  nonconvected  wave 
equation  of  Lighthill  has  the  form 


Spectral 

Density 


4>0  Oilxiw) 


~ 2  2  /  *  (5>  w)eiwT  d3£ 

16^2x2  « 

phase  factor 


(11) 


Equations  (10)  and  (11)  are  narrow  band  directivities  of  the  jet  noise 
emitted  from  unit  volume  at  y.  For  the  case  of  the  Lighthill  derived  format, 
(11),  was,  in  effect,  postulated  as  the  Fourier  transform  of  a  Gaussian 
2-point  correlation  function  in  earlier  work  of  Ribner  and  of  Ffowcs  Williams. 
(Ribner  used  a  "moving"  correlation  function  in  a  stationary  frame,  as 
required  in  the  present  formulation.  Ffowcs  Williams  used  a  correlation 
function  that  appeared  stationary  in  a  frame  moving  with  the  jet  flow). 

Doth  approaches  gave  identical  results  upon  carrying  out  the  integration 
in  (11).  The  results,  when  normalized  to  unity  at  0  =  90°  take  the  form 


4>0(9»^)  ^  C  ^ (8)  Convection  Factor 

4>o (90°  ,to* )  c“4(90)  (Narrowband) 


where 


C(8) 

u)C(e) 


2  2 

[(1-Mccos0)  +  a  Mc  ]  ;  Doppler  Factor 


o)'C(90°)  =  Source  frequency 


ui  =  observed  frequency 


(12) 


r# 


(13) 


I 
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(Although  (1-Mc  cos  6)  Is  the  true  Doppler  shift  of  frequency  owing  to  source 
motion,  the  constituent  frequencies  of  jet  noise  behave  as  if  there  were  an 
effective  Doppler  shift  C  ^  between  source  frequencies  and  observed  frequencies 
The  reasons  are  fully  explained  in  Ribner  (1960),  (1962),  (1963),  etc.) 

The  corresponding  broad  band  correction  factor  was  obtained  as  C  ^(0)/ 

C  (90°)  by  Ffowcs  Williams  and  by  Ribner. 

The  corresponding  ratio  of  narrow  band  directivities  resulting  from  the 
corrected  wave  equation  is  derived  from  equation  (10) .  It  Is  convenient  to 
multiply  and  divide  by  C  ^(0)/C  ^(90°). 


Correctior 


9/  dB 


NARROWBAND 

DIRECTIVITY 


Monopoie 

Directivity 


CONVECTION 


REFRACTION 


Figure  51 


(Somewhat  elliptic  "Basic  Directivity"  is  overlooked  in  this  analysis). 


•^*^> rswr^ »’r#r*'' , •“wjwf^r  «*<?»*<»«■'«*  jm^+m --» •< 


■,» *x$m 


*wr: 


Thus  Che  net  effect  of  replacing  the  Lighthill  equation  by  the  convected 

wave  equation  is  to  multiply  (or  to  add  on  a  dB  basis)  the  Lighthill -based 

-A 

convection  factor  C  by  a  refraction  factor.  The  refraction  factor  consists 

2  2 

of  a  dominant  part  F  (6,  id),  which  is  the  (amplitude)  of  the  Green's  function, 

2 

and  a  correction  (integral).  The  F  (0,  to)  is  simply  the  directivity  that 
could  be  (and  has  been)  observed  by  injecting  a  monopole  pure  tone  point 
source  into  the  Jet.  The  correction  -  it  is  thought  (cf.  MacGregor,  Ribner, 

Law  (1973)  -  accounts  for  the  difference  between  the  refraction  experienced 
by  a  spherical  initial  pattern  (monopole  pattern)  and  that  which  would  be 
experienced  by  the  nonspherical  pattern  labelled  "Convection"  in  the  sketch. 
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APPENDIX  II- 1 

(E-01,  E-02  and  E-03  denote  10_1,  10_J  and  10" 5 ) 

M  =  0.50 
o 

$  (degrees)  =  10.0  0  (degrees)  =  48.7 


6/1 

M 

arq(R)  (deqrees) 

0. 

4 . 871E-01 

180.00 

0.05 

4 . 856E-01 

178.35 

0.10 

4.81 3E-01 

176.80 

0.18 

4 . 702E-01 

174.76 

4>  (degrees)  =  20.0  6  (degrees)  =  50.3 


i/1 

luL 

arg(R)  (deqrees) 

0. 

2. 094E-01 

180.00 

0.05 

2.081E-01 

177.18 

0.10 

2.048E-01 

174.59 

0.18 

1 . 964E-01 

171.28 

4>  (degrees)  =  30.0  6  (degrees)  =  52.8 


i/1 

ill 

arg(R)  (degrees) 

0. 

5. 301E-02 

180.00 

0.05 

5 . 299E-02 

179.82 

0.10 

5 . 291E-02 

179.66 

0.18 

5 . 274E-02 

179.50 

0  (degrees)  =40.0  8  (degrees)  =  56.4 


6/1 

M 

arq(R)  (deqrees) 

0. 

3 . 276E-02 

0. 

0.05 

3. 253E-02 

-19.09 

0.10 

3. 191E-02 

-37.13 

0.18 

3.056E-02 

-62.16 

<>  (degrees)  =  50.0  0  (degrees)  =  60.9 


6/1 

M 

arg(P)  (deqrees) 

0. 

7.346E-02 

0. 

0.05 

7 . 191E-02 

-15.38 

0.10 

6 . 784E-02 

-29.73 

0.18 

5.912E-02 

-49.21 

94 


M  -0.50  (Cent 
o 

(dl  irc-l’:'.  ) 

-  CO.O 

6  (di  r  ••  .)  -  6  6.4 

<■/> 

JRJ. 

a  i  g  ( E )  ( daqrocs  ) 

0. 

8.  3G5E-02 

0. 

0.05 

8.08] K-cn 

-16.64 

0.10 

7.0108-02 

-31.87 

0.18 

6. 360E-02 

-51 . 80 

(degrees ) 

=70.0  0  (degrees) 

=  73.0 

6/)- 

M 

arq(R)  {deqrees) 

0. 

7.001E-02 

0. 

0.05 

6 . 772E-Q2 

-18.52 

0.10 

6.208E-02 

-35.18 

0.18 

5. 143E-02 

-56.36 

(degrees ) 

=  80.0  0  (degrees) 

=  80.0 

in 

111 

arq(R)  (deqrees) 

0. 

4.042E-02 

0. 

0.05 

3. 888E-02 

-20.51 

0.10 

3.521E-02 

-38.67 

0.18 

2. 869E-02 

-61.20 

(degrees ) 

=  100.0  0  (degrees) 

=  101.0 

6/1 

M 

arq(R)  (degrees) 

0. 

4 . 383E-02 

180.00 

0.05 

4 . 1 86E-02 

156.21 

0.10 

3. 738E-02 

135.56 

C .  18 

3.010E-02 

110.67 

(degrees ) 

=  120.0  0  (degrees) 

=  131.8 

6/1 

M 

arq(R)  (deqrees) 

0. 

6 . 871E-02 

180.00 

0.05 

6. 573E-02 

157.55 

0.10 

5. 887E-02 

137.97 

0.18 

4 . 737E-02 

114.19 
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0.70 


M 


o 


0  (degrees)  -  10.0 


0  (degrees )  =•  54.3 


6/X 

M 

arq(R)  (deqrocs) 

0. 

4.695E-01 

180.00 

0.05 

4.686L-01 

179.20 

0.10 

4.664E-C1 

179.47 

0.18 

4 . 606E-01 

177.56 

(degrees) 

=  20.0  0  (degrees) 

=  55.5 

iLL 

M 

arg(R)  (degrees) 

0. 

1. 847E-01 

180.00 

0.05 

1 . 844E-01 

179.50 

0.10 

1. 838E-01 

179.06 

0.13 

1 . 82  3E-01 

178.54 

4>  (degrees)  =  30.0 


6/X 


0. 

A —  0.05 

\ -  0.10 

0.18 


<{)  (degrees)  =  40.0 


6/X 


0. 

0.05 

0.10 

0.18 


<;  (degrees)  =  50.0 


6/X 

0. 

0.05 

0.10 

0.18 


0  (degrees)  - 


M 

2 . 376E-02 
2 . 626E-02 
3. 206E-02 
4 . 238E-02 


0  (degrees)  = 


M 

6 . 508E-02 
6 . 431E-02 
6 . 224E-02 
5 . 757E-02 


0  (degrees)  = 


M 

1.068E-01 
1.047E-01 
9. 901E-02 
8.663E-02 


57.4 


arg(R)  (degrees) 

180.00 

-162.57 

-152.96 

-149.48 


60.1 


arg(R)  (degrees) 

0. 

-16.25 

-31.72 

-53.63 


63.7 


arg(R)  (degrees) 

0. 

-14.48 

-28.03 

-46.53 
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St 

p 


M  -  0.70 

(Coni . ) 

o 

(deart  s) 

-  00.0 

0  (degrees)  - 

68.3 

6/X 

M 

arq(R)  (degrees) 

0  . 

1. 145K-01 

0. 

0.0  5 

l.  nor.-oi 

-75.71 

0.10 

1.041E-0] 

-30.16 

0.18 

8. 87J.E-02 

-49.20 

(degrees ) 

=  70.0 

0  (degrees)  = 

74.0 

6/X 

M 

arg(R)  (deqrees) 

0. 

9.573E-02 

0. 

0.05 

9 . 276E-02 

-17.68 

0.10 

8 . 5  36E-02 

-33.66 

0.18 

7. 107E-02 

-54.10 

(degrees ) 

O 

o 

00 

It 

9  (degrees)  = 

81.1 

6/1 

M 

arg(R)  (degrees) 

0. 

5 . 571E-02 

0. 

0.05 

5. 366E-02 

-19.96 

0.10 

4. 871E-02 

-37.69 

0.18 

3. 980E-02 

-59.80 

(degrees ) 

=  100. 

0  0  (degrees)  = 

101.4 

6/X 

IlL 

arg(R)  (deqrees) 

0. 

6 . 241E-02 

180.00 

0.05 

5. 949E-02 

155.41 

0.10 

5 . 295E-02 

134.19 

0.18 

4 . 256E-C2 

108.78 

(degrees ) 

=  120.0  0  (degrees) 

=  140.3 

6/X 

M 

arq(R)  (deqrees) 

0. 

6.328E-02 

180.00 

0.05 

5 . 986E-02 

160.17 

0.10 

5.219E-02 

143.77 

0.18 

3. 987E-02 

125.81 

•  — -.v 
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M  =  0.90 
o 


<})  (degrees)  =  10.0 


0  (degrees)  =  58. 5 


S/A 
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III.  EXPERIMENTAL  INVESTIGATION  OF  HEATED  SUPERSONIC  JET 


ACOUSTIC  CHARACTERISTICS 
P.R.  Knott,  K.R.  Bilwakesh,  J.  Brausch 

In  order  to  augment  and  provide  the  theoretical  aero  and  acoustic  efforts 
with  a  body  of  experimental  acoustic  results  at  extended  temperature  ranges, 
a  series  of  detailed  far-field  and  near-field  acoutic.  parametric  studies  were 
performed.  Additionally  a  series  of  photographic  studies  was  also  performed 
for  further  documentation  of  the  exhaust  plume  characteristics.  In  the  para¬ 
graphs  which  follow  some  of  the  preliminary  results  of  these  experiments  will 
be  reported. 

1.0  FAR-FIELD  ACOUSTIC  EXPERIMENTS 

1.1  Experimental  Apparatus.  Test  Set-up  and  Conditions 

The  experiments  were  carried  out  at  General  Electric's  jet  engine  noise 
outdoor  test  site  (JENOTS) .  Sixteen  far-field  microphones  at  a  height  of 
15.93  ft  from  the  ground  were  used  on  a  forty  foot  arc,  at  ten  degree  incre¬ 
ments  starting  from  an  angle  of  20°  with  respect  to  the  jet  axis.  Ten  near¬ 
field  microphones  were  used  to  obtain  representatives  near-field  acoustic 
characteristics  (Section  2  will  discuss  more  about  additional  near  field  tests 
performed).  Figures  52  and  53  schematically  show  the  far-field  and  near 
field  microphone  locations.  The  nozzle  centerline  height  was  55-inches.  The 
first  25  foot  radius  of  the  test  site  is  concrete,  and  the  rest  is  large 
crushed  gravel. 

Acoustic  measurements  were  taken  with  1/4-inch  Bruel  and  Kjaer  condenser 
microphones  (4135) .  The  signals  were  recorded  on  a  twenty-eight  channel  tape 
recorder  (80  KH  flat  response  with  appropriate  corrections),  and  the  microphone 
were  calibrated  with  a  piston  phone  and  oscillator  calibrations  before  each 
test. 


1,2  Nozzles  Tested 

Two  basic  nozzle  configurations  were  tested.  One  nozzle  was  a  convergent  - 
divergent  (C/D)  conical  nozzle  designed  for  parallel  shock  free  flow  at 
Mj  =  1.5  for  stagnation  temperature  operation  between  1500°  R  and  2500°  R. 
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The  second  nozzle  was  a  convergent  conical  thin  lip  nozzle  with  a  1/2  Inch 
thick  lip  adapter.  The  C/D  nozzle  design  was  performed  using  a  method  of 
characteristics  computer  program  with  corrections  for  boundary  layer  displace¬ 
ment.  Both  nozzles  were  water  cooled  and  designed  to  withstand  continuous 
and  non-continuous  testing  at  gas  stream  temperatures  of  3200°  R  at  internal 
total  pressures  to  ambient  pressure  ratios  of  4.0.  The  C/D  nozzle  had  static 
pressure  insi...  ^mentation  provided  for  measuring  wall  static  pressures  along 
the  divergent  section  jf  the  nozzle  and  near  the  exit  plane  of  the  nozzle. 

The  conical  nozzle  had  an  exit  diameter  of  4.3  inches.  The  C/D  nozzle  had  a 

4.3  inch  throat  diameter.  Figures  54  and  55  are  photographs  of  the  nozzles. 

1.3  Experimental  Conditions 

The  experimental  test  conditions  were  so  designed  as  to  examine  the 
velocity  and  temperature  dependence  of  high  velocity  and  high  temperature  Jets. 
Table  I  shows  the  number  of  test  points  and  flow  conditions  taken  for  the  C/D 
nozzle.  Tables  2  and  3  show  the  flow  conditions  tested  for  the  conical  thin 
lip  and  thick  lip  nozzles.  Table  4  shows  the  flow  conditions  tested  for 
shock  free  Jet  exhaust  operation  at  elevated  test  conditions.  The  range  of 
conditions  are  seen  to  be  from  1000  fps  to  3000  fps  each  for  stagnation  tempera- 
tions  ranging  from  1000°  R  to  3000°  R. 

1.4  Preliminary  Results 

1.4.1  Effect  of  Temperature  on  Jet  Noise 

The  work  of  Hoch^  et  al.  first  illustrated  the  influence  of  exhaust  tempera¬ 
ture  on  the  acoustic  characteristic  of  jet  noise.  Following  their  example, 
except  for  an  extended  range  of  conditions.  Figures  56  and  57  illustrate  the 
effects  of  temperature  on  overall  power  level.  Shown  are  results  for  constant 
acoustic  Mach  numbers  ranging  from  .6  to  2.7  for  a  wide  range  of  jet  densities. 
Clearly  shown  is  the  decrease  in  acoustic  power  with  increasing  jet  temperature 
while  keeping  jet  velocity  constant  (the  results  at  MQ  ■  .6  show  an  opposite 

effect).  Also  shown  is  the  deviation  of  density  power  law  dependency  from  the 
2 

classical  p  law  of  Llghthill.  The  solid  lines  drawn  through  the  curves  and 
based  on  Pao’s  theory  are  discussed  in  Section  II-2. 
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TABLE  1 


C/D  NOZZLE  -  FLOW  CONDITIONS  -  ACOUSTIC  TEST  -  Tt/Vj  MATRIX 


Rdg  No 

TPS  Pt  No 

Pt8/Po 

Tt8-°R 

Vj-fps 

1 

1 

1.366 

1081 

1057 

2 

2 

2,045 

1067 

1553 

3 

3 

2.691 

1061 

1795 

4 

4 

3.862 

1052 

2046 

5 

5 

4.513 

1062 

2153 

6 

6 

1.228 

1484 

1011 

7 

7 

1.595 

1469 

1494 

8 

8 

2.417 

1490 

2020 

9 

9 

3.029 

1479 

2227 

10 

10 

4.143 

1474 

2476 

11 

11 

4.516 

1475 

2538 

12 

16 

4.757 

1934 

2947 

13 

15 

2.816 

1904 

2453 

14 

14 

1.891 

1921 

1976 

15 

13 

1.411 

1925 

1478 

16 

12 

1.156 

1925 

970 

1 

17 

1.120 

2362 

954 

2 

18 

1.302 

2308 

1423 

3 

19 

1.633 

2329 

1925 

4 

20 

2.234 

2302 

2407 

5 

21 

3.268 

2330 

2379 

6 

27 

2.643 

2949 

2968 

7 

26 

1.929 

2889 

2458 

8 

25 

1.518 

2857 

1975 

9 

24 

1.313 

2882 

1718 

10 

23 

1.254 

2846 

1468 

11 

22 

1.105 

2885 

987 
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TABLE  2 


CONE.  THIN  LIP  -  FLOW  CONDITIONS  -  ACOUSTIC  TEST  -  Tt/Vj 
MATRIX  +  SHOCK  FREE  DESIGN  LINE  OF  C/D  NOZZLE 


Rdg  No 

TPS  Pt  No 

Pt8/Po 

Tt8-°R 

Vj-f] 

1 

22 

1.109 

2815 

995 

2 

23 

1.243 

3270 

1543 

3 

24 

1.368 

3355 

1865 

4 

25 

1.500 

3316 

2098 

5 

26 

1.949 

3271 

2634 

6 

27 

1.654 

3276 

3134 

7 

21 

3.263 

2576 

3024 

8 

20 

2.231 

2339 

2424 

9 

19 

1.639 

2398 

1960 

10 

18 

1.288 

2375 

1416 

11 

17 

1.126 

2343 

969 

12 

19 

3.829 

1417 

2369 

13 

20 

3.890 

1591 

2522 

14 

21 

3.954 

1755 

2663 

15 

22 

3.990 

1946 

2812 

16 

23 

4.017 

2150 

2961 

17 

24 

4.020 

2309 

3070 

18 

16 

4.797 

1853 

2892 

19 

15 

2.818 

1850 

2418 

20 

9 

3.075 

1409 

2187 

21 

10 

4.120 

1435 

2438 

22 

3 

2.706 

1018 

1762 

23 

8 

2.426 

1444 

1992 

24 

4 

3.888 

1003 

2003 

1 

1 

1.347 

1033 

1010 

2 

2 

2.060 

1019 

1525 

3 

7 

1.661 

1441 

1539 

4 

6 

1.222 

1436 

983 

5 

13 

1.393 

1908 

1445 

6 

12 

1.150 

1883 

942 

7 

24 

1.373 

3000 

1774 

8 

25 

1.488 

2925 

1967 

9 

27 

2.657 

2947 

2973 

10 

15 

2.792 

1876 

2426 

11 

14 

1.891 

1908 

1969 

12 

7 

4.512 

1453 

2519 

13 

5 

4.510 

1018 

2108 
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TABLE  3 


CONE,  THICK  LIP  -  FLOW  CONDITIONS  -  ACOUSTIC  TEST  -  SHOCK  FREE  nESTfifJ  LINE 


Rdg  No 

TPS  Pt  No 

Pt8/Po 

Tt8-°R 

Vj/fps 

4 

19 

3.876 

1411 

2373 

5 

20 

3.916 

1624 

2554 

6 

21 

3.977 

1766 

2676 

7 

22 

4.010 

1946 

2816 

8 

23 

4.035 

2129 

2951 

9 

24 

4.047 

2347 

3102 

14 

19 

3.863 

1447 

2401 

15 

20 

3.903 

1633 

2558 

16 

21 

3.947 

1805 

2699 

17 

22 

4.014 

1979 

2841 

18 

23 

4.019 

2179 

2982 
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TABLE  4 


C/D  NOZZLE  -  FLOW  CONDITIONS  -  ACOUSTIC  TEST  -  SHOCK  FREE  DESIGN  LINE 


Rdg  No 

TPS  Pt  No 

Pt8/Po 

Tt8-°R 

Vj-fps 

12 

19 

3.837 

1415 

2369 

13 

20 

3.925 

1575 

2517 

14 

21 

3.967 

1752 

2664 

15 

22 

4.022 

1927 

2805 

1 

22 

4.022 

1927 

2950 

2 

23 

4.031 

2129 

2950 

3 

24 

4.065 

2326 

3092 
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Figure  58  Illustrates  the  density  dependence  of  peak  OASPL  for  Ideal  exhaust 

speeds  of  1000  fps  to  3000  fps  over  a  vide  range  of  jet  densities.  At  these 

angles,  shock  wave  noise  Is  usually  negligible  and  the  Influence  of  the  jet 

density  on  acoustic  radiation  indicates  that  these  influences  correspond  to  the 

jet's  turbulent  mixing.  The  results  shown  here  indicate  a  somewhat  lesser 

exponent  dependence  than  observed  by  Reference  1;  however,  the  limited  results 
2 

of  recent  Boeing  work  seem  to  be  in  closer  agreement  with  the  data  presented 
here.  For  the  same  range  of  flow  conditions.  Figure  59  shows  a  study  of  the 
Influence  of  jet  density  on  jet  noise  at  6  *  90°.  The  results  shown  are  for 
subsonic  as  well  as  supersonic  jet  exhaust  Mach  numbers.  At  super-critical 
pressure  ratios,  the  noise  levels  may  contain  shock  noise  influence  as  well  as 
turbulent  mixing  noise  characteristics.  These  results  as  well  as  the  velocity 
dependent  nature  of  the  jet  at  other  angles  will  be  studied  further  during  the 
course  of  the  program. 

1.4.2  Velocity  Dependence  of  Jet  Noise 

Preliminary  results  are  shown  only.  The  full  analysis  of  these  results 
and  how  they  compare  with  theory  predictions  will  come  later.  Figure  60 
illustrates  the  velocity  dependent  nature  of  overall  power  level.  The  results 
shown  were  •.  .rmelized  with  respect  to  their  density  power  law  dependence.  A 
velocity  eighth  power  law  dependence  is  shown  as  a  reference  line  in  this 
figure.  Further  resu1  .s  illustrating  the  velocity  dependence  of  high  tempera¬ 
ture  model  scale  jets  are  shown  in  Figures  61  through  68.  These  figures  show 
the  velocity  dependence  of  1/3  octave  band  sound  pressure  for  a  wide  range  of 
Stroahal  numbers:  .1,  .3,  1.,  3.,  5.0,  and  10.0.  Figures  61,  62  and  63  are 
figures  of  the  velocity  dependence  at  0  t  =  150°  at  T*p  *  2000°  R,  2500*  R  and 
3000°  R.  Figures  64  ,  65  and  66  are  for  ®jet  at  90°  at  TT  ■  2000*  R,  2500*  R, 
and  3000°  R.  Figures  67  and  68  illustrate  the  velocity  dependences  per  frequency 

band  f or  6 .  -  40*  for  Tt  ■  2000°  R  and  2500°  R. 

jot  T 

1.4.3  Temperature  Influences  on  Jet  Directivity 

In  subsections  1.4.1  and  1.4.2  results  of  the  parametric  acoustic  test 
results  illustrated  the  velocity  and  temperature  dependent  nature  of  heated 
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(Tg  »  1000*  R  -*•  3000*  R)  air  jets  over  a  velocity  range  of  1000  fps  -►  3000  fps. 
Figures  69  and  70  Illustrate  the  temperature  effect  on  a  jet’s  one-third  octave 
band  directivity  pattern  for  source  Strouhal  number  of  1.0.  Figures  71  and 
72  show  a  similar  directivity  pattern  at  a  constant  frequency  of  3150  Hz  (close 
to  a  Strouhal  number  of  1.0  based  on  the  observed  acoustic  frequency). 

For  the  data  at  a  jet  speed  of  1000  fps  (Figures  69  and  71)  increase  in 
temperature  is  seen  to  move  the  peak  of  noise  to  larger  jet  angles  as  well  as 
to  decrease  the  acoustic  radiated  at  the  small  jet  angles.  The  character  of 
the  directivity  pattern  is  seen  to  be  very  similar  when  the  results  are  plotted 
based  on  source  Strouhal  number  or  observed  Strouhal  number.  For  the  jet 
nozzle  at  an  Ideal  jet  speed  of  2500  fps  (Figures  70  and  72)  the  peak  noise 
stays  constant  at  50*  to  the  jet  axis.  At  the  larger  jet  angles  the  lower 
temperature  tests  show  a  substantial  Increase  in  noise  over  the  higher  tempera¬ 
ture  jets.  These  Increases  may  be  primarily  due  to  increased  shock  noise 
effects.  Further  data  reduction  of  test  results  at  different  Strouhal  numbers 
will  be  pursued  in  future  studies  on  this  program. 

2.0  NEAR  FIELD  ACOUSTIC  EXPERIMENTS 

2.1  Experimental  Apparatus,  Test  Set-up.  Conditions.  Results 

Tests  were  conducted  to  obtain  detailed  near  field  maps  (simultaneously 
with  far  field  measurements)  for  shock-free  and  shocked-flow  hot  jets  for 
source  location  studies.  The  4.3"  throat-diameter  C/D  nozzle  and  the  4.3" 
exit  dia.  thin  lip  conical  nozzle  were  used  respectively  for  the  shock-free 
and  shocked  flow  conditions. 

Simultaneous  with  the  near  field  acoustic  measurement  a  far  field  micro¬ 
phone  array  consisting  of  eight  B&K  4136  microphones  at  a  height  of  15.93  feet 
from  the  ground  on  a  40-ft.  arc  at  angles  of  20",  30“ ,  40*,  50*,  60*,  70*,  80* 
and  90“  to  the  jet  axis  were  also  taken.  The  near  field  mi.ro phone  array 
consisted  of  ten  microphones  (B&K  4136)  mounted  along  a  line  parallel  to  the 
jet  axis  on  a  stand,  which  in  turn  was  mounted  on  a  tracked  cart  capable  of 
being  traversed  in  a  direction  perpendicular  to  the  jet  axis.  The  setup  is 
shown  in  Figure  73  and  74.  This  system  made  possible  obtaining  near  field 
data  at  fifty  microphone  locations,  as  shown  in  Figure  75. 
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The  test  conditions  were  primarily  chosen  to  repeat  the  shock-free  cases 
of  the  C/D  nozzle  for  which  Schlieren  pictures  and  acoustic  far  field  data 
had  been  taken  as  described  in  Section  1.0.  These  detailed  near  field  tests 
covered  a  velocity  range  of  2400-3100  ft.  per  second,  and  in  a  total  tempera¬ 
ture  range  of  1400-2500°  R.  The  same  total  temperature  and  total  pressure 
ratio  conditions  were  then  specified  for  the  thin-lip  conical  nozzle.  A 
correlation  of  the  acoustic  data  with  the  aerodynamic  data  can  thus  be 
expected  to  aid  in  the  source  location  efforts.  In  addition,  a  subsonic 
(Mj  =  O.d)  and  an  Mj  =  1  condition  were  prescribed  for  the  conical  nozzle, 
and  an  Mj  =  1  point  for  the  C/D  nozzle.  Table  5  summarizes  the  principal 
test  conditions. 


Detailed  analysis  is  continuing  on  the  data  obtained  from  these  near 
field  acoustic  measurements.  Plots  of  variation  of  overall  sound  pressure 
level  with  axial  distance  at  fixed  radial  distances  from  the  jet  axis  are 
being  made  for  both  the  conical  and  the  C/D  nozzle  at  various  operating  condi¬ 
tions.  Such  plots  should  be  expected  to  identify  zones  of  increased  acoustic 
output  corresponding  to  localized  sources.  In  addition,  at  supersonic  Mach 
numbers,  comparison  plots  between  a  conical  and  a  shock-free  C/D  nozzle  can 
be  expected  to  shed  light  on  the  influence  of  shocks. 

Sample  plots  of  OASPL  vs.X  (X  «  axial  distance  from  exit  plane)  at  different 

D 

radial  locations  are  shown  in  Figures  76  through  78  for  tht  conical  nozzle 
and  -28  and  -29  for  the  C/D  nozzle.  Figures  76  through  78  correspond  to 
a  constant  temperature  (1500°  R)  and  Increasing  velocity  thus  showing  also 
the  influence  of  velocity  at  constant  temperature  on  the  acoustic  output.  The 
increase  in  OASPL  with  X/D  within  the  limits  of  X/D  shown  appears  to  show  the 
effect  of  the  increasing  turbulence  zone. 

The  influence  of  shock  waves  on  the  acoustic  output  may  be  seen  from 

Figures  81  in  tfiich  OASPL  vs.  X  has  been  plotted  at  different  radial 

D 

locations  for  both  the  conical  and  the  C/D  nozzles  at  ”  1.55  at  the 

same  temperature.  These  results  show  that 

AdB  -  OASPL.  .  .  -  OASPL— 

Conical  CD 

is  considerable,  close  to  the  nozzle  exit  plane  (small  X)  and  drops  off  at 

D 

higher  X  values. 

D 
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In  order  to  examine  this  phenomenon  more  closely,  one-third  octave  band 

spectra  of  SPL  are  described  in  Figures  82(a)  through  82 (j)  corresponding  to 

the  ten  axial  positions  of  the  microphones  at  radial  location  ^  =  4.  These 

spectra  correspond  to  T_  =  2000°  R,  M  -  1.55  and  V.  ideal  =  2800  fps, 

1 8  J 

that  is,  the  same  conditions  represented  in  Figure  81«  It  is  seen  that  at 
the  closer  microphones,  there  is  a  distinctly  higher  broad  band  content  above 
1000  Hz  for  the  conical  nozzle  as  compared  to  the  C-D  nozzle,  but  as  you  move 
downstream  these  effects  are  minimized. 

On  method  of  noise  source  location  is  to  study  the  detailed  near  field 
radiation  field.  To  do  this  contour  plots  of  0ASPL  are  being  obtained  for 
selected  test  conditions.  In  addition,  contour  plots  of  SPL  at  chosen  1/3 
octave  band  frequencies  are  also  being  generated  from  the  data  measured  at 
each  near  field  microphone  locations. 

An  example  contour  plot  of  OASPL  in  the  near-field  region  for  a  shock 
free  flow  is  shown  in  Figure  83  for  a  jet  velocity  of  3100  ft/sec  and  jet 
total  temperature  of  2400°  R.  Comparisons  of  results  are  shown  here  as 
well  as  others,  and  examination  of  the  Schlieren  data  for  corresponding  test 
runs  obtained  earlier  is  also  in  progress  and  will  be  reported  at  a  later  date. 

3 . 0  PHOTOGRAPHIC  &  SOME  SHOCK  NOISE  EXPERIMENTS 

3 . 1  Photographic  Results 

In  addition  to  the  acoustic  tests  discussed  above  a  series  of  Schlieren 
still  and  high  speed  movie  pictures  were  taken  to  further  document  the 
exhaust  plume  charactertist  ics .  Figures  84  through  88  are  photographs  of  the 
Schlieren  still  and  high  speed  movie  systems  used  for  photographic  work 
recently  completed  at  JENOTS.  Figure  84  shows  a  couplets  setup  for  the 
Schlieren  system.  Figures  85  and  86  show  closeups  of  the  Schlieren  still 
and  high  speed  movie  systems.  Figures  87  and  88  show  the  Schlieren  light 
source  systems  for  the  still  and  high  speed  systems. 

Figures  54  and  55  are  photographs  of  the  C/I)  nozzle  and  the  conical 
thick  and  thin  lip  nozzles  for  which  Schlieren  and  high  speed  Schlieren  movies 
were  taken.  Figures  89  through  96  show  typical  still  Schlierens  of  the  C/D 
nozzle  and  the  conical  thick  and  thin  lip  nozzles  operating  at  high  temperature 
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and  velocity  conditions.  Figure  89  shows  the  C/D  nozzle  operating  at 
shock  free  conditions.  Figure  90  and  91  show  the  conical  thin  lip 
nozzle  plume.  Figures  90  shows  the  immediate  vicinity  of  the  nozzle  exit, 
while  Figure  91  shows  the  jet  plume  up  to  the  second  shock  bottle.  Figures 
92  and  92  show  the  exhaust  plume  of  the  conical  thick  lip  nozzle  at  the 
same  operating  conditions  as  shown  for  the  C/D  and  the  conical  thin  lip  nozzles. 
Figures  94,  95,  and  96  show  spark  Schlieren  (previous  photos  were  with 
a  continuous  light  source)  photos  of  the  C/D  and  conical  thick  and  thin  lip 
nozzles.  The  C/D  nozzle  is  operating  at  shock  free  conditions.  A  radiation 
pattern  eminating  from  the  nozzle  lip  is  clearly  visible.  The  frequency  of 
acoustic  radiation  for  these  waves  are  however  around  100,000  Hz,  and  are  not 
considered  to  be  of  acoustic  significance,  but  may  be  indicators  of  jet  insta¬ 
bility. 


3 . 2  Preliminary  Shock  Noise  Results  for  Thin  and  Thick  Lip  Nozzles 

Narrowband  farfield  SPL  spectral  comparisons  of  the  convergent -divergent 
shock  free  jet  with  the  conical  thin  lip  and  thick  lip  nozzles  operating  at  flow 
conditions  of  1400°  F  and  2500  fps  are  shown  in  Figure  97,  98  and  99  for 
jet  angles  of  30°,  130°  and  140°  respectively.  At  near  peak  jet  noise  angles 
the  narrow  band  noise  levels  for  the  three  nozzles  are  nearly  identical.  At 
the  larger  jet  angles  the  conical  thin  and  thick  lip  noz.sles  show  a  shock 
screech  well  above  the  shock  free  C/D  nozzle  noise.  Increased  levels  of  broad 
band  noise  at  frequencies  higher  than  the  screen  tone  are  also  observed.  The 
thick  lip  nozzle  is  seen  to  enhance  shock  screech  over  the  thin  lip  nozzle 
('6  dB  at  the  130°  location)  as  well  as  the  over  all  broad  band  noise  (~5  dB 
over  the  frequency  range  of  1.5  KHz  to  5.0  KHz).  Of  particular  interest  here 
is  the  increased  broad  band  noise  due  to  shock  waves.  Future  work  efforts  will 
be  carried  on  to  study  this  behavior  in  conjuction  with  the  shocked  aero  acoustic 
models  discussed  in  Section  I. 
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TABLE  5  Test  Conditions  for  Detailed  Near  Field  Tests 
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Figure  52  Acoustic  Far-Field  Measuring  Stations 
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Figure  55  4.3  In.  Exit  Diaaeter  Water-Cooled  Conical  Convergent  Nozzle 
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Figure  57  Effect  of  Temperature  on  OAPWL 


146 


1000 


Figure  58  Effect  of  Jet  Density  on  Peak  OASPL 
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Figure  60  Velocity  Dependence  of  Overall  Power  Level 
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Figure  62  Velocity  Dependence  of  l/3  Octave  Band  SFL  syMB( 
at  Source  Strouhal  Nuabera  fron  0.1  to  10.0  j 
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Figure  63  Velocity  Dependence  of  l/3  Octave  Bud  SPL  SYMBOL 
at  Source  Strouhal  Kuabere  frou  0.1  to  10.0  p 
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Figure  64  Velocity  Dependence  of  l/3  Octave  Band  SPL  SYMBOL  fD/V 

at  Source  Strouhal  Suaben  fro*  0.1  to  10.0  Q  0.1 
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Figure  66  Velocity  Dependence  of  1/ 3  Octave  Band  SPL 
at  Source  Strouhal  Numbers  front  0.1  to  10.0 
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Ugure  68  Velocity  Dependence  of  l/3  Octave  Band  SPL  SYMBOL  fD/V 
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1/3  OCTAVE  BAND  SPL,  dB 


SYMBOL 


D  O  <3  O  C 


1/3  OCTAVE  BAND  SPL,  dB 


1/3  OCTAVE  BAND  SPL,  dB 


ANGLE  FROM  JET  AXIS,  DEG. 


Figure  71  Effect  of  Jet  Temperature  on  l/3  Octave  Band 
Directivity  Pattern  at  a  Constant  Source 
Strouhal  Number  of  1.0 
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Figure  72  Effect  of  Jet  Temperature  on  l/3  Octave  Band 
Directivity  Pattern  at  a  Constant  Source 
Strouhal  Number  of  1.0 
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Figure  76  Variation  of  OASPL  with  Axial  Distance  froa  Nozzle  Exit  Plane 

at  Different  Radial  Locations  (4.31  in.  Diameter  Conical  Nozzle) 
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Figure  77  Variation  of  OASPL  with  Axial  Distance  froa  Nozzle  Exit  Plane 

at  Different  Radial  Locations  (4.31  in.  Diaaeter  Conical  Nozzle) 
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Figure  79  Variation  of  OASPL  with  Axial  Distance  fro*  Nozzle  Exit  Plane 

at  Different  Radial  Location*  (4.31  In. Throat  Dlaaeter  C/D  Nozzle) 
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Figure  80  Variation  of  OASPL  with  Axial  Distance  fro*  Nozzle  Exit  Plane 

at  Different  Radial  Location*  (U.31  in.  Throat  Diameter  C/D  Nozzle) 
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Figure  81  Comparison  of  OASPL  vs.  X/D  Profiles  for  Conical  and  C/D 
Jloztles  (M  I.55,  Vj  2800  fps,  T0  2000*R) 
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Figure  02  a-J  l/3  Octave  Band  Spectra  from  Conical  4  C/n  Nozzles 
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Figure  86  JEROTS-Schlieren  n  -Speed"  Movie  System 


Figure  89  Schlieren  Ftaotograph  of  Shock  Free  Nozzle  at  Exit  Plane 
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Figure  91  Schlieren  Photograph  of  Thin  Lip  Conic  Nozzle  Downstream 
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Figure  92  Schlieren  Riotograph  of  Thick  Lip  Conic  Nozzle  at  Exit  Plane 
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Figure  93 


Schlieren  Photograph  of  Thick  Lip  Conic  Nozzle  Downstream 


Figure  95 


Spark  Schlieren  Photograph  of  I 


184 


186 


Figure  98  Far-Field  Barrow  Band  Spectra  C< 
©  •=  130*,  V,  =  2500 fps,  T  -  14 
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IV.  IN-JET  NOISE  SOURCE  LOCATION  STUDIES 
1.0  LASER  VELOCIMETER  DEVELOPMENTS 
P.  Scott,  P.  Mossey,  P.  Knott 


1 . 1  Theoretical  Analysis  for  Turbulent  Spectra 

A  method  has  been  developed  for  determining  the  turbulent  spectra  of 
hot  ^  .J  cold  jets  using  the  LV  device.  Due  to  the  random  arrival  of  the  velo- 
ci ; y  information  in  the  LV  device,  conventional  spectrum  analyzers  could  not 
he  used.  Thus  special  computational  techniques  and  hardware  had  to  be  realized 
before  successful  hot  and  cold  jet  spectra  could  be  obtained.  Future  work  in 
the  area  of  LV  spectra  analysis  consists  of  finalizing  the  analysis  procedures, 
improvement  of  the  LV  processor  to  tailor  it  to  the  requirements  of  spectra 
analysis,  and  construction  of  a  laboratory  LV  spectra  system. 

The  difficulty  in  obtaining  spectra  from  the  LV  device  arrises  from  the 
fact  that  the  velocity  estimate  available  at  its  output  is  not  continuous. 
Successive  estimates  of  the  velocity  are  only  available  when  particles  traverse 
the  probe  volume.  Conventional  spectra  estimation  techniques  assume  that  all 
values  of  the  input  signal  are  known  in  the  analysis  interval.  Such  knowledge 
is  not  available  at  the  LV  output.  The  effects  of  this  are  easily  shown.  If 
we  assume : 

1.  Particles  arrive  at  the  probe  volume  with  equal  probability  (Poisson 
probability  distribution)  in  each  time  interval. 

2.  The  arrival  of  a  particle  at  the  probe  volume  is  statiscally  indepen¬ 
dent  the  value  of  the  jet  velocity  at  the  prtbe  volume,  we  may 
signify  the  LV  output  signal.  y(t),  as  the  product  of  the  time  velocity 
v(t),  and  a  sampling  sequence,  s(t),  consisting  of  impulses  occurring 
at  the  particle  arrival  times. 

Thus  our  model  for  the  LV  output  becomes 
y(t)  *  v (t)  s(t) 
with  auto  correlation  function 


R  (t) 

yy 


R  (t)  R  (t). 
vv  ss 


This  follows  from  the  assumption  of  independence  of  particle  arrival  time  and 
velocity  value. 
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If  the  particles  arrive  in  a  completely  random  manner,  but  at  constant 
rate,  X,  then  the  arrivals  constitute  a  Poisson  process  and  we  may  write 

Rss  (T)  “  Xwo(T)  +  *2 

Using  the  above  and  the  Weiner  relation  which  states  that  the  auto¬ 
correlation  function  and  spectrum  and  wide  sense  stationary  process  are  Fourier 
transform  pairs,  we  obtain 

CO 

S  (a>)  “  X2  S  (u)  +  7  /  S(u)  dp  . 

yy  w  2tt_<o  w 

Thus  the  process  at  the  output  of  the  LV  device  is  the  sum  of  the  original 
velocity  spectrum  and  a  noise  term.  If  the  noise  term  is  large  compared  to  the 
velocity  spectrum  this  spectrum  will  be  obliterated.  The  ratio  of  the  spect.rum 
to  the  noise  term,  D,  is 


V“sw(u) 

From  this,  one  sees  that  there  are  two  ways  to  reduce  the  effects  of  the 

noise  term;  either  analyze  narrow  band  processes  (which  increases 

S  (cj)/  /“ S  (y)  dy)  or  increase  the  mean  sampling  rate  X.  For  spectra 

w  •.  w 

from  know  turbulent  jets,  it  may  be  shown  that  the  mean  sampling  rate  required 
2  6 

is  10  to  10  times  that  achievable  with  the  LV  device.  Thus  it  is  impossible 
to  obtain  spectra  from  turbulent  jets  by  conventional  techniques. 


Method  for  Obtaining  Spectra 


We  are  interested  in  determining  the  spectrum  of  the  jet  velocity  and  not 
the  actual  time  history  of  the  velocity  at  some  point  in  the  jet.  Thus  if 
we  can  determine  a  technique  for  obtaining  spectra  without  reconstruction  of 
the  time  history,  the  problem  associated  with  conventional  techniques  may  be 
circumvented.  A  method  for  doing  this  may  be  based  on  the  fact  that  the  spectr  u121 
of  a  signal  contains  information  about  how  a  point  relates  to  its  neighbors. 

The  "exact"  tine  at  which  this  point  occurs  is  unimport  nt,  an  assumption  which 
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f; 

f 

1 


is  required  if  a  signal  is  to  have  a  unique,  singely  dimensioned  spectrum.  If 
the  jet  velocity  is  sampled  in  such  a  way  that  information  about  the  adjacent 
relationships  is  not  lost  then  there  will  be  sufficient  information  in  the 
sampled  signal  to  reconstruct  the  original  spectrum.  We  will  develop  a 
method  for  performing  the  reconstruction,  stating  in  a  more  rigorous  fashion 
the  constraints  on  the  sampling  sequence.  We  will  then  present  a  model  for 
the  sampling  sequence  corresponding  to  particle  arrivals  in  the  LV  device  and 
show  that  it  has  the  required  properties. 

The  spectrum  estimator  may  be  derivated  one  first  considers  the  auto¬ 
correlation  function  of  the  sampled  process.  If  y(t)  is  the  process  available 
at  the  output  of  the  LV  device,  we  may  model  it  as  the  product  of  a  sampling 
sequence,  s(t),  consisting  of  impulses  at  the  sampling  (particle  arrival)  times 
and  v(t),  the  original  velocity  signal.  We  will  assume  that  v(t)  and  s(t)  are 
statistically  independent  so  that  if 

y(t)  =  s(t)  v(t) 


then 


R  (t)  = 

yy 


R 

ss 


<*> 


where  the  R(r)  ' s  are  the  appropriate  auto  correlation  functions.  This  equation 
suggests  a  method  for  determining  R^Cx),  if  we  write 


R  (t)  = 
w 


R  (x) 

_H _ 

R  (T> 

S3 


We  may  determine  R  (x)  if 

vv 

R  (T)  *  o 
ss 

in  the  region  of  interest.  This  is  the  basic  restriction  on  the  sampling  sequence. 
It  can  be  shown  that  this  is  equivalent  to  stating  that  s(t)  must  contain  samples 
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t  apart  for  all  required  values  of  t.  Knowledge  of  R^Cx)  allows  us  to  deter¬ 
mine  the  velocity  spectrum  SVv  (to)  by  the  Weiner  relation. 


S 

vv 


<U>) 


_  ■  .  -iuT, 

R  \.x)  e  dx 
w 


i.e.  the  auto-correlation  function  and  spectrum  0f  a  signal  are  Fourier  trans¬ 
form  pairs. 

We  are  now  faced  with  the  problem  of  determining  Ryy(x)  and  RgS(x).  To 
faciltate  this,  we  will  change  our  measure  of  x  from  continuous  to  a  finite 
grid,  so  that 


x  ■+>  nT 


where  T,  a  constant,  is  the  grid  spacing.  The  impulses  of  our  sampling  sequence, 
s(t),  will  become  pulses  with  width  T  and  height  1/T  (to  preserve  an  area  of 
one).  These  pulses  will  occur  at  times  nT  with  a  finite  probability,  P.  We 
will  assume  that  T  is  so  small  that  v(t)  is  constant  in  these  Intervals.  Thus 


y(t0)  s  y  (nT) 


for 


(n  -  1/2)  T  <  t  <  (n  +  1/2)  T. 
—  o 


To  estimate  R  (nT) ,  we  will  use  the  conventional 

yy 

*  1  M 

R  (n)  -  Z!  E  y(i  T)  y((i  +  n)  T) 

yy  m  i=1 


where  for  convenience  we  will  assume  the  availability  of  2M  values  of  y(i). 

a 

To  determine  R  (n) ,  we  observe  that 
88 


Ras 


■  E  (s(i)  s(i  +  n)} 
1  1 

-  Z  Z 

s(i)-0  s(i+n)»0 


•W  p{s(1)  .  ilii  ,  s(l+n)  - 

TZ 


or 
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•vawjrfi1.  *’**!'*#?**>  <  .■»-<»  yrn'f  ■'W -v r+r 


>'***^  ^^f*rr*iFT»*T^r>mvmv^rTwrT,’T8W?ty^^ 


R88  (n)  "  ^2  PtS(i)  "  T  »  S(i+n)  "  T}- 

We  will  use  the  last  equality  to  estimate  R  (n)  by  writing 

8S 


R  <n> 

ss 


~2  P  *S(i)  -  £  .  S(i  +  n)  - 


We  now  subsimte  8  e  conventional  estimate  of  P{ S (i)  ■  S(i+n)  ■  ^-}, 


to  obtain 


R  Cn)  -  i?. 

88  t2m 


where  L(n)  is  the  number  of  times  there  is  a  sample  pair  n  apart  after  observing 
M  successive  positions  in  the  record.  We  combine  these  two  estimates  to  obtain 

2  M 

R  (n)  -T  z  y (i)  y(i+n)/L(n) 

W  j  .  i 


In  examining  the  sum  above  we  note  that  M  -  L(n)  of  the  terms  will  be 
zero,  corresponding  to  the  positions  where  s(i)  *  0,  s(i  +  n)  »  0  or  both. 
Thus  we  may  sum  only  the  non-zero  terms  so  that 

2  L(n) 

R  (n)  -  T  E  y(i)  y(i+n)  /  L(n) 

**  i  -  1 


What  the  above  equation  says  is  "sum  all  the  product  pairs  which  occur  in 
the  sampled  data  at  lag  n  and  divide  by  the  number  of  pairs  summed":  i.e., 
determine  the  sample  mean  of  the  v(i)  v(i  +  n)  product  for  each  lag  n.  This 
is  easily  related  to  the  definition  of  the  auto-correlation  function  as  the 
expected  value  of  v(i)  v(i+  n). 

In  determining  the  spectrum  of  a  turbolet  jet  we  use  the  above  to  estimate 

the  auto-correlation  function  on  a  fine  grid.  We  then  use  the  techniques  of 
1  2 

Blackman  &  Tukey  and  Parzen  to  take  the  weighted,  discrete  Fourier 
transform  to  obtain  the  spectrum  estimate. 
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1.3  Error  Analysis 


In  the  previous  section  we  have  developed  a  method  for  estimating  the 
auto-correlation  function  of  velocity  turbulence  data  obtained  with  a  laser 
velocimeter.  This  estimate,  under  the  assumption  of  an  infinitely  fine  time 
grid,  is  bias  free,l*e*: 


M 

I 


1X1 

M 


R 


xx 


(n) 


The  introduction  of  a  finite  width  time  grid,  however,  introduces  two 
bias  errors.  The  reason  for  these  errors  may  be  seen  by  examining  Figure  100 . 
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Figure  100 


2T  3T  4T 

3.5T  4.5T 


5T 


6T 


If  the  time  between  arrivals  is  measured  in  discrete  intervals  of  length 
T,  then  our  estimate  of  really  consists  of  the  average  of  lag  products 

x(t)x(t+NT+ro)  where  t0  is  a  random  variable  on  interval 


-T/2  <  t  <  T/2 

—  o 


These  "jitters"  in  the  lag  of  each  product  pair  causes  a  bias  error  whose  effect 

3 

may  be  estimated  using  the  techniques  of  Balakrishnan.  This  analysis  appears 
as  Appendix  IV-1.  It  is  shown  that  if  the  jitters  have  nearly  the  same  density 
function  then: 

1.  The  spectrum  of  weighted  by  the  characteristics  function  of  the 
jitter's  probability  density. 

2.  This  weighted  spectrum  will  then  be  aliased  about  frequency  which  is 
the  reciprocal  of  the  grid  spacing. 
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3.  Differences  in  the  jitter  density  function  from  interval  to  interval 
cause  additional  biases  to  be  added  to  the  spectrum  these  terms 
disappear  at  f  -  o,  and  increase  with  frequency. 

Experiment  shows  that  the  grid  spacing  for  turbulent  velocity  spectra 
measurement  should  be  at  least  a  fifth  of  what  would  be  chosen  as  the  Nyqulst 
spading  in  convention  "sample  data"  analysis.  In  the  system  used  at  General 
Electric  a  factor  of  one  tenth  was  used. 

One  may  wonder  why  an  equally  spaced  grid  used  if  it  introduces  these 
errors.  We  use  it  because: 

1)  It  is  required  to  get  a  finite  probability  of  arrival  for  products 
at  a  given  lag  v^lue. 

2)  It  allows  the  use  of  the  FFT  in  performing  the  spectrum  calculation 
which  makes  this  computation  tractable  In  a  small  computer. 

3)  It  allows  the  use  of  all  of  the  arriving  product  pairs. 

We  will  now  consider  the  other  parameters  affecting  the  spectrum  estimate. 
It  is  useful,  since  aliasing  occurs  about  the  grid  Interval,  to  consider 
the  spectrum  in  the  Z-domain  here  the  grid  interval  is  used  as  the  sampling 
interval.  If  we  define  f^  as  the  frequency  above  which  there  is  no  significant 
power,  then  at  one  tenth  grid  spacing,  the  grid  width,  T,  should  be: 


For  turbulence  spectra  f^  lo  assumed  to  be  20  kHz  so  that  for  our  system 

T  -  2.5  x  10-6  sec. 

On  the  unit  circle  in  the  Z-domain,  we  have: 
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Figure  101 

Hence  we  are  only  interested  in  determining  the  spectrum  on  the  corres¬ 
ponding  to  the  pie  shaped  shaded  region  in  Figure  101.  Let  us  now  consider  the 
bandwidth  of  the  analysis.  This  is  determined  by  the  length  (total  number  of  lags 

computed)  of  the  auto-correlation  function  (a.c.f.)  and  the  "window"  used  in 

1  2 

the  spectrum  estimate.  The  window  function  *  is  used  to  reduce  undesirable 

effects  induced  by  the  finite  a.c.f.  length.  The  window  function  used  in  the 

2 

General  Electric  system  is  attributed  to  Parzen.  The  bandwidth  of  the  spectrum 

3 

analysis  may  to  shown  to  be 


B  W  ^  dB 


V 

(m  -  1)T 


where  p  is  a  function  of  the  window  (p  'v  1.4  for  Parzen)  j  and  m  is  the  number 
of  lags  computed.  If  a  bandwidth  of  about  100  Hz  is  required  then  a  value  of 
m  -  5120  -  512  x  10  may  be  used.  Thus  we  must  calculate  5120  averaged  lag 
products  to  obtain  a  spectrum  with  100  Hz.  Band  width.  The  spectrum  may  be 
calculated  by  taking  the  discrete  Fourier  transform  of  the  weighted  auto-correla¬ 
tion  function.  One  way  of  doing  this  would  be  to  use  an  FFT  routine  on  a  5120 
point  block.  Computation  time  and  computer  storage  may  be  reduced  if  we  compute 
the  transform  output  points  only  over  the  segment  shaded  in  Figure  10L.  l't 
may  easily  shown  that: 

M-l 

S(u)  -  T  E 

1-0 

-F1  (2  Rn  ((£  +  kN))T)»  Sin  (tutT)l -TR^O) 


[FR  (2  R  ((£  +  kN)T)»Cos  (utT) 
xx 
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where  F  and  F  are  the  real  and  imaginary  parts  of  an  m  point  FFT  routine.  In 
our  case,  N  =  10,  M  ■  512  and  to(J)  =  2ti(J)/NWT  where  T  is  the  grid  spacing. 

The  auto-correlation  function  and  spectra  equations  were  implemented 
in  a  combination  of  machine  language  and  Fortran  on  a  PDP8/E  mini- computer. 

This  machine  is  equipped  with  8K  bytes  of  core  memory  and  two  Dectape  magnetic 
tapa.  drives.  A  TZ  digital  interface  was  designed  to  link  the  laser  velocimeter 
processor  to  the  computer.  This  Interface  also  measures  the  arrival  time  be¬ 
tween  successive  particles  and  provides  sixteen  words  of  buffer  memory.  This 
high  speed  memory  is  required  because  the  arrival  time  between  successive 
particles  can  be  faster  than  the  computer's  ability  to  absorb  the  data  (even 
though  the  computer  mean  data  read  rate  is  considerably  faster  than  the 
particle  arrival  rate) .  Analysis  shows  that  the  16  word  memory  should  over¬ 
flow  less  than  one  in  one  billion  arrivals  if  the  mean  particle  arrival  rate 
is  less  that  20  KHZ.  A  display  scope  was  also  interfaced  to  the  computer  and 
software  was  written  to  allow  the  display  of  the  turbulence  spectrum  and  auto¬ 
correlation  function  together  with  the  sampling  auto- correlation  function, 
and  the  interarrival  and  velocity  amplitude  histograms. 


1.4  Experimental  Set-Up  Used  for  Laser  Velocimeter  Turbulent  Spectra 
and  RMS  Turbulent  Plume  Surveys 

1.4.1  The  Laser  Velocimeter  Arrangement 

The  arrangement  used  was  similar  to  that  of  the  1971  A.F.  Program  as 

4 

reported  in  AFAPL-;TR-72-52  Chapter  V  .  The  details  of  that  report  will  not  be 
repeated  here.  Instead,  changes  or  new  features  and  improvements  will  be 
described  that  have  been  used  in  this  current  reporting  period.  The  basic 
optics  system  used  was  again  the  differential  Doppler,  backscatter,  single 
package  arrangement  that  has  the  proven  feature  of  ruggedness  for  our  rather 
severe  environment.  A  substantial  improvement  in  temperature  stability  of 
the  optics  alignment  was  made  (on  an  internal  G.E.  Program)  that  also  allowed 
greater  range  without  sacrificing  special  resolution  of  the  measurement  volume. 
Figure  102  shews  a  schematic  arrangement  of  the  laser  package  used  on  this 
program.  The  laser  beams  no  longer  are  coaxial  with  the  receiving  lens,  but 
instead  are  projected  from  belcw  the  lens,  forming  an  angle,  a,  that  keeps 
the  major  axis  of  the  control  volume  ellipsoid  to  a  minimum.  Compare  this 
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figure  to  Figure  V-B-4  (page  550)  of  the  AFAPL-TR-72-52.  The  dimensions 
of  the  control  volume  are:  major  axis  0.25  inch,  minor  axis  0.020  inch.  The 
range  has  been  extended  from  40  inches  to  85  inches  (fixed,  as  before).  The 
extension  has  allowed  backing  off  the  package  by  45  inches  from  the  jet  and  the 
acoustic  protection  enclosure  as  used  in  1971  is  no  longer  necessary.  The 
three  steering  mirrors  and  the  beam  splitter  were  remounted  on  adjustable 
supports,  all  the  same  aluminum  alloy,  and  this  eliminated  the  temperature- 
alignment  problem  seen  in  1971  when  we  used  some  brass  and  stainless  steel 
pieces.  Figure  103  shows  the  Laser  Velocimeter  setup  at  GE's  JENOTS  facility. 

The  range  extension  to  85"  brought  about  two  problems,  however.  The  first 
of  these  is  the  Increased  sensitivity  to  temperature  gradients  in  the  air 
intervening  between  the  package  and  the  sensitive  volume.  This  has  caused  a 
small  decrease  in  accepted  data  rate  in  cold  weather  where  convection  gradients 
occur  right  at  the  front  of  the  package,  where  the  laser  beams  emerge. 

The  data  rate  was  estimated  to  drop  by  20%  to  30%  when  the  outdoor  tempera¬ 
ture  was  25°  F.  The  package  temperature  is  usually  held  to  60  to  80°  F  by 
a  thermostat. 

The  second  problem  encountered  in  the  range  extension  is  that  of  a  reduc¬ 
tion  in  the  return  ligjht  power.  This  reduction  in  return  light  power  is  because 
a  6-inch  diameter  lens  is  being  used  at  a  greater  distance  (smaller  solid  angle 
of  light  collection).  This  reduction  in  light  power  collected  was  partly  com¬ 
pensated  for  by  increasing  the  electrical  gain  of  the  photomultiplier  through  use 
of  higher  supply  voltage.  This  resulted,  however,  in  greater  shot  (light 
quantum)  noise  reaching  the  LV  processor.  As  a  result,  the  percentage  of  laser 
Doppler  bursts  validated  dropped  by  perhaps  a  factor  of  two. 

1.4.2  Actuator  and  Seeding 

A  remotely  actuated  platform  was  again  used  but  this  time  on  all  three 
axes:  vertical,  horizontal  and  axial.  Travel  capabilities  were  32,  32,  and 
240  inches,  respectively,  all  remote  read-out.  Resolution  was  +1/16  inch  for 
each  axis,  except  the  last  208  Inches  of  the  axial:  +  1/8  inch. 


198 


Seeding  was  by  injection  of  AL^  powder,  nominal  1  micron  diameter,  into 
the  supply  air  to  the  burner  and  also  into  the  region  of  the  nozzle,  so  as  to 
seed  the  entrained  air.  The  power  feeder  equipment  used  was  as  reported  in 
AFAPL-TR-72-52  Chapter  V  section  3,  except  that  the  fluidized  bed  column  supply 
air  was  heated  to  about  250°  F  to  prevent  power  aggregation  by  moisture  absorption. 

1.4.3  Signal  Processing  and  Recording 

The  laser  velocimeter  signal  processor  used  was  the  G.E.  developed  direct- 
counter  (time  domain)  type  similar  to  that  reported  in  AFAPL-TR-72-52  Chapter 
V,  but  with  some  improvements  made  during  1973  on  a  G.E.  internal  program. 

These  improvements  resulted  in  a  lowered  rate  of  false  validations  and  Improved 
linearity  and  resolution.  The  recording  of  turbulent  velocity  probability 
distributions  (histograms)  was  again  by  the  NS633  Pulse  Height  Analyzer,  256 
Channel,  dumped  into  an  X-Y  plotter.  The  recording  of  spectra  was  by  the 
same  equipment  as  used  as  discussed  in  Subsection  1.5. 

1.5  Laser  Velocimeter  Turbulent  Spectra  Measurements 

In  the  above  subsection  descriptions  of  the  analysis  and  laser  equipment 
used  for  estimating  turbulent  spectra  were  given.  Descri7  «*d  here  are  a  set  of 
demonstration  experiments  performed  to  measure  the  longitudinal  component  of 
turbulence  spectra  in  the  exhaust  plume  of  a  4.3  inch  diameter  convergent  noz¬ 
zle  (See  Figure  54)  at  General  Electrics  JENOTS  facility.  Experiments 
were  performed  on  a  cold  subsonic  jet  and  in  a  sonic  heated  jet.  These  tests 
were  performed  in  late  December  of  1973  and  are  preliminary. 

1.5.1  Cold  Jet  LV  Turbulence  Spectra  Results 

The  conical  nozzle  was  set  at  M,  »  0.5,  T„  ■  ambient.  The  location  of 

J  ° 

the  laser  velocimeter  measuring  volume  was  at  an  x/d  *  6,  r/rQ  *  1.  Figure 
104  shows  the  measured  autocorrelation  function  and  the  interparticle-arrival 
histogram.  Figure  105  shows  the  constructed  histogram.  The  LV  estimated 
turbulence  spectrum  for  this  jet  using  the  analysis  of  Section  1.1  is  shown  on 
Figure  106.  Also  shown  on  Figure  106  is  a  hot  film  spectrum  taken  on  the 
same  jet  at  the  same  location.  The  agreement  between  the  Laser  Velocimeter 


199 


results  analyzed  to-date  were  for  a  limited  data  sample.  The  limited  data 
sample  is  expected  to  give  a  larger  variance  at  the  high  frequencies.  To 
reduce  the  variance  at  the  high  frequencies,  more  data  product  pairs  of  the 
data  samples  will  be  processed. 


1.5.2  Hot  Jet  LV  Turbulence  Spectra  Results 

Figures  l07  to  109  illustrate  LV  measurements  for  a  sonic  jet  (Vj~1500 
fps)  with  total  temperature  at  15008  R.  The  LV  was  set  at  X/D  -  10,  r/rQ  ■  0, 
and  was  measuring  the  longitudinal  component  of  turbulence.  Figure  107  shows 
the  constructed  autocorrelation  function  for  the  heated  high  velocity  jet. 

Figure  108  shows  the  constructed  histogram,  and  Figure  109  shows  the  LV 
measured  turbulence  velocity  spectrum.  There  are  no  other  measuring  devices 
that  can  measure  the  turbulence  velocity  spectra  in  this  type  of  environ¬ 
ment,  and  therefore,  no  comparisons  will  be  available.  As  was  Indicated  with 
the  cold  jet  turbulent  measurements  the  higher  frequency  data  is  expected  to 
be  in  some  error  due  to  the  limited  data  sample  analyzed.  Currently  more  data 
product  pairs  of  the  data  samples  are  being  processed  to  reduce  the  variance 
at  the  high  frequencies  and  these  results  will  be  reported  on  at  a  later  date. 

The  above  results  demonstrate  the  feasibility  of  using  the  laser  velocimeter 
to  construct  turbulent  veloicty  spectra  in  realistic  exhaust  jet  environments. 

In  addition  to  performing  measurements  on  an  ambient  temperature  *»  0.5  jet, 
turbulent  velocity  measurements  were  performed  on  a  high  temperature,  high  vel¬ 
ocity  jet.  The  laser  velocimeter  system  used  to  perform  the  measurements  is 
one  which  can  be  used  for  small  or  for  the  large  model  scale  tests  as  per¬ 
formed  above,  or  on  an  actual  jet  engine  nozzle.  These  preliminary  tests  are 
a  prelude  to  demonstration  experiments  planned  for  in-jet  noise  source  location 
studies  at  hot  supersonic  exhaust  jets. 


1.6  Laser  Velocimeter  Jet  Plume  Surveys 

As  part  of  the  development  and  evaluation  of  the  laser  velocimeter  as  a 
useful  non-contact  aero-acoustic  probe,  extensive  surveys  of  the  mean  velocity 
and  turbulent  rms  velocity  characteristics  were  performed  on  subsonic  and  sup¬ 
ersonic  shock  free  and  shocked  hot  jet  exhaust  plumes.  Nearly  1000  data  points 
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were  taken.  At  the  time  of  this  writing  the  bulk  of  this  data  is  being  processed 
for  data  analysis.  Reported  here  are  some  comparisons  of  LV  measured  flow 
properties  with  measurements  taken  with  a  hot  film  anemometer  for  a  Mach  0.5 
cold  jet.  Preliminary  LV  results  are  also  shown  illustrating  the  mean  and 
turbulent  flow  characteristics  of  a  shock-free  (M^  =  1.55)  heated  (T^  -1500°  R) 
jet. 

The  nozzles  tested  were  the  4.3-inch  convergent  conical  nozzle  (See  Figure 
54)  and  the  4.3  inch  throat  diameter  convergent-divergent  nozzle  (Figure 
55.  The  laser  velocimeter  system  used  the  one  described  in  subsection  1.4. 

The  comparison  hot  film  anemometer  probe  for  the  cold  jet  test  was  a  Thermo  Systems 
model  1210.  Figure  110  shows  the  arrangement  used  for  the  hot  film  measure¬ 
ments.  An  extention  beam  was  attached  to  the  automatic  LV  traversing  cart  so 
that  accuracy  of  probe  location  could  be  obtained. 


1.6.1  Cold  Jet  Laser  Velocimeter/Hot  Film  Measurements 

1.6. 1.1  Axial  Variation  of  Mean  Velocity 


Figure  111  shows  the  axial  variation  distribution  of  the  centerline  mean 
velocity  as  measured  by  the  laser  velocimeter  and  the  hot  film  anemometer.  For 
reference  purposes  hot  wire  measurements  of  Wooldridge  and  Wooten"*  at  M^  -  0.3 
are  shown  as  a  solid  triangle.  The  LV/Hot  film/hot  wire  comparisons  as  seen 
to  be  in  good  agreement.  The  results  show  that  the  mean  velocity  in  the  core 
began  to  decrease  at  an  axial  location  of  X/R  ■  6.  Between  X/D  ■  6  and  15  the 
flow  undergoes  a  complicated  readjustment  toward  that  of  a  self  slmiliar  region. 


1.6, 1.2  Radial  Variation  of  Mean  Velocity 

A  mean  velocity  profile  comparison  between  the  laser  velocimeter  and  the 
hot  film  measurements  is  shown  in  Figure  112.  The  radial  profile  shown  is 
at  an  axial  of  X/D  ■  2.  This  station  was  choosen  so  that  the  results  of 
Wooldridge  and  Wooten  could  also  be  use  as  a  comparison.  Good  agreement  be¬ 
tween  the  LV,  hot  film,  and  hot  wire  measurements  is  observed. 
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1.6. 1.3  Axial  Variation  of  Turbulence  Intensity 


Figure  113  shows  comparisons  between  laser  velocimeter  measurements  and 
hot  film  anemometer  measurements  for  axial  variations  of  turbulence  intensity 
along  the  jet  centerline  and  at  the  jet  lip  (r/D  *  .5).  Wooldridge  and 
Wooten's  measurements  are  again  shown  for  comparison  purposes.  The  agreement 
between  all  three  instruments  is  again  found  to  be  good.  The  centerline  var¬ 
iation  of  turbulence  intensity  is  seen  to  vary  smoothly  from  a  very  low  value 
ct  the  exit  to  a  peak  value  around  an  X/D  of  about  10.  The  lip  region  var- 
i*Mon  of  turbulence  intensity  is  seen  to  rapidly  Increase  from  a  low  value 
at  the  exit  plane  to  a  rather  large  value  at  X/D-2.  From  X/D~2->8  the  intensity 
levels  appear  rather  uniform. 

1.6.2  Hot  Superronic  Jet  Laser  Velocimeter  Measurements 

1.6. 2.1  Axial  Variation  of  Mean  Velocity 

The  axial  variation  distribution  of  the  centerline  mean  velocity  is  shown 
in  Figure  114.  A  uniform  distribution  is  observed  until  about  10  diameters 
downstream.  Thereafter  the  usual  axial  velocity  decay  is  observed. 

These  results  are  in  agreement  with  what  would  be  expected  for  a  supersonic 
shock-free  exhaust  jet. 

1.6. 2. 2  Radial  Variation  of  Mean  Velocity  and  Turbulent  Intensity 

Figure  115  illustrates  the  LV  measured  radial  distribution  of  mean 
velocity  and  turbulent  intensity  of  the  Mj  “  1.55,  Tj  -15000  R  shock  free 
exhaust  at  an  axial  location  of  X/D  =•  9.6.  The  mean  velocity  distribution  is 
as  expected  and  the  turbulence  intensity  profile  shows  a  peak  in  turbulence 
intensity  near  the  nozzle  lip  region  (r/D'. 5). 

1.6.2. 3  Axial  Variation  of  Turbulence  Intensity 

The  axial  variations  of  the  supersonic  jet  turbulence  intensity  along 
the  jet  centerline  and  the  nozzle  lip  is  shown  in  Figure  116.  The  centerline 
distribution  slows  a  very  gradual  increase  in  the  turbulence  intensity  level 
out  to  about  ter  diameters.  Past  the  10  diameter  location  an  abrupt  increase 
in  level  is  observed.  The  nozzle  lip  axial  variation  of  the  turbulence  level 
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is  also  seen  to  be  more  gradual  than  was  observed  for  the  cold  subsonic  jet 
results  (Figure  113).  A  peak  value  is  observed  to  occur  between  10  to  12 
diameters  downstream. 


The  above  comparison  between  General  Electric's  Laser  Velocimeter,  the  hot  film 
anamometer,  and  other  published  results  have  shown  good  agreement  for  mean  velocity 
and  turbulent  intensity  variations.  Preliminary  laser  velocimeter  measurements  of 
mean  velocity  and  turbulent  intensity  have  been  shown  for  a  shock  free  (M^  -1.55) 
high  temperature  (T^  -1500°  R)  jet.  These  results  are  a  first  of  a  kind  type 
of  measurements  and  more  of  the  measurements  that  were  taken  on  this  jet  will 
be  reported  at  a  future  date.  A  number  of  other  LV  measurements  have  been 
obtained  on  cold  subsonic  and  heated  supersonic  exhaust  jet  over  a  varied  range 
of  velocities  (1000  3000  fps)  and  temperatures  (1000°  R  -*•  3000°  R) .  These 

results  are  currently  being  reduced  and  analyzed  for  future  reporting.  The  above 
results  and  the  turbulence  spectrum  results  discussed  in  section  1.4  clearly 
demonstrate  the  General  Electric's  Laser  Velocimeter  as  a  viable  non-contact 
type  of  measurement  device  capable  of  measuring  the  needed  detailed  fluid  dynamic 
turbulent  flow  characteristics  of  high  velocity  high  temperature  jets  necessary 
for  performing  noise  source  location  type  of  experiments. 
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2.  STATIC  PRESSURE  PROBE  DEVELOPMENT  FOR  JET 
NOISE  SOURCE  CORRELATION  STUDIES 


R.A.  Kantola,  J.F.C.  Wang 

Recent  research  activity  in  the  area  of  jet  noise  has  been  heavily  involved 
in  cross-correlation  measurements.  These  have  included  correlations  of  the 
unsteady  pressure  and  velocities  in  the  jet  plume  to  the  far-field  acoustic 
signal  and  in-jet  two  point  correlations  of  pressure,  and  pressure  and  velocity. 
This  work  is  aimed  at  extending  the  techniques  to  high  speed  flows,  particularly 
the  in-jet  correlations  and  resolving  the  discrepancies  reported  in  earlier 
investigations. 


2. 1  Theoretical  Background 

The  analytical  basis  for  the  correlation  technique  has  been  well  documented 
and  will  only  be  briefly  presented  here  as  a  point  of  reference.  The  starting 
point  for  this  is  the  "dilatation  equation",  which  has  been  shown  to  be  equiva¬ 
lent  to  Llghthill's  original  formulation  using  acoustic  quadrupoles.  This 
dilatation  equation  is  given  as: 


i?£(1) 

at2 


+  V  p 


2  (1)  _  _1_ 
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3t2 
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where  p  ^  -  propagating  part  of  the  pressure  and 
portion  of  the  unsteady  jet  pressure,  p^\  is  defined 


the  non-propagating 
as 


P 


(0) 


(2) 


where  p  *  instantaneous  static  pressure 
and  p0  «  average  static  pressure 

The  right  side  of  equation  1  is  viewed  as  a  forcing  function  of  the  wave 
equation  for  the  acoustic  pressure  p^.  This  wave  equation  assumes  that  the 
propagation  is  taking  place  in  still  air,  not  in  a  jet  flow,  and,  therefore, 
the  effects  of  refraction  and  convection  are  not  accounted  for.  For  farfield 
microphones  sufficiently  far  from  the  sources  and  for  the  origin  of  the 
coordinate  system  sufficiently  close  to  the  sources,  equation  1  can  be  solved 
in  terms  of  Kirchoff's  retarded  potentials  as 
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'  volume 

where  X  =  microphone  arc 

Cc  =  isentropic  wavespeed 
y  =  source  location 


dy3 


(3) 


t  *  retarded  time  ■  t  - 


Now  assuming  that  the  far-field  fluctuating  static  pressure  p(x,t)  Is 
represented  by  p^(x,  t)  and  then  multiplying  both  sides  of  equation  2  by 
p(x,t)  evaluated  at  a  new  time  t  +  t,  the  result  can  then  be  time-averaged  to 
yield 


p(x, t )p(x, t+t) 


[p‘ V0^2.,t)J*  p(x,t  +t)d^J 


volume 


(4) 


With  statistically  stationary  variables  only  the  time  delay  t  is  important 
and  results  in  the  equation  4  reduce  to: 


_  -1 

p(x)p(x,t)  ™  4^fCo2x 


[p‘(°)(^)p(x,t)]  ?  (d£3) 

volume  * 
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This  equation  says  that  the  rms  value  of  the  far-field  acoustic  signal 

a 

can  be  found  by  setting  t  ■  0,  and  evaluating  the  integral  at  t  ■  x/C0  and 
then  taking  the  square  root.  This  will  give  the  far-field  acoustic  pressure 
at  one  point  due  to  all  the  sources  in  the  jet  plume.  The  time  derivative 
shown  here  can  be  represented  for  stationary  random  variables,  in  other  forms, 
as: 

~2_  ^  *~2  _____ 

p(o)p  (?)  "  2^  p(o)  P  (t)  P  (o)  p(?)  (6) 


This  development  follows  very  closely  that  given  in  Reference  6. 
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2.2.  PROBE  SELECTION 


2.2.1  In- Jet  Pressure 

2. 2. 1.1  Pressure  Probe  Selection  Requirements 

To  adequately  measure  the  fluctuating  pressure  within  the  Jet  plume  and 
not  cause  extraneous  readings  in  the  in-jet  pressure  measurement  and  more 
critically  in  any  correlation  of  the  pressure  signal  to  far-field  measurement 
demands  that  the  probe  selection  be  done  very  carefully.  Of  course  the  inser¬ 
tion  of  any  probe  will  cause  some  disturbance  to  the  jet  plume  and  radiate 
some  sound,  however,  this  must  be  minimized  by  the  proper  selection  of  probe 
size,  jet  diameter  and  probe  type.  For  cold  jets  the  criteria  can  be  given 
as  follows: 

o  Probe  Frequency  Response 

o  Probe  Contamination  Factor 

o  Repeatability 

A.  Probe  Frequency  Response 

Probe  amplitude  and  phase  response  must  be  reasonably  flat  to  frequencies 
greater  than  the  signal  content  in  the  jet  plume  and  the  far-field.  Fuchs3 
has  measured  the  fluctuating  pressure  in  the  jet  and  found  that  at  3  diameters 
downstream,  on  the  jet  centerline,  the  power  spectral  density  is  down  by  40 
dB  at  a  Strouhal  number,  Sj,  of  2.5  compared  to  the  peak  value  at  Sj«0.45. 

Rakl^  and  Siddon7  made  measurements  of  the  in-jet  pressure  field  with  a  small 
alrfoll-llke  probe  and  these  indicate  that  the  major  portion  of  the  correlated 
"pseudo  sound"  is  below  a  Strouhal  number  of  2.0.  Rakl  used  a  far-field 
microphone  at  45s  from  the  jet  center  line  and  cross-correlated  this  signal  with 
the  in- jet  sensor,  or  pseudo  sound  measurement. 

g 

Lee  and  Rlbner  correlated  an  ln-jet' hot-wire  measurement  of  the  fluctua¬ 
ting  velocity  to  a  far-field  microphone  and  their  data  indicates  that  little 
information  existed  above  Sj-2.5.  Scharton  and  Meecham^  also  measured  the 
fluctuating  pressure  in  the  jet  plume  and  found  that  the  1/3  octave  band  sound 
pressure  level  was  20  dB  down  at  Sq<~4  from  the  maximum.  Based  on  these 
previous  investigations  it  can  be  said  that  a  Strouhal  number  of  no  greater 
than  4  is  required  to  adequately  measure  a  cold  subsonic  jet.  For  a  2  inch 
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diameter  jet  with  a  1000  fps  exit  velocity  the  maximum  frequency  of  interest 
is  then  24  KHz. 

Several  probe  possibilities  exist  and  will  be  discussed  briefly  below: 

1)  Conventional  microphones,  of  the  1/8"  size,  have  an  adequate  response 

greater  than  50  KHz.  With  an  aerodynamlcally  shaped  nose  cone,  this 
* 

is  a  strong  candidate. 

11)  Siddo&'s  probe  has  a  minimized  dipole  radiation,  but  in  Its  present 
form  has  a  linear  amplitude  response  range  of  less  than  2.5  KHz  and 
an  unspecified  phase  response. 

Ill)  Long  snouted  adaptors  can  be  applied  to  conventional  microphones  to 
reduce  the  volume  of  the  probe,  at  the  sensing  point,  and  still  have 
adequate  amplitude  response.  However,  the  phase  response  suffers 
due  to  the  long  tansaisslon  path.  For  a  probe  with  a  3  Inch  trans¬ 
mission  tube  the  phase  shift  at  24  KHz  Is  approximately  2000  degrees. 


B,  Probe  Contamination  Factor 

Rakl6  and  Siddon7  have  defined  a  probe  contamination  factor  which  is  the 
ratio  between  the  mean  squared  acoustic  pressure  due  to  dipole  type  radiation 
from  the  probe  surface  to  the  mean  squared  acoustic  pressure  due  to  quadrupole 
type  radiation  from  the  adjacent  eddy  of  the  fluid.  For  cylindrical  probes 
this  factor,  C,  can  be  expressed  as 


C  = 


4420  /  d  \  4  /  dpt 

(u'  /U)^2  '  x  '  '  d  • 

rms 


where:  u’ 

rms 

M 

d 

d 

P 

X 

U 


RMS  value  of  the  fluctuating  velocity 

Local  Mach  number 

Jet  diameter 

Probe  diameter 

Downstream  distance 

Local  velocity 
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A  sample  calculation  is  given  below: 


For 


Mexit  " 
X 


/U 

rms 


d 

C 


0.8 

5d,  M  ■ 
«  0.25 
1/8" 

2" 

0.0027 


Mexit 


This  calculation  indicates  that  a  1/8"  microphone  is  probably  adequate 
for  hi.jh  subsonic  jets  and  subsonic  portions  of  supersonic  jets  with  reasonable 
exit  diameter.  The  probe  used  by  Rakl  and  Slddon  had  a  C  ■  0.1. 


C.  Repeatability 

It  is  important  that  accurate  means  exist  to  calibrate  the  probe  to  an 
absolute  scale.  In  this  respect  the  commercially  available  probes  generally 
will  be  better. 

2 . 2 . 1 . 2  Pressure  Probe  Selection 

The  1/8"  condensor  microphone  with  an  aerodynamical ly  shaped  nose  cone  and 

a  long  90°  probe  support  to  align  the  probe  with  the  flow  was  selected  as  in-jet 

6 

acoustic  probe.  It  was  felt  that  an  airfoil  probe  such  as  used  by  Rakl  and 
Siddon  7  could  not  have  been  developed  in  time,  with  an  adequate  frequency 
response.  Consultation  with  Prof.  T.E.  Slddon  revealed  that  their  present 
probe  resonated,  Indicating  an  inductive-capacitance  type  action,  and  it  was 
necessary  to  fill  the  microphone  cavity  with  cotton  to  dampen  this  oscillation. 
Therefore,  enlarging  the  hole  size  would  have  little  effect  on  the  dynamic 
response.  Slddon  felt  that  this  present  probe  is  not  the  ultimate,  but  rather, 
a  first  model  of  this  type,  and  he  is  actually  pursuing  a  refined,  more 
miniaturized  version  of  this  probe.  Using  conventional  1/8"  microphones,  it 
should  be  possible  to  develop  a  Siddon-type  probe  with  a  5-6  KHz  response,  but 
this  is  still  too  low. 


2.2.2  Probe  Selection 

Velocity  mapping  in  the  2"  diameter  cold  jet  exhaust  at  Mach  number  of 
0.8  was  performed  using  hot  film  anemometers  with  single  and  "x”  probes.  The 
special  "x"  hot  film  probe,  mounted  on  a  wedge  shape  stem,  was  purchased  from 


Therm-System  Inc.  for  two  dimensional  velocity  measurements.  Figure  117 
shows  a  schematic  picture  of  this  "x"  hot  film  probe  which  consists  of  two  hot 
films  on  the  wedge  at  45°  with  respect  to  the  jet  axis  and  perpendicular  to 
each  other  viewed  from  side.  Thus,  by  taking  the  sum  and  the  difference  of  the 
signals  from  these  two  films,  the  axial  and  the  radial  components  of  the  jet 
flow  can  be  obtained,  respectively  (Reference  9).  A  single  film  probe, 

TSI  model  1210-20,  war.  chosen  to  provide  injet  velocity  measurements  with 
minimum  disturbance  in  the  flow  for  the  velocity  and  the  in-jet  pressure  cross 
correlation  investigation. 

2.3  Experiments 

Jet  noise  experiments  were  conducted  on  an  outdoor  test  facility  with 
cold  air  exhausting  through  a  2  inch  diameter  convergent  nozzle.  Thin  nozzle 
is  fed  by  a  12  inch  plenum  equipped  with  straightening  screens  and  honey¬ 
comb,  which  is  in  turn  fed  by  a  long  4  inch  line.  The  jet  axis  is  parallel  to 
a  grass  covered,  ground  plane  at  a  height  of  5.5  feet. 

Far-field  acoustics  are  measured  with  a  1/2"  Bruel  and  Kjaer  microphone  on 
a  10  foot  arc  (in  the  same  horizontal  plane  as  the  jet  axis)  and  ranging  from 
19°  to  146°  from  the  jet  axis.  In-jet  and  farfield  acoustic  signals  are  fed  to 
a  large  indoor  Sangamo  (Sabre  IV)  tape  recorder  and  processed  through  a  Hewlett 
Packard  mini-computer.  General  Radio  1921  Real  Time  Analyser,  and  a  Salcor 
Correlator  (400  point). 

The  temperature  difference  between  the  jet  total  temperature  and  the 
ambient  temperature  ranged  from  22s  F  to  a  maximum  of  40°  F  with  the  ambient 
temperature  ranging  from  a  low  of  43®  F  to  63*  F. 

The  ambient  noise  level  for  this  facility  is  quite  low,  65  to  67dB, 

(re.  0.0002pbar),  and  the  far-field  data  ranged  from  100  to  115  dB  so  that  the 
ambient  contamination  is  very  low. 

Several  types  of  experimental  data  were  gathered  and  are  listed  below: 

Far-field  acoustics 

In-jet  unsteady  pressure  measurements 

In-jet  unsteady  velocity  measurements 

Simultaneous  measurements  of  unsteady  jet  pressure  and  far-field 
acoustics 


Two-point  unsteady  jet  pressure  measurement 

Simultaneous  measurements  of  unsteady  jet  velocity  and  far-fleld 
acoustics . 

2.3.1  In-Jet  Measurements 

2. 3. 1.1  Static  Pressure  Fluctuations 

The  fluctuating  static  pressure  within  the  plume  of  the  2"  jet,  has  been 
mapped  with  a  rozzle  pressure  ratio  of  2  (slightly  under-expanded).  Initial 
trials  at  a  pressure  ratio  of  3  showed  that  the  pressure  field  was  very  intense, 
and  except  for  very  large  downstream  distances,  greater  than  15  diameters,  the 
1/8"  B&K  microphone  could  not  tolerate  the  high  amplitudes  and  therefore,  a 
nearly  sonic  jet  was  used.  Before  taking  these  measurements  it  was  necessary 
to  determine  the  response  of  the  1/8"  B&K  microphone  and  90°  angle  support  to 
the  jet-induced  vibration.  This  was  accomplished  by  covering  the  static 
pressure  opening  of  the  microphone  nose  cone  with  silicone  rubber.  This  micro¬ 
phone  was  then  subjected  to  an  intense  acoustic  field,  which  would  have  produced 
a  rms  signal  of  100  times  the  quiescent  signal  if  the  opening  had  not  been 
sealed.  The  sealing  technique  prevented  the  acoustic  signal  as  the  output 
signal  did  not  change  in  the  presence  of  the  field.  The  micophone  was  then 
submerged  in  the  jet  and  the  vibration  induced  "rms"  signal  was  mapped. 

Comparing  this  rms  signal  to  the  signals  achieved  with  an  unblocked  nose  cone 
showed  that  the  acoustic  plus  vibration  signals  were  from  14  dB  to  20  dB  greater 
than  that  of  the  vibration  alone.  Since  this  is  a  random  signal,  the  correction 
to  obtain  a  pure  acoustic  signal  is  less  than  0.1  dB. 

Four  downstream  stations  were  measured  at  1,3,5  and  10  diameters  downstream. 
At  each  of  these  points  a  radial  traverse  of  the  overall  average  (rms)  jet 
pressure  signal  was  plotted  on  a  X-Y  recorder.  These  plots  are  shown  on  Figures 
118  through  121.  Measurements  of  this  type  have  been  made  before  by  Scharton 
and  White. ^  The  measurements  reported  here  show  more  definite  peaks  and  valleys 
at  the  shorter  downstream  distances;  this  is  to  be  expected  since  the  probe 
to  jet  diameter  ratio  was  smaller  for  the  current  measurements.  Scharton  'and 
White  show  that  maximum  "rms"  pressure  on  the  jet  center  line  occurs  at  5 
diameters  downstream  while  these  measurements  indicate  that  the  overall  jet 
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pressure  rises  uniformly  along  Che  center  line  from  one  diameter  to  10 

diameters  downstream.  The  maximum  jet  pressure  does  not  occur  on  the  center 

line  for  x/d  of  5  and  less.  The  radial  position  and  the  maximum  value  of  the 

jet  pressure  for  each  axial  station  are  shown  on  Figure  122.  The  absolute 

maximum  jet  pressure  measured  was  43  of  the  plenum  gage  pressure  or  5.73  of  the 

jet  dynamic  pressure.  This  level  agrees  well  with  the  Scharton  and  Meecham^ 

results  on  a  sonic  jet,  formed  by  the  exhaust  of  turbojet  and  with  the  low 

12  6 

velocity  results  (m  -  0.35)  of  Planchon  and  Jones  and  Rakl  .  These  results 
indicate  that  the  maximum  overall  jet  pressure  level  is  rather  insensitive  to 
jet  diameter,  temperature  and  upstream  turbulence.  The  radial  locus  of  the 
jet  pressure  maximum  lies  somewhat  inside  corresponding  measurements  of  the 
fluctuating  axial  velocity. 

The  shape  of  the  .aximum  jet  pressure  level  versus  axial  distance  Is  very 

similar  to  the  noise  source  distribution  for  a  2.8  Inch  diameter  hot  jet 

14 

(1150°  F)  given  by  MacGregor  and  Slmcox.  An  even  better  fit  could  probably 
be  obtained  by  making  a  comparison  based  on  a  cross-sectional  area  weighting 
of  the  jet  pressure  level. 

Along  this  locus  of  maxima,  a  1/3  Octave  band  analysis  was  made  of  the 
signal.  An  example  of  this  is  shown  on  Figure  123;  the  spectral  shape  is  very 
similar  for  all  axial  locations  with  the  closer  axial  stations  having  a  slightly 
more  peaked  spectrum.  The  Strouhal  number  of  the  peak  1/3  octave  band  level 
is  shown  on  the  insert  to  Figure  123. 

2 . 3 . 1 . 2  RMS  Velocity  Mapping 

Figure  117  shows  a  schematic  diagram  of  the  apparatus  for  obtaining  the 
hot  film  velocity  data.  Signals  from  both  films  of  the  "x"  probe  were  linearized 
and  calibrated  to  an  accuracy  of  +0.253  of  the  velocity  measured  between  400  to 
1000  ft/sec.  Two  dual  channel  Hewlett  Packard  differential  amplifiers  with 
unity  gain  were  used  to  obtain  the  sum  and  difference  from  the  hot  film  signals, 
and  give  the  axial  and  radial  velocity  components,  l.e.,  u'  and  v*.  The  mean 
axial  and  radial  velocities,  i.e.,  u  and  v  were  obtained  from  the  sum  and  the 
difference  of  the  average  voltage  readings  at  the  outputs  of  the  two  hot  film 
linearizers.  A  TSI  Model  1060  RMS  meter  was  employed  at  the  out-puts  of  the 
HP  differential  amplifiers  to  give  the  turbulence  levels  of  the  axial  and 
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radial  velocities.  The  time-varying  axial  and  radial  velocities  from  the  HP 
differential  amplifiers  were  recorded  on  the  Sangamo  Tape  recorder  simultan¬ 
eously  with  the  far-field  pressure  signal  from  a  microphone  located  10  ft  away 
from  the  nozzle  exit  and  at  43.8°  from  the  jet  axis.  The  mean  velocity  and  the 
turbulence  level  distributions  in  the  Mach  0.8  jet  exhaust  are  shown  in 
Figure  124.  Good  agreements  of  the  mean  axial  velocity  distribution  and 
turbulence  level  distribution  are  found  between  the  present  measurements  and 
those  presented  in  References  15  and  16.  The  radial  velocity  of  the  jet  was 
found  to  be  about  102  of  the  jet  exit  velocity  and  decaying  slightly  in  the 
downstream  direction.  The  turbulence  levels  in  the  axial  direction  were 
found  to  be  about  two  times  of  those  in  the  radial  direction  at  most  locations 
in  the  jet  exhaust. 

2.3.2  Pressure  Probe  Contamination 

Probe  contamination  of  the  far-field  signal  is  always  of  concern  in  this 
type  of  measurement  and  can  be  determined  by  comparing  the  far-field  signal 
with  and  without  the  inlet  probe.  On  Figure  125  it  can  be  seen  that  the  "rms" 
value  of  the  far-field  SPL  is  severely  affected  by  the  presence  of  the  in-jet 
probe  at  the  higher  angles.  This  result  is  due  to  the  different  radiation 
patterns  of  moving  (jet)  and  fixed  (probe  contamination,  dipole-like)  noise 
sources,  with  the  moving  sources  having  a  radiation  pattern  biased  towards 
the  shallower  angles  and  therefore,  less  contaminated.  Subsequent  results 
show  clearer  correlations  at  the  shallow  angles.  The  contamination  is  also 
seen  to  be  greater  as  the  downstream  kistance  is  increased;  part  of  this  is 
due  to  the  greater  insertion  required  to  reach  the  maximum  "rms"  jet  pressure 
level.  One-third  octave  band  analysis  of  the  far-field  signals  (with  and 
without  the  in-jet  probe)  reveals  that  the  contamination  appears  at  frequencies 
higher  than  the  peak  intensity  (Strouhal)  frequency. 

2.3.3  In-Jet  to  Far-Field  Correlations 
2. 3. 3.1  Jet  Pressure  to  Far-Field 

Cross-correlations  of  an  1/8"  B&K  microphone,  buried  in  the  plume  of  a  2" 
convergent  jet,  with  a  far-field  microphone  have  been  found  to  be  sensitive 
to  the  far-field  microphone  position  and  in-jet  axial  location.  At  the  four 
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axial  stations  of  x/d  ■  1,  3,  5  and  10  the  in-jet  microphone  was  placed  on  the 
radial  point  of  maximum  overall  (rms)  jet  pressure  level  and  then  cross-corre¬ 
lated  with  several  microphone  positions,  ranging  from  20°  to  60c  from  the  jet 
axis.  These  correlations  were  for  the  most  part  processed  using  a  400  point 
Saicor  correlator  and  checked  with  a  Hewlett  Packard  Fourier  Analyser  and 
mini-computer. 

Figure  126  shows  these  results  for  the  unfiltered  cross-correlations. 

Band  pass  filtering  of  the  signals  will  undoubtedly  improve  the  correlation, 

and  may  allow  a  discernible  correlation  at  the  shallow  angles.  In  particular, 

note  that  larger  downstream  distances  have  much  higher  correlations.  When  the 

probe  is  placed  at  the  larger  axial  distances  it  encounters  a  larger  scale  of 

turbulence  and  should  yield  a  better  correlation.  Work  not  reported  here,  but 

currently  carried  out  in  suppressor  models,  with  much  finer  turbulent  scales, 

shows  a  much  lower  level  of  correlation  with  a  similar  directivity  effect. 

18 

Similar  measurements  were  made  by  Meecham  and  Hurdle  with  in-jet  probe  at 
5.2  diameters  downstream  of  a  6.5"  jet  (M*0.99)  and  one  diameter  off-axis. 

The  correlation  levels  are  in  good  agreement  and  the  directivity  only  slightly 
poorer  agreement  with  the  results  presented  here  (see  Figure  10  of  Reference  18) . 

Siddon  and  Rakl**’7  measured  a  correlation  coefficient  of  0.02  with  the 

far-field  at  90°  and  the  injet  probe  at  x/d  =  5  and  r/d  =  1.  This  result  is 

18 

in  agreement  with  the  present  findings  and  those  of  Meecham  and  Hurdle. 

18  6 

tloise  source  location  methods  suggested  by  Meecham  and  Hurdle,  Rakl  and 
Siddon^  require  the  complete  mapping  of  the  correlation  function  for  the  entire 
jet  plume  and  then  taking  the  second  time  derivative  of  these  correlations. 

The  coi.relograms  themselves,  without  differentiation,  range  from  fairly 
definitive  to  very  marginal,  so  that  this  process  is  tenuous,  at  best. 

Figures  127  through  129  show  the  unfiltered  cross-correlation  shapes  for  one 
axial  position  and  they  illustrate  the  range  of  clarity  that  exists.  The  shape 
of  the  correlation  functions  shown  here  tends  to  more  anti-symmetric  than 
symmetric  shapes  predicted  by  others.  When  closer  to  the  jet  and  at  larger 
far-field  angles  the  correlation  shape  is  completely  obscured  by  noise.  The 
use  of  cross-correlations  of  in-jet  pressure  to  the  far  field  acoustic  pressure 
as  a  noise  source  location  method  does  not  appear  to  be  very  attractive  as  it 
requires  a  very  tedious  and  marginal  procedure.  An  improved  correlation  method 
to  determine  the  axial  source  distribution  would  be  to  place  a  small  probe 
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close  to  the  jet  flow  (to  obtain  good  definition  of  source  location)  but  not 
in  the  plume  to  eliminate  probe  contamination  of  the  far-field  and  probe 
pick-up  of  self-noise.  In  this  manner,  the  distribution  along  the  jet  axis  of 
the  near-field  radiation  to  the  far  field  could  be  determined. 

2 . 3 . 3 . 2  In-Jet  Velocity  to  Far-Fleld  Pressure 

If  r  is  the  directional  vector  connecting  the  injet  sound  source  point 
and  the  far-field  microphone  location,  the  cross  correlation  between  the  square 
of  the  injet  velocity  component,  Vr,  and  the  far-field  pressure,  P(«S) ,  can  be 
obtained  by  the  cross-correlation  of  the  axial  and  radial  velocity  components, 
i.e.  u  and  v,  with  the  far-field  pressure  in  the  following  manner  (Reference  12): 
where  i  is  the  angle  between  r  and  the  jet  axis  and  t  is  the  time  delay.  The 
advantages  of  employing  the  above  formula  are  twofold:  (1)  velocity  components 

•f 

u  and  v  are  usually  easier  to  obtain  than  Vr;  and  (2)  for  an  experimental  set 
up  with  more  than  one  far-field  microphone,  the  cross-correlation  of  the 
velocity  vs.  the  pressure  signal  from  every  microphone  can  be  calculated  with 
one  simultaneous  measurement.  The  measurement  of  the  velocity  component,  Vr, 
is  usually  made  by  using  a  single  hot  wire  probe  oriented  in  the  direction 
parallel  to  r  (Reference  4) .  This  kind  of  measurement  can  not  be  easily 
justified  as  an  absolute  Vr  measurement,  because  a  single  hot  wire  probe  can 
not  provide  one  dimensional  velocity  measurement  unless  the  flow  field  is  purely 
one  dimensional.  For  an  axisymmetric  flow  field  such  as  jet  exhaust  from  a 
round  nozzle,  the  transverse  velocity  perpendicular  to  Vr  and  the  wire  will 
always  be  measured  along  with  Vr.  On  the  other  hand,  If  a  proper  "x"  hot  wire 
probe  is  used  (Reference  9) ,  the  velocity  components  u  and  v  can  be  resolved 
and  measured  with  no  difficulty.  The  objective  of  using  the  special  "x"  hot 
film  probe  for  the  2"  Diameter  Jet  Velocity  Mapping  (Section  2.2.2)  is  to 
provide  the  u  and  v  measurements  for  the  cross-correlation  of  the  in-jet 
velocity  to  the  farfield  microphone  pressure  signal. 

Because  the  correlations  between  velocity  and  pressure  are  very  weak, 

(about  5%),  the  narrow  band  analysis  procedure  has  to  be  adopted  (Reference  13). 
Figure  130  shows  a  block  diagram  of  the  instruments  used  for  the  cross  correla¬ 
tion  study.  A  typical  narrow  band  cross  correlation  ARy-+  2p^  *8  8^l0Wn  *n 
Figure  131  where  f  is  the  filter  central  frequency.  Due  to  some  mismatch 
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between  the  two  narrow  bandpass  filters,  the  noise  associated  with  the  cross 
correlation  output  is  high  and  masks  the  correlation.  Better  in-jet  velocity 
and  farfield  pressure  cross  correlation  results  should  be  obtained  in  the  near 
future.  The  in-jet  pressure  correlations  to  the  farfield  acoustic  signal  are 
much  clearer  as  evidenced  by  the  fact  that  a  discernible  correlation  could  be 
obtained  at  these  same  conditions  without  filtering  (see  Figure  126). 

2 . 4  In-jet  Correlations 

A.  Radial  Coherence 

Simultaneous  measurements  of  fluctuating  static  pressure  In  the  jet  plume 

were  made  at  two  points  and  cross-correlated  to  determine  the  tangential  and 

radial  coherence.  Cross  correlation  of  the  jet  pressure  and  fluctuating 

velocity  are  used  to  determine  the  axial  coherence.  The  primary  goal  of  these 

19 

experiments  were  to  see  if  the  results  of  Fuchs  would  be  valid  at  higher  jet 
Mach  numbers.  Fuchs  used  a  4  inch  diameter  jet  and  1/4"  Bruel  and  Kjaer  micro¬ 
phones  at  a  jet  Mach  number  of  0.12,  while  these  experiments  used  a  2"  jet  with 
1/8"  Bruel  and  Kjaer  microphones  so  that  the  principal  difference  will  be  due 
to  the  effects  of  higher  speed.  The  high  degree  of  radial  coherence  found  by 
Fuchs  (at  zero  time  delay)  does  not  hold  at  this  high  flow  velocity  (M=l). 
Figure  132  shows  the  zero  time  delay  cross  correlation  as  well  as  that  for 
the  optimum  time  delay  with  one  microphone  on  the  jet  centerline  and  the  other 
microphone  at  various  radii.  The  structure  of  the  jet  pressure  field  (at  this 
Mach  number)  is  quite  unlike  that  of  Fuch's  results.  The  zero  time  delay  point 
is  just  one  point  on  a  cross-correlation  function  and  is  not  necessarily  the 
most  interesting.  It  does  not  yield  any  information  on  propagating  disturbances 
and  will  only  indicate  the  level  of  coherence.  The  correlations,  however,  show 
strong  indications  of  repetitive  wave-like  motion  in  the  radial  direction. 
Figure  133  shows  these  function  for  r  -  0.711";  the  regularity  of  the  peaks 
and  the  asymmetric  behavior  rear  the  origin  are  the  predominant  features.  From 
these  peaks  and  valleys  a  reflection  length  r  -  1.54"  was  determined.  Another 
point  of  interest  is  the  appearance  of  a  sign  inversion  at  the  origin  (zero 
time  dilay),  as  can  be  seen  on  Figure  134.  On  the  right  side  of  Figure  134  the 
outer  microphone  signal  is  delayed  so  that  inward  going  pressure  waves  are 
being  correlated,  and  on  the  left  side  the  centered  microphone  is  delayed  so 
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that  outward  going  waves  are  being  correlated.  For  no  sign  inversion  the 
correlation  function  should  be  symmetric  about  the  origin  and  going  to  unity 
as  r  tends  to  zero.  The  sign  inversion  seen  here  appears  also  on  the  other 
correlation  functions  for  the  larger  radial  separations.  A  source  ring  can 
also  be  postulated  since  the  time  delay  for  the  correlation  peak  closest  to 
the  origin  changes  sign  as  the  radial  separation  increases;  this  can  be  seen 
on  Figure  132.  The  location  of  this  source  ring  is  found  to  be  Vs**0.67" 
connecting  the  time  delays  of  two  radial  separation  points  with  a  Mach  line. 

This  radial  source  location  is  used  on  Figure  133  to  predict  the  location  of 
the  3econd  minimum  and  the  corresponding  repetitions.  The  agreement  of  these 
predicted  peaks  and  valleys  with  the  actual  correlation  function  is  quite  good 
and  tends  to  verify  the  model  employed  here.  A  complete  Mach  wave  diagram  was 
constructed  using  another  observation,  that  the  signals  reverse  sign  as  they 
pass  through  the  reflection  radius,  Vf.  Although  this  is  a  most  unusual 
phenomenon,  agreement  of  the  predicted  position  and  the  signs  of  the  correlation 
peaks  with  the  data  is  quite  good  particularly  since  a  very  simple  acoustic 
propagation  model  (Mach  line)  was  employed.  This  experimentally  deduced 
reflective  model  of  the  jet  agrees  in  principle  with  some  of  the  later 
theoretical  models  (Reference  20)  that  suggest  that  the  jet  flow  acts  like  a 
wave  guide  to  internal  sources.  Also  the  model  of  a  large  scale  pulsating 
cell  as  proposed  by  Scharton  and  White^  does  not  seem  to  be  supported  by  this 
evidence;  it  would  require  the  correlation  functions  to  be  symmetrical  about 
the  time  origin  and  show  peaks  and  valleys  on  both  sides .  The  peak  and  valley 
seen  here  could  be  due  to  a  harmonic  content  of  the  signal,  as  both  the 
1/3  octave  specra  of  the  jet  pressure  (see  Figure  123)  and  the  far-field 
(at  90°)  microphone  peak  near  the  4.4  KHz  frequency  that  corresponds  to  the 
peak  to  peak  wave  length,  but  auto-correlations  of  the  jet  pressure  did  not 
show  any  sign  of  harmonic  content. 

The  cross-correlation  technique  has  been  seen  to  be  a  powerful  technique 
to  define  the  nature  of  the  jet  pressure  structure  rather  than  just  to  deter¬ 
mine  the  zero  time  delay  coherence  as  has  been  the  primary  result  of  previous 
investigators . 
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B.  Circumferential  Coherence 


Circumferential  coherence  measurements  were  taken  at  two  different  radii, 

one  on  the  outside  edge  of  the  jet  (r/d  =  1.43)  and  one  buried  in  the  jet 

(r/d  =  0.36).,  For  the  outer  edge,  these  results  (Figure  135)  are  of  the  same 

19 

magnitude  as  those  of  Fuchs  and  confirm  that  the  outer  edges  of  the  jet  do 

21 

undergo  some  soit  of  toroidal  motion  as  observed  by  Crowe  and  Champagne. 

In  the  interior  of  the  jet,  however,  the  coherence  is  down  by  a  factor  of  ten 
with  very  little  coherence  except  at  A0  =  45°.  This  indicates  that  even  if 
the  outer  layers  are  exhibiting  coherent  (toroidal)  motion  the  interior  portion 
has  a  much  lower  level  of  coherence. 


C.  Axial  Coherence 

To  determine  the  axial  characteristics  of  the  jet  pressure  field  would 

require  two  microphones  to  be  placed  on  an  axial  line.  This  would  cause  a 

probe  interference  problem  on  the  downstream  microphone  and  was  not  attempted. 

A  hot  film  probe  was  used  instead  as  the  upstream  probe,  since  the  upper 

velocity  limit  on  the  hot  film  probe  is  about  M  *=  0.9,  (the  jet  Mach  number 

was  reduced  to  0.8)  to  assure  probe  survival.  The  microphone  was  fixed  at  the 

axial  station  and  the  hot  film  probe  traversed  upstream  away  from  the  micro- 

2 

phone.  Both  the  unsteady  axial  velocity  u'  and  (u')  were  corrrelated  to  the 
jet  pressure.  A  typical  correlation  function  is  shown  on  Figure  136  with  the 
corresponding  convection  velocity  shown  on  Figure  137.  The  convection 
velocity  at  Ax  =  0.23,  is  about  50%  of  the  local  mean  velocity  (at  the  micro¬ 
phone)  with  an  increase  to  80%  of  the  local  mean  velocity  at  the  maximum 
separation. 

19 

Contrary  to  Fuchs  low  speed  observations,  no  evidence  of  wave-like 
disturbances  appeared.  The  correlation  monotonically  decreases  as  the 
separation  increased  as  shown  on  Figure  138.  Fuchs'  data  was  taken  on  the 
jet  center-line  and  is  not  directly  comparable;  however,  the  wave-line 
character  appears  in  nearly  all  of  his  other  correlations.  Some  were  taken 
with  ore  microphone  on  the  jet  axis  and  the  other  slightly  off-axis  (0.5d)  and 
moved  parallel  to  it.  Apparently  this  axial  wave-like  character  does  not  exist 
at  the  higher  jet  velocities  and  is  a  low  speed  phenomenon!.  Fuchs  used  the 
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zero  time-delay  correlation.  If  that  point  were  used  here,  the  level  of 

correlation  would  be  reduced  to  about  one-tenth  of  the  optimum  time  delay 

2 

value.  By  squaring  the  linearized  hot  film  signal,  the  (u')  signal  was 
obtained  and  correlated  to  the  jet  pressure.  Figure  139  shows  that  the 
level  of  correlation  is  only  about  1/3  to  1/2  of  that  of  velocity  to  pressure 
correlation  indicating  that  the  cross  product  of  unsteady  velocity  u’  and  mean 
velocity  (shear  noise)  contributes  more  strongly  to  the  jet  static  pressure 
fluctuation  than  the  square  of  the  unsteady  velocity  (self  noise). 


APPENDIX  IV  -  1 


Effects  of  Grid  Spacing  on  LV  Spectra  Estimation 

LV  spectra  will  be  estimated  by  reconstructing  the  auto-correlation  fun¬ 
ction  by  evaluating  the  average  of  each  lag-product  value.  When  this  is  done, 
an  error  is  introduced  because  lag  time  values  are  quantized  on  a  grid.  If 
the  grid  has  spacing  T,  then  the  lag  estimate  for  T  =  NT  +  T  +  T/2.  We  must  deter¬ 
mine  how  small  to  make  T  so  that  the  induced  error  is  sufficiently  small.  This 
error  may  be  analyzed  in  the  manner  of  Balakr ishnan+:  Consider  the  auto-cor¬ 
relation  function  estimate: 

Rxx(T)  =  E  {X(T)  X  (t  +  T  (T))} 

where  t  is  the  nominal  lag  time  value 

f,  is  the  random  variable  describing  the  variation  in  the  lag 
product  of  t 

x  is  the  process  under  analysis 

If  x  and  t,  are  independent,  and  if  all  processes  are  stationary  and  ergodic, 

then 

***  (I)  'EC  <Rx*  +5(T))) 
where  F.  is  the  expectation  over  £. 

If  the  auto-correlation  function  is  sampled  on  a  grid  of  T,  then  -T/2 
T 

C  0  )l  and 

R  (t)  =  E  {T  R_  (t  +  £  (nT))  +  l  y0(x  -  nT)} 

**  C  xx  n-  — 

The  spectrum  is  then  the  Fourier  transform  of  Equation  (3),  which  is 

c  v  l  r°“c  /  %  K(nT)(w  -  y)  “  ,  27m. 

SKX  =  (u  -  b)eJ  * J,  Voiv  -  —  )dy 

We  now  consider  cases  of  the  above.  First  if  the  £(x)  are  identical  (which 

would  be  the  case  of  Poisson  sampling,  neglecting  (0)  ),  then 

C  (t)  =  £ 
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Interchanging  Integration  and  expectation,  we  find 


Sxx^  “  -J  Sxx  ^  (u-v)  u0  (y  -  ^jr)  dy 

where  (w)  *  E  [e^uJ  -  characteristic  function  of  the  r.v.£. 
Thus,  the  spectrum  becomes 


S 

xx 


(u) 


S 

XX 


(oi 


which  Is  the  aliased  version  of  the  original  spectrum  times  the  characteristic 
function.  If  we  assume  that  £  is  uniformly  distributed,  then 

,  („)  .  2  si.  m 

E  uT 


Thus,  the  error  is  zero  at  u  ■  0  (a  general  result  since  $  (0)  -  1  for  all 
characteristic  functions)  and  Increases  with  frequency.  Reasonable  values  of 
T  (i.e.,  complying  with  the  Nyquest  criterion)  limits  this  error  at  high  fre¬ 
quencies  to  about  3dB. 

This  model,  however,  neglects  some  important  effects.  First,  the  estimate 
of  (0)  will  contain  no  jitter.  Second,  near  T  ^  0,  there  will  be  different 
distributions  for  £  (t)  (i.e.  £)  will  be  function  of  t)  due  to  the  finite  par- 
tical  transit  and  processor  dead  times.  In  a  similar  manner,  we  may  show  that 


S 

xx 


(u) 


+  E  f  _Jm  [Sxx(u-p)J[^  (iD-p)  -  $g(u»-y)]  e  ^yT  dp 


where  "most"  of  the  jitters  are  Identical  with  4>  (w)  their  characteristic  func- 

8 

tlons  and  the  exception  lag  values,  t^,  have  jitters  with  characteristic  func¬ 
tions  <f>^  (w)  • 
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From  the  above  we  may  determine: 


1.  A  "noise"  error  (second  term)  is  added  to  the  spectrum  if  the  jitters  have 
different  distributions. 

2.  Since  all  4(ui)-  1  at  u  ■  0,  the  smaller  the  gird  spacing  and  frequency, 
the  smaller  the  noise  error,  since 

[$  (w)  -  4  (u))3  -►  0 

VL  8 

Several  simulations  were  conducted  under  various  assumptions  simulating 
expected  spectra  conditions.  From  these  it  was  concluded  that  if  T<Jn/5, 

where  is  the  Nyqulst  period,  the  induced  errors  were  at  least  30  dB  below 
the  peak  of  the  original  spectra.  Thus  a  grid  spacing  of 

T  1  Tn/5 

should  be  used  in  LV  spectra  estimates.  Thus,  if  40  KHz  is  assumed  as  the 
required  Nyquist  frequency,  then  a  grid  spacing  of 

T  <_  5y  sec 

is  required  in  the  analysis  system. 
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Figure  102  Laser  Velociaeter  Optics  Bsckags 
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b)  INTER  PARTICLE  ARRIVAL  HISTOGRAM 
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Figure  106  Comparison  of  Laser  Velocimeter  and  Hot  Film  Turbulence  Spectra 
M  *  0-5,  X/D  =  6,  r/r  =1,  80  Hz  Bandwidth 


a)  AUTOCORRELATION  FUNCTION  (1st  RECORD) 


M  =  1.0  TT  =  1500°R 

x/d  =  10  r/rQ  =  0 


b)  AUTOCORRELATION  FUNCTION  (2nd  RECORD) 


Figure  107  Autocorrelation  Function 
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Figure  112  Mean  Velocity  Profile  Comparison  of  LV  with  Hot  film 
X/D  -  2.0,  M,  -  0.5 
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Figure  115  Laser  Velocineter  Measured  Mean  Velocity  and  Turbulence  Intensity 
Radial  Profiles;  H,  -  1-55,  Tm  =  1500*R,  x/D  •  9-6 
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Figure  118  In-Jet  Static  Pressure  "RMS"  Profile 
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Figure  119  In-Jet  Static  Pressure 
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Figure  120  In-Jet  Static  Pressure  "RMS"  Profile 


Figure  121  In-Jet  Static  Pressure  "RMS”  Profile 
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Figure  129  Unfiltered  Cross -Correlation  of  In-Jet  Microphone  and  Far-Field  Microphone: 

2  in.  Diameter  Convergent  Jet  at  Pr  «=  2,  In-Jet  at  X  «=  20  in.  on  Centerline 
Far-Field  at  R  *  10  ft.  19»1°  from  Centerline 
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Figure  133  Radiol  Cross -Correlation  Function 
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Figure  135  Circumferential  Cross -Correlation  Unfiltered  Microphone 
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Figure  139  Axial  Cross -Correlation  Coefficient  of  Unsteady  Velocity  Squared, 
and  Jet  Pressure,  p' 
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